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1.  INTRODUCTION 


The  United  States  Coast  Guard  maintains  approximately  50,000 
Federal  aids  to  navigation  on  United  States  waterways.  The  Coast 
Guard  publishes  information  on  these  navigational  aids  in  Light 
Lists  and  Notices  to  Mariners.  For  approximately  4100  lighted 
buoys,  this  information  includes  the  flash  characteristic  of  the 
light  (Morse  A,  quick  flash,  etc.)  and  nominal  range. 

A  typical  omnidirectional  aid  to  navigation  lantern  consists  of 
fresnel  drum  lens  with  a  vertical  filament  lamp  at  its  center. 
The  drum  lens  refracts  light  toward  the  optical  plane  (plane 
^  perpendicular  to  lantern's  vertical  axis) ,  thus  increasing 

^  optical  plane  intensity  over  that  provided  by  a  bare  lamp  alone. 

Intensity  peaks  in  the  optical  plane,  and  rapidly  falls  off  to 
zero  above  and  below  the  optical  plane  with  an  approximately 
Gaussian  profile,  as  shown  in  Figure  1.  Profile  full  width  at 
half  maximum,  called  vertical  divergence,  is  about  5  degrees  for 
the  Coast  Guard  155mm  buoy  lantern. 

Currently,  nominal  range  is  the  standard  performance  measure  of 
all  navigational  lights  and  indicates  the  distance  at  which  most 
observers  will  detect  a  light.  Nominal  range  calculations  are 
based  on  optical  plane  peak  intensity.  These  calculations  are 
appropriate  for  fixed  structures  where  the  optical  plane  is 
constantly  aligned  with  the  mariner's  line  of  sight.  Each  flash 
of  the  light  produces  the  same  peak  intensity  in  the  mariner's 
direction  and  results  in  a  constant  detectability. 

However,  nominal  range  calculations  are  inappropriate  for 
floating  aids  to  navigation  where  the  optical  plane  pitches  and 
rolls  with  the  buoy.  As  the  buoy  rolls,  intensity  in  the 
mariner's  direction  continually  changes  between  its  peak  value 
and  zero,  as  illustrated  in  Figure  2.  This  variation  in 
intensity  over  time  causes  the  detection  distance  of  the  buoy 
light  to  change  over  time.  The  detection  distances  may  be 
represented  as  a  probability  distribution  function  with  the  most 
probable  detection  distance  being  significantly  less  than  the 
calculated  nominal  range  value. 

After  detecting  a  lighted  buoy,  the  mariner  must  identify  its 
characteristic.  This  requires  detecting  enough  flash  pulses  to 
distinguish  it  from  other  possible  lighted  buoys.  Buoy  motion 
may  cause  some  flashes  to  go  undetected,  making  the  recognition 
task  more  difficult. 

Solutions  to  the  problem  of  detecting  and  recognizing  flashed 
buoy  lights  are  not  apparent,  particularly  when  no  method  exists 
to  quantify  the  problem  or  evaluate  proposed  solutions.  Greater 
vertical  divergence  produces  increased  intensities  at  angles  off 
the  optical  plane,  but  decreased  intensity  on  and  near  the 
optical  plane.  Greater  flash  characteristic  duty  cycle  produces 
increased  probability  of  seeing  flashes,  but  decreased  battery 
life.  Redesigned  buoys  with  natural  roll  frequencies 


Intensity  Pattern 


Observation 

Plane 


Optical  Plane 


FIGURE  2.  Effect  of  Buoy  Angle  on  Observed  Intensity 


significantly  different  from  predominant  wave  frequencies  noy 
achieve  reduced  buoy  motion  and  increased  detection  range  o:  “ihe 
light.  A  method  to  measure  the  degree  to  which  these  may  :  ’ ve 
the  problem  has  not  been  available  until  now. 

This  study  assesses  the  impact  of  buoy  motion  on  a  marii  r's 
ability  to  detect  and  recognize  standard  Coast  Guard  buoy  li(  ts. 
It  proposes  a  new  measure  of  buoy  light  performance  ind 
introduces  a  new  method  to  evaluate  different  buoy  des  ns, 
lantern  designs,  and  flash  characteristics.  The  study  shows  uoy 
light  detectability  depends  on  vertical  divergence,  sh 
characteristic,  buoy  roll  amplitudes,  and  other  variables. 


2.  HISTORY 

Feingold,  Ritter,  and  Tozzi  (1977)  proposed  a  theoretical  :  odel 
for  predicting  probabilities  of  detecting  a  buoy  light  give;  the 
probability  density  function  of  buoy  roll  angle.  They  de  :ned 
probability  of  detection  as  the  integral  (over  all  angles)  o  the 
product  of  two  probability  density  functions:  a)  the  probab  ity 
of  detection  distribution  versus  angle  and  b)  the  buoy  roll  ngle 
probability  distribution.  The  probability  of  dete  tion 
distribution  versus  angle  was  derived  from  simple  geom  ric 
considerations . 

Hirsch  (1982)  used  the  formulation  of  Feingold  et.  al.  (197  '  to 
assess  the  influence  of  vertical  divergence  on  POD.  He  reported 
a  significant  increase  in  probability  of  detection  with 
increasing  vertical  divergence,  but  did  not  take  into  accoun  the 
necessary  decrease  in  optical  plane  intensity  due  to  spre  ling 
the  flux  over  wider  divergence  angles. 

There  have  been  other  attempts  by  Coast  Guard  engineers  to  odel 
detection  probability  of  buoy  lights,  but  none  have  been  w  lely 
accepted.  Most  efforts  and  debates  have  centered  on  ver  ical 
divergence  of  buoy  lanterns.  Flash  characteristic  is  gene  ally 
not  considered  in  probability  of  detection  models,  although  t  is 
the  simplest  and  least  expensive  design  variable  to  change  i  the 
lighted  buoy  system. 


3.  A  NEW  APPROACH 

Previous  efforts  to  predict  probabilities  of  detection  of  ro  ling 
buoys  did  not  consider  observation  time,  distance  from  buo  or 
flash  characteristic.  Intuitively,  one  expects  probabilii  .  of 
detection  to  increase  the  longer  a  mariner  looks  in  the  dire  Tion 
of  the  buoy.  One  also  expects  distance  from  the  buoy  and  Lash 
characteristic  duty  cycle  to  affect  detection  probability. 


This  author's  approach  was  to  determine  the  number  of  detected 
flashes  of  a  light  on  an  actual  rolling  buoy  for  various 
distances  and  observation  times.  Detection  and  recognition 
probabilities  were  calculated  from  a  time  series  history  of 
intensity  in  the  horizontal  plane  for  a  flashed  light  on  a 

rolling  buoy.  Probability  of  detection  (POD)  and  probability  of 
recognition  (POR)  were  modeled  as  functions  of  flash 

characteristic,  lantern  divergence,  observation  time,  observation 
distance,  and  buoy  roll. 

3.1  Definitions 

The  following  definitions  apply  to  this  report. 

Effective  Luminous  Intensity.  -  Luminous  intensity  is  the 
luminous  flux  per  unit  solid  angle,  measured  in  units  of  candela. 
The  effective  luminous  intensity  (or  equivalent  fixed  intensity) 
of  a  flashing  light  is  the  intensity  of  a  fixed-on  light  that 
produces  the  same  luminosity  as  the  flashing  light.  The  apparent 
reduction  of  intensity  for  flashing  lights  is  due  to  the  eye's 
time  response.  In  this  report  effective  intensity,  I  ,  is 
calculated  from  a  flash's  instantaneous  intensity,  I(t),  ana  peak 
intensity.  Ip,  using  the  Schmidt-Clausen  method: 

Ip  /  Ig  =  1  +  0.2  (  Ip  /y*  I(t)  dt  )  (1) 

where  the  integral  is  over  the  flash  duration. 

Nominal  Range .  -  Maximum  distance,  D,  a  light  can  be  detected 

based  on  a  light's  effective  luminous  intensity,  a  meteorological 
visibility  of  10  nautical  miles  (19  kilometers) ,  and  a  threshold 
of  illuminance  at  the  eye  of  0.67  sea-mile-candela  (0.2 
microlux) .  It  is  calculated  from  Allard's  Law: 

0.67  =  Ig  T°  /  D^  (2) 

where  the  atmospheric  transmissivity,  T,  is  determined  from  a 
meteorological  visibility,  V,  of  10  nautical  miles  by: 

T  =  0.05^/^  =  0.74  (3) 

Flash  Period.  -  The  length  of  time  to  complete  one  cycle  of  a 
flash  characteristic.  Example:  an  FL  6  (.6)  characteristic  has  a 
flash  period  of  6  seconds. 

Flash  Pulse.  -  Refers  to  a  single  period  of  time  when  the  lamp  is 
lit.  Example:  an  I  QK  FL  characteristic  has  six  flash  pulses  per 
flash  period.  A  steady  fixed  on  light  that  is  fading  on  and  off 
due  to  buoy  motion  has  only  one  flash  pulse. 

Observation  Event .  -  A  single  event  wherein  a  mariner  at  a 
specified  distance  looks  in  the  direction  of  the  buoy  for  a 
finite  but  continuous  time  period. 


Observation  Time.  -  The  period  of  time  in  seconds  a  mariner  looks 
in  direction  of  buoy. 

Probability  of  Detection  (POD) .  -  The  proportion  of  observation 
events  in  which  mariner  detects  at  least  one  flash  pulse.  More 
precisely,  it  is  the  probability  the  light  produces  an 
illuminance  of  at  least  0.67  sea-mile-candela  at  the  observer's 
position  in  any  randomly  selected  contiguous  time  interval. 
Illuminance  is  determined  from  Allard's  Law  using  effective 
intensity  (EFI) ,  distance,  and  a  visibility  of  10  nautical  miles. 

Probability  of  Recognition  (POR) .  -  The  proportion  of  observation 
events  in  which  a  mariner  detects  all  flash  pulses  within  at 
least  two  adjacent  flash  periods.  More  precisely,  it  is  the 
probability  that  in  any  randomly  selected  time  interval,  the 
light  produces  flash  pulses  for  at  least  two  contiguous  flash 
periods  which  all  have  illuminances  of  at  least  0.67  sea-mile- 
candela  at  the  observer's  position.  Illuminance  is  determined 
from  Allard's  Law  using  effective  intensity  (EFI) ,  distance,  and 
a  visibility  of  10  nautical  miles. 

90%  POD  Distance.  -  The  distance  at  which  a  mariner  has  a  90% 
probability  of  detecting  a  specified  flashing  light  on  a 
specified  rolling  buoy,  while  observing  for  a  specified  length  of 
time. 

90%  POR  Distance.  -  Same  as  90%  POD  Distance  except  mariner 
recognizes  characteristic  according  to  definition  of  POR. 

90%  POD  Nominal  Range  Factors.  -  A  light's  90%  POD  Distance 
divided  by  its  rated  nominal  range.  For  a  particular  lantern 
and  flash  characteristic,  POD  nominal  range  factor  is  a  measure 
of  range  degradation  of  a  buoy  light  due  to  buoy  motion. 

3.2  Assimptions  and  Conditions  of  Experiment 

a.  Buoy  motion  was  recorded  in  only  one  plane.  The 
mariner  is  located  in  this  plane.  Effects  of  buoy  motion  in  this 
plane  of  observation  are  representative  of  effects  of  buoy  motion 
in  all  other  planes  of  observation. 

b.  The  mariner  knows  where  to  look  for  the  buoy.  He  looks 
continuously  in  the  correct  direction  for  a  finite  period  of 
time.  Thus  if  a  flash  produces  at  least  0.67  sea-mile-candela  at 
the  observer's  position  within  the  observation  time,  it  will  be 
detected. 

c.  The  lamp  flashes  with  instantaneous  rise  and  fall 
times.  Lamp  filament  nigrescence  time  effects  are  considered 
minimal  compared  to  buoy  motion  effects.  Thus  instantaneous 
intensity  in  the  mariner's  direction  is  determined  from  buoy 
angle,  lantern  vertical  beam  profile,  and  whether  or  not  the  lamp 
is  lit.  The  ideal  square  flash  pulse  is  modified  only  by  buoy 
roll  and  vertical  divergence  of  the  lantern. 


d.  The  Schmidt-Clausen  method  (lALA  (1978))  determines 
effective  intensity  (EFI)  for  each  flash  pulse.  Flash  pulse 
duration  defines  limits  of  integration  for  this  method,  not  the 
actual  instantaneous  intensity  in  the  horizontal  plane. 

e.  If  a  mariner  detects  a  flash  pulse  at  a  certain 
distance  in  a  given  observation  time,  he  will  detect  the  same 
pulse  in  greater  observation  time  at  that  distance.  He  may  also 
detect  other  flashes  in  the  greater  observation  time. 

f.  If  a  mariner  detects  a  flash  pulse  at  a  given  distance 
within  a  certain  observation  time,  he  will  also  detect  the  same 
pulse  at  shorter  distances  within  the  same  observation  time.  He 
may  also  detect  other  flash  pulses  at  these  shorter  distances. 

g.  In  order  to  recognize  a  flash  characteristic,  the 
mariner  must  detect  all  pulses  within  two  contiguous  flash 
periods.  A  detected  flash  pulse  is  considered  a  recognized  flash 
pulse:  apparent  pulse  width  or  fading  in  and  out  of  a  single 
pulse  does  not  affect  recognition.  It  is  assumed  a  mariner  can 
distinguish  fading  of  a  pulse  due  to  buoy  roll  from  electronic 
switching  of  a  pulse.  (Note  that  this  is  similar  to 
distinguishing  rotating  beacons  from  electronically  flashed 
lights. ) 

h.  Visibility  is  assumed  constant  at  10  nautical  miles 
(0.74  atmospheric  transmissivity). 


4.  DATA  ACQUISITION 

Buoy  angle  was  sampled  every  0.2  seconds  and  recorded  using  an 
optical  buoy  motion  recording  system,  shown  in  Figure  3.  Light 
from  the  buoy  was  polarized  in  a  fixed  direction  with  respect  to 
the  buoy  axis.  As  the  buoy  tilted,  the  plane  of  polarization 
rotated.  The  buoy  motion  recording  system  detected  polarization 
rotation  in  a  plane  perpendicular  to  line  of  observation  from  the 
recording  system.  Software  was  developed  (Appendix  B 
RECORD_BUOY_MOTION)  to  control  the  recording  system  and  collect 
data. 

A  brief  hardware  description  is  as  follows;  The  buoy  motion 
recorder  uses  a  rotating  Polaroid  and  phase  sensitive  detection 
to  determine  buoy  rotation  (or  polarization  rotation)  angle.  A 
Compumotor  stepper  motor  (25,000  steps  per  revolution)  rotates 
the  Polaroid  at  about  960  RPM.  The  rotating  polaroid  is  located 
between  a  600mm  F/4.5  Canon  lens  and  an  EMI  9558  photomultiplier 
tube.  A  signal  from  an  optical  shaft  encoder  on  the  stepper 
motor  locates  a  reference  position  of  the  polaroid.  The  shaft 
encoder  signal,  after  passing  through  a  divide-by-N  counter, 
provides  two  pulses  per  revolution  (32  Hz)  of  the  polaroid. 
Polarized  light  incident  on  the  lens  results  in  a  sinusoidal 
signal  from  the  photomultiplier  tube  with  a  DC  offset.  The  DC 
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FIGURE  3.  Optical  Buoy  Motion  Recording  System 
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offset  IS  due  to  unpolarized  light.  An  Edmac  bandpass  filter 
removes  noise  and  DC  level  from  detector  signal.  An  HP  5334 
universal  counter  determines  phase  between  shaft  encoder  and 
photomultiplier  signals.  Phase  is  sampled  by  the  computer  every 
0.2  seconds  and  converted  to  instantaneous  buoy  angle. 


Laboratory  measurements  showed  the  system  to  have  a  detection 


threshold  of  0.73  sea-mile-candela.  Field  tests  showed 
atmospheric  scintillation  and  polarized  moonlight  reflected  off 
the  water  significantly  reduced  sensitivity.  Atmospheric 
scintillation  was  partially  overcome  by  extending  the  source  so 
that  light  arrived  at  the  lens  over  multiple  paths.  Data 
collection  was  avoided  during  moonlit  nights. 


Two  buoys,  an  8x26  and  5x11  (see  Figure  4) ,  were  moored  off  Avery 
Point,  Groton,  CT  to  provide  actual  buoy  motion.  The  8x26  buoy 
was  installed  with  a  battery  powered  polarized  light  source.  A 
30-inch  fluorescent  lamp  worked  best  for  overcoming  atmospheric 
scintillation.  A  2-foot  cylinder  of  retroreflective  sheeting 
was  wrapped  with  Polaroid  sheeting  and  installed  on  the  5x11  buoy 
(see  Figure  5) .  Retroreflective  material  served  to  reflect  light 
from  a  1000-watt  xenon  arc  searchlight  to  the  detector. 


Five  30-minute  continuous  recordings  of  buoy  roll  angle  versus 
time  were  obtained  on  two  buoys.  More  recordings,  spanning  a 
larger  range  of  conditions,  were  planned;  time,  atmospheric 
conditions,  and  hardware  failures  limited  data  collection.  These 
recordings  are  listed  in  TABLE  1. 


TABLE  1 


BUOY  DATA  RECORDINGS 


RMS  ROLL 


WAVE  HT  rest 


12/4/86 

11/13/86 

11/13/86 

2/13/87 

2/13/87 


8x26 


8x26 


8X26 


5x11 


5X11 


3 . 3  deg 
8 . 0  deg 
9.9  deg 
8.2  deg 

8.4  deg 


20-25  kts 


20-25  kts 


5  -10  kts 


1  -  2  ft 


10-20  kts 


1  -  2  ft 


No  effort  was  made  to  quantify  the  frequency  distribution  of  sea 
states.  Long  term  daily  monitoring  would  be  required  to  collect 
this  information.  It  is  the  author's  opinion  that  conditions 
experienced  during  the  period  of  data  collection  are  typical  of 
exposed  and  semi-exposed  locations. 
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A)  5  X  11  BUOY  B)  8  X  26  BUOY 

WEIGHT  3004  LBS  WEIGHT  11,382  LBS 


FIGURE  4.  Subject  Buoy  Dimensions  and  Weights 
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FIGURE  5.  5  X  11  Buoy  with  Polarizing  Retroreflector 
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5.  DATA  ANALYSIS 


5.1  Probability  Calculations 


Software  was  developed  to  calculate  POD  and  POR  as  a  function  of 
observation  time,  observation  distance,  flash  characteristic, 
lens  vertical  divergence,  and  buoy  motion  (Appendix  B 
BUOY_PROBABILITIES)  .  The  program  required  two  data  file  inputs, 
a  recording  of  buoy  angle  versus  time  and  vertical  intensity 
profile  data  for  the  lantern.  Vertical  intensity  profile  was 
provided  in  0.2-degree  increments.  It  was  used  to  convert  buoy 
roll  angle  to  intensity  in  the  horizontal  plane.  The  ten 
standard  flash  characteristics  listed  in  USCG  (1983)  were  used  in 
the  analysis.  These  characteristics  are  defined  in  Figure  6 
(taken  from  USCG  (1983)). 


Table  2  lists  lanterns,  with  respective  file  names,  used  to 
determine  detection  and  recognition  probabilities.  File  names 
are  given  for  convenience  in  identifying  the  data.  The  155mm  and 
200mm  lantern  profiles  were  measured  in  the  photometric 
laboratory  at  the  U.S.  Coast  Guard  Research  and  Development 
Center.  The  last  five  lantern  profiles  are  simulated  profiles. 
These  gaussian  profiles  provide  the  same  total  flux  as  the  yellow 
155mm  lantern.  Software  was  written  to  generate  these  profiles 
(Appendix  B  -  LanternProfileMaker)  to  study  the  effect  of 
vertical  divergence  on  detection  and  recognition  probabilities. 


TABLE  2 


LANTERN  DATA  FILES 


Buoy  angle  versus  time  was  converted  to  intensity  as  a  function 
of  time  in  the  horizontal  plane  for  a  fixed-on  light.  This  was 
done  by  mapping  buoy  roll  angles  into  the  appropriate  vertical 
divergence  file.  Prior  to  mapping,  additional  buoy  angle  data 
points  were  interpolated  to  increase  intensity  resolution. 
Interpolation  gave  an  effective  intensity  sampling  rate  of  35  Hz. 

A  nonaalized  flash  characteristic  was  then  superimposed  onto  the 
intensity  versus  time  function.  Five  hundred  30-second  windows 
were  randomly  selected  from  the  fixed-on  light  intensity  versus 
time  record  and  multiplied  by  a  normalized  flash  characteristic. 
Phase  between  selected  windows  and  normalized  flash 
characteristic  was  also  randomized.  Software  testing  with 
generated  sinusoidal  and  cosinusoidal  buoy  roll  data  assured 
proper  randomization  in  the  algorithm. 

Effective  intensities  were  calculated  (using  the  Schmidt-Clausen 
method)  for  each  flash  pulse  in  the  30-second  window.  In  each 
window,  time  periods  (corresponding  to  observation  times)  varying 
from  2  to  30  seconds  were  searched  for  minimum  and  maximum 
intensities.  Minimum  intensities  determined  POR  and  maximum 
intensities  determined  POD  for  various  combinations  of 
observation  distance  and  observation  time. 

Figure  7  illustrates  the  procedure  for  converting  buoy  roll 
versus  time  to  flashed  intensity  versus  time  for  one  30-second 
data  window.  The  top  graph  is  a  typical  30-second  recording  of 
buoy  angle  versus  time.  The  middle  graph  is  this  recording 
converted  to  intensity  versus  time  for  a  fixed-on  light.  The 
bottom  graph  is  the  product  of  the  middle  graph  and  a  normalized 
FL  6  (0.6)  flash  characteristic.  For  this  particular  window,  the 
observer  (depending  on  distance)  might  detect  the  first  two  flash 
pulses  but  not  the  others.  This  window  would  allow  detection  for 
all  observation  times  from  2  seconds  to  30  seconds.  Had  these 
two  pulses  occurred  at  24  and  30  seconds,  detection  would  not 
occur  until  24  seconds  of  observation  time  had  elapsed. 

Matrices  of  detection  and  recognition  probabilities  as  functions 
of  observation  time  and  distance  were  calculated  for  each  of  ten 
standard  flash  characteristics  and  for  various  combinations  of 
buoy  roll  data  and  lantern  types.  Results  presented  here  are 
extracted  from  a  collection  of  280  such  matrices,  far  too  many  to 
be  reproduced  in  their  entirety  .  Figures  8  and  9  are  examples 
of  these  matrices.  Figure  8  shows  POD  of  a  perfectly  vertical 
light  (zero  degree  roll  amplitude).  The  light,  a  155  mm  lantern, 
has  a  peak  intensity  of  232  candela.  With  a  0.6-second  contact 
closure  time,  this  light  produces  a  nominal  range  of  6  nautical 
miles  according  to  USCG  (1984).  Indeed,  Figure  8  shows  this  is  a 


The  280  data  matrices  can  be  obtained  by  writing  to  Commanding 
Officer,  USCG  Research  and  Development  Center,  Avery  Point, 
Groton,  CT  06340-6096,  Attention:  Chief,  Physics  Branch. 
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FIGURE  8.  Typical  POD  Matrix  For  Zero  Roll 
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FIGURE  9.  Typical  POD  Matrix  For  Non-zero  Roll 


6-mile  light.  Figure  9  shows  the  probability  of  detecting  he 
same  light  with  only  10  degrees  rms  buoy  roll.  It  is  low 
something  less  than  a  6-mile  light. 

5.2  Frequency  Spectra  of  Buoy  Roll 

Since  detection  and  recognition  probabilities  depend  on  buoy 
roll,  it  is  desirable  to  characterize  the  rolling  motion  of  a 
buoy.  Spectral  analysis  is  one  tool  for  characterizing  motion. 
A  frequency  spectra  or  power  spectral  density  (PSD)  gives  the 
roll  amplitude  (or  energy)  at  various  frequencies.  The 
usefulness  of  this  tool  is  demonstrated  by  the  fact  that  although 
waves  and  buoys  move  semi-randomly ,  their  power  spectral 
densities  are  constant  over  long  periods  of  time. 

Additionally,  if  a  buoy  is  assumed  to  be  a  linear  system,  its  PSD 
may  be  determined  from  wave  PSD  by  a  simple  multiplication  of  two 
frequency  curves.  If  sea  spectra  were  flat  (white  noise) ,  buoy 
PSD  peaks  would  correspond  exactly  to  a  buoy's  natural  roll 
frequency.  But  sea  spectra  are  not  flat.  They  contain 
predominant  frequencies  which  produce  actual  roll  frequencies 
slightly  different  from  natural  roll  frequencies.  Fortunately, 
researchers  have  measured  wave  PSD  in  waterways  all  over  the 
globe.  Thus  one  should  be  able  to  predict  the  roll  spectra  of  a 
buoy  for  practically  anywhere  in  the  world,  given  the  mechanical 
response  of  the  buoy  to  waves  (called  the  transfer  function)  is 
known . 

This  is  analogous  to  predicting  the  output  of  an  electrical 
network  or  filter  for  random  noise  input.  The  time  interval  over 
which  the  output  signal  of  a  filter  is  correlated  to  a  single 
frequency  approximately  equals  the  reciprocal  of  its  full 
bandpass.  Likewise  for  a  buoy,  the  time  interval  over  which  a 
buoy  rolls  at  a  constant  frequency  is  indicated  by  the  width  of 
roll  PSD  peaks.  This  correlation  time  was  discovered  to  affect 
POD  and  POR,  thus  the  importance  of  buoy  roll  spectral  analysis. 
These  results  are  discussed  in  section  6.5. 

Efforts  to  measure  buoy  transfer  functions  directly  and  ver.fy 
response  linearity  were  thwarted  by  failure  of  the  wave-rider 
buoy.  It  was  not  discovered  until  late  in  the  experiment  that 
marginal  operation  of  the  wave-rider  buoy  produced  invalid  wave 
height  data.  However,  responses  to  step  inputs  on  both  buoys 
were  measured.  This  step  input  is  used  to  determine  the  natural 
roll  frequency  of  the  buoy,  as  well  as  the  buoy  transfer 
functions  and  frequency  response.  Appendix  A  discusses  the 
mathematics  of  spectral  analysis,  the  linearity  assumption,  and 
derivation  of  transfer  functions  from  step  responses^ . 


Software  for  making  these  calculations  can  be  found  in  Appendix 
B  -  FastFourierTransform. 


A  typical  step  response  recording  for  the  5x11  buoy  is  shown  in 
Figure  10.  Figures  11  and  12  show  relative  magnitudes  of 
transfer  functions  derived  from  5x11  and  8x26  step  responses 
according  to  Appendix  A.  Figure  11  shows  the  natural  roll 
frequency  of  the  5x11  buoy  to  be  0.273  Hz,  giving  a  roll  period 
of  3.7  seconds.  Figure  12  shows  the  natural  roll  frequency  of 
the  8x26  buoy  to  be  0.156  Hz,  giving  a  roll  period  of  6.4 
seconds.  The  8x26  buoy  has  a  narrower  bandpass  of  frequencies, 
meaning  less  damping  and  longer  correlation  times. 


6.  RESULTS 

6.1  Observation  Time  and  Dist2mce  Effects 

As  expected  POD  and  POR  are  strongly  affected  by  observation 
distance  and  length  of  time  spent  looking  in  the  direction  of  the 
buoy.  Figures  13  and  14  show  the  relationship  of  POD  and  POR  to 
observation  time  at  a  position  of  4  nautical  miles  from  an  8x26 
buoy  with  10  degrees  rms  roll.  POD  approaches  maximum  faster  for 
high  flash  rate  characteristics  than  for  low  flash  rate 
characteristics . 

An  unexpected  result  shows  POR  is  higher  for  FL  2.5  than  for  QK 
FL.  This  result  is  found  consistently  with  other  rms  amplitudes 
and  with  the  5x11  buoy.  This  is  due  to  correlation  between  flash 
period  and  buoy  roll  period  and  is  discussed  in  section  6.5 

Figures  15  and  16  show  the  relationship  of  POD  and  POR  to 
observation  distance  at  a  constant  12-second  observation  time. 
These  data  are  from  the  8x26  buoy  with  10  degrees  rms  roll.  As 
expected  POD  drops  with  distance  from  the  buoy  and  drops  more 
slowly  for  the  higher  flash  rate  characteristics.  Again  POD  is 
higher  for  the  QK  FL,  but  POR  is  higher  for  the  FL  2.5. 

6.2  Performance  Measurement 

Variation  of  POD  and  POR  with  observation  time  and  distance  makes 
it  difficult  to  measure,  specify,  or  compare  lantern  performance. 
Comparing  POD  and  POR  performance  of  different  buoys,  lanterns 
and  flash  characteristics  is  easier  if  the  probability  matrices 
are  reduced  to  a  single  number,  similar  to  currently  used  nominal 
range.  This  was  done  by  searching  the  matrices  for  the  distance 
at  which  90%  POD  occurs  with  a  12-second  observation  time.  Since 
higher  PODs  occur  at  shorter  distances  (or  longer  observation 
times)  than  this  90%  POD  distance,  it  is  a  reasonable  measure  of 
system  performance.  A  90%  POD  level  and  12-second  observation 
time  were  arbitrarily  chosen.  (See  section  3.1  definitions  for 
90%  POD  and  POR  distances.) 

Ninety  percent  (90%)  POD  and  POR  distances  for  12-second 
observation  were  determined  from  the  probability  matrices  and  are 
given  as  Tables  3  and  4.  These  data  are  for  a  yellow  I55mm 
lantern  with  a  1.15-ampere  lamp.  Distances  are  given  in  nautical 
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miles  as  a  function  of  buoy  roll  and  flash  characteristic. 
Column  labels  give  buoy  type  /  rms  roll  amplitude  (degrees) .  Row 
labels  give  flash  characteristic.  Ninety  percent  (90%)  FOR 
distances  are  not  available  for  the  I  QK  FL  and  MO  (A) 
characteristics  since  two  periods  of  these  characteristics 
exceeds  the  12 -second  observation  time. 

It  is  apparent  that  probabilities  of  detecting  and  recognizing 
flashing  lights  on  rolling  buoys  are  a  strong  function  of  flash 
characteristic  and  buoy  motion.  Characteristics  with  low  duty 
cycles  and  low  repetition  rates  have  lowest  probabilities  of 
detection  and  recognition.  Characteristics  with  high  repetition 
rates  or  high  duty  cycles  have  highest  probabilities  of  detection 
and  recognition. 


TABLE  3 

90%  POD  DISTANCES  FOR  12 -SECOND  OBSERVATION 
VERSUS  BUOY  NOTION 


CHAR. 

•k 

N.R. 

8x26/3° 

8x26/8° 

8x26/10° 

5x11/8.2° 

5x11/8.4 

FL  6 (.6) 

6.3 

2.8 

1.8 

1.8 

2.1 

2.1 

FL  4 (.4) 

6.1 

5.1 

2.2 

2.2 

2.1 

2.1 

FL  2.5(.3) 

5.9 

5.2 

2.3 

2.2 

2.8 

3.0 

QK  FL 

5.9 

5.5 

5.1 

5.1 

5.1 

5.3 

GP  FL  5 

6.1 

5.5 

3.0 

3.0 

3.8 

3.8 

GP  FL  6 

6.4 

5.9 

5.3 

5.1 

5.3 

5.1 

I  QK  FL 

5.9 

5.5 

4.5 

4.5 

4.4 

4.6 

MO  (A) 

6.5 

6.0 

5.2 

5.2 

5.5 

5.5 

E  INT  6 

6.6 

6.1 

5.7 

5.6 

5.8 

5.7 

OCC  4 

6.6 

6.2 

5.9 

5.8 

6.0 

5.8 

*N.R.  =  Nominal  Range  determined  for  lamp/ lantern  combination 

measured  in  photometric  laboratory. 


Ninety  percent  (90%)  POD  distances  for  each  light  may  be  compared 
with  its  nominal  range  rating  by  calculating  90%  POD  Nominal 
Range  Factors,  the  ratio  of  90%  POD  distance  to  nominal  range. 
Nominal  Range  Factors  for  the  yellow  155mm  lantern  are  provided 
by  Table  5.  Comparing  flash  characteristics  of  identical  duty 
cycle  (e.g.  FL  6  versus  FL  4,  QK  FL  versus  MO  (A))  indicates  that 
repetition  rate  is  more  important  than  duty  cycle  for  preventing 
nominal  range  degradation.  As  would  be  expected,  QK  FL  produces 
the  highest  nominal  range  factors  and  FL  6  produces  the  lowest. 


It  is  clear  that  nominal  range  is  not  a  good  measure  of  buoy 
lantern  performance.  Nominal  range  may  be  more  than  300%  in 
error.  For  example,  an  FL6  (0.6)  characteristic  provides  a 
greater  nominal  range  (due  to  its  longer  contact  closure  time) 
than  the  QK  FL  characteristic  (US  Coast  Guard;  1984) ,  but  results 
reported  here  showed  that  QK  FL  performed  significantly  better 
than  FL6  (0.6).  In  some  cases  QK  FL  produced  more  the  3  times 
the  90%  POD  distance  as  FL6  (0.6). 


TABLE  4 

90%  FOR  DISTANCES  FOR  12-SECOND  OBSERVATION 
VERSUS  BUOY  NOTION 


FLASH  CHAR. 
FL  6 (.6) 

FL  4 (.4) 

FL  2. 5 (.3) 
QK  FL 
GP  FL  5 
GP  FL  6 
I  QK  FL 
MO  (A) 

E  INT  6 
OCC  4 


Higher  POD  does  not  always  imply  higher  POR,  as  shown  by 
comparing  the  two  group  flash  characteristics  (GP  FLs)  with  FL 
2.5  in  Tables  3  and  4.  Both  GP  FL  5  and  GP  FL  6  have  higher  POD 
distances  than  FL  2.5  because  of  longer  flash  durations,  but 
lower  POR  distances  due  to  longer  flash  periods.  A  similar 
anomaly  occurs  with  QK  FL  and  FL  2.5  due  to  correlation  between 
flash  period  and  roll  period  (section  6.5). 

6.3  Vertical  Divergence  Effects 

As  vertical  divergence  increases,  optical  plane  intensity  (peak 
intensity)  decreases,  but  the  chance  that  a  light  will  be  visible 
through  the  duration  of  a  pulse  increases.  Effects  of  increasing 
vertical  divergence  were  investigated  by  calculating  detection 
and  recognition  probabilities  for  lanterns  with  various  vertical 
divergence  profiles.  Vertical  divergence  is  represented  by  the 
full  vertical  width  of  the  beam  at  half  peak  intensity.  Profiles 


were  simulated  by  distributing  light  in  a  Gaussian  form  with  full 
divergence  widths  of  4.2,  5,  10,  15,  and  20  degrees  (see  Table 
2)  .  All  simulated  profiles  contained  the  same  total  light  flux, 
that  of  a  yellow  155mm  /  1.15A  lantern.  Calculations  show  the 
effects  of  changing  the  lens,  but  maintaining  the  same  lamp  and 
power  consumption. 

Buoy  roll  data  were  taken  from  the  8x26  buoy  with  10  degrees  rms 
roll.  Ninety  percent  (90%)  POD  distances  were  determined  from 
probability  of  detection  matrices.  Table  6  gives  90%  POD 
distances  for  these  Gaussian  profiles.  90%  POD  distances  for  the 
yellow  155mm  are  also  repeated  for  convenience  under  "LANT”. 
These  data  are  plotted  in  Figure  17  for  four  different 
characteristics . 

Increasing  vertical  divergence  increases  detection  distance  for 
flash  characteristics  with  low  repetition  rates  and  low  duty 
cycles,  but  slightly  degrades  detection  distance  for 
characteristics  with  high  repetition  rates  or  high  duty  cycles. 
The  FL  6  characteristic  benefits  most  from  increasing  vertical 
divergence.  The  OCC  4  characteristic  benefits  least  from 
increasing  vertical  divergence. 


TABLE  5 

90%  POD  NOMINAL  RANGE  FACTORS  FOR 
12-SECOND  OBSERVATION  VERSUS  BUOY  MOTION 


FLASH  CHAR. 

8x26/3° 

8x26/8° 

8x26/10° 

5x11/8.2° 

5x11/8.4 

FL  6 (.6) 

.44 

.29 

.29 

.33 

.33 

FL  4 (.4) 

.84 

.36 

.36 

.34 

.34 

FL  2. 5  (.3) 

.88 

.39 

.37 

.47 

.51 

QK  FL 

.93 

.86 

.86 

.86 

.90 

GP  FL  5 

.90 

.49 

.49 

.62 

.62 

GP  FL  6 

.92 

.83 

.80 

.83 

.80 

I  QK  FL 

.93 

.76 

.76 

.75 

.78 

MO  (A) 

.92 

.80 

.80 

.85 

.85 

E  INT  6 

.92 

.86 

.85 

.88 

.86 

OCC  4 

.94 

.89 

.88 

.91 

.88 

A  figure  of  merit  for  each  profile  can  be  obtained  by  calculating 
90%  POD  distances  weighted  by  frequency  of  occurrence  of  each 
characteristic  (Table  6) .  A  random  sampling  of  62  lighted  buoys 
from  a  current  Coast  Guard  Light  List  gave  the  following 


Vertical  Divergence  (degrees) 


FIGURE  17. 


90%  POD  Distance  Versus  Vertical  Divergence 
For  12-Second  Observation  Time 


distribution  of  flash  characteristics: 


FL  6  (.6) 

4 

6% 

FL  4  (.4) 

43 

69% 

FL  2.5  (.3) 

5 

8% 

QK  FL 

5 

8% 

MO  (A) 

5 

8% 

Applying  these  weights  gave  weighted  average  90%  POD  distances 
listed  in  Table  6.  These  averages  predict  that  a  lens  with  10  to 
15  degrees  vertical  divergence  is  optimum.  This  lens  would  yield 
a  73%  average  increase  in  90%  POD  distance.  In  the  case  of  the 
most  commonly  used  characteristic,  FL  4  (.4),  doubling  the 
divergence  nearly  doubles  the  detection  range. 


TABLE  6 

90%  POD  DISTANCES  FOR  12-SECOND  OBSERVATION 
VERSUS  LANTERN  DIVERGENCE 


FLASH  CHAR. 

LANT 

GAU4 

GAU5 

GAUIO 

GAU15 

GAU20 

FL  6( .6) 

1.8 

0.4 

0.8 

3.0 

3.9 

3.5 

FL  4 (.4) 

2.2 

1.8 

2.5 

3.7 

3.7 

3.5 

FL  2.5(.3) 

2.2 

1.8 

2.3 

3.6 

3.6 

3.5 

QK  FL 

5.1 

5.1 

5.0 

4.5 

4.0 

3.5 

GP  FL  5 

3.0 

3.5 

4.0 

4.1 

4.0 

3.5 

GP  FL  6 

5.1 

5.1 

5.1 

4.6 

4.1 

4.0 

I  QK  FL 

4.5 

4.5 

4.5 

4.2 

4.0 

3.5 

MO  (A) 

5.2 

5.1 

5.1 

4.6 

4.2 

4.0 

E  INT  6 

5.6 

5.6 

5.5 

5.0 

4.5 

4.0 

OCC  4 

5.8 

5.8 

5.6 

5.0 

4.5 

4.0 

Wt .  Average 

2.7 

2.2 

2.8 

3.8 

3.8 

3.5 

In  order  to 

resolve  the 

long 

standing 

200mm  versus 

155mm 

lantern 

debate  among  Coast  Guard  engineers,  POD  calculations  were  made 
using  vertical  profiles  of  these  two  lanterns.  The  profiles  were 
measured  in  the  USCG  R&D  Center  photometric  laboratory  using 
identical  lamps.  Weighted  averages  were  calculated  as  above  and 
are  given  in  Table  7.  Some  characteristics  are  slightly  better 


on  the  200inin  lantern  and  others  are  better  on  the  ISSmin  lantern. 
Results  show  the  200min  lantern  provides  only  3%  more  average 
range  than  the  155mm  lantern.  Clearly,  there  is  very  little 
difference  between  the  two  lanterns. 


TABLE  7 


COMPARISON  OF  1551BB  AND  200lBm  CLEAR  LANTERNS 
(90%  POD  DISTANCES) 


FLASH  CHAR. 


200mm/l. 15A 


2.4 

2.9 

2.7 


155mm/l. 15A 


1.9 

2.6 

2.5 


6.4  Buoy  Roll  Amplitude  and  Frequency  Effects 

Strong  dependence  of  POD  and  POR  on  roll  amplitude  is  immediate  ly 
apparent  from  the  data  presented  thus  far.  As  Table  5  shovs, 
significant  nominal  range  degradation  results  with  less  than  lo 
degrees  rms  roll  using  standard  aids-to-navigation  hardware.  in 
general,  as  rms  roll  amplitude  increases,  POD  and  POR  decrease, 
but  nonlinearly.  Although  insufficient  data  exist  to  determine 
this  exact  relationship.  Table  5  indicates  nominal  range 
degradation  may  approach  a  limiting  value  beyond  about  8  degrees 
rms  roll. 

Buoy  roll  power  spectral  densities  were  determined  directly  from 
angle  versus  time  recordings.  Roll  PSDs  are  given  in  Figures  I8 
through  22.  Prior  to  applying  the  Fourier  transform,  time  domain 
data  were  adjusted  in  three  ways:  a)  any  DC  offset  was  removed  or 
minimized,  b)  high  irequency  noise  was  removed  by  convolution, 
and  c)  leakage  of  high  frequencies  into  the  low  frequency  range 
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was  minimized  by  multiplying  with  either  Hamming,  Hanning,  or 
Hanning  squared  windows  (Ramirez  (1985)). 

Figures  18  through  22  show  power  spectral  densities  (PSDs)  of  the 
five  buoy  roll  recordings.  Ordinate. values  on  PSD  plots  are  in 
units  of  degrees  squared  per  Hertz"^.  Rms  roll  amplitude  is 
proportional  to  the  square  root  of  the  area  under  the  PSD  curve 
(Appendix  A) .  Thus  roll  amplitudes  are  reduced  if  bandwidth  and 
peak  response  of  the  buoy  transfer  function  are  reduced.  This  is 
analogous  to  electrical  bandpass  filters.  Reducing  bandwidth 
reduces  the  amount  of  noise  (wave)  power  transferred  to  system 
output . 

PSDs  show  strong  components  near  natural  roll  frequencies. 
Figure  18  shows  a  peak  roll  component  at  0.18  Hz  giving  a  roll 
period  of  5.6  seconds.  This  conforms  closely  to  Milne,  Hsu,  and 
Smith  (1975)  who  show  roll  periods  for  the  8x26  buoy  of  about  5.5 
seconds. 

Roll  amplitudes  are  reduced  if  buoys  are  designed  with  natural 
roll  frequencies  different  from  frequencies  of  sea  spectra  peaks. 
A  cursory  review  of  the  literature  shows  that  both  the  5x11  and 
8x26  buoys  have  natural  roll  periods  near  sea  spectra  peaks. 
Colburn  (1974)  showed  PSD  sea  spectra  peaks  at  wave  periods  from 

5.5  to  15.5  seconds  off  the  California  coast.  Yamazaki  and 
Herbich  (1985)  show  spectral  density  peaks  at  wave  periods  from 

3.5  to  7.6  seconds  in  the  Gulf  of  Mexico  near  Galveston,  Tx. 
Schwab,  Bennett,  and  Liu  (1984)  observed  wave  periods  between  2 
to  6  seconds  in  Lake  Erie. 

6.5  Roll  Period  /  Flash  Period  Correlation  Effects 

Roll  amplitude  effects  were  expected  and  are  easily  understood. 
Less  apparent  is  the  dependence  of  POD  and  POR  on  actual 
frequency  components  present  in  buoy  motion.  Compare  the  results 
of  the  8x26  with  8  degrees  roll  and  two  5x11  buoy  recordings  in 
Tables  3  and  4.  The  two  5x11  buoy  recordings  both  have  rms  roll 
amplitudes  greater  than  8  degrees,  yet  they  provide  greater 
average  detection  and  recognition  distances  than  the  8x26  buoy 
with  8-degree  rms  roll.  This  indicates  POD  and  POR  depend  on 
buoy  roll  spectra  frequencies. 

Figures  18  through  22  show  differences  in  frequency  composition 
between  the  two  buoys.  Bandwidths  indicate  buoy  roll  is 

^Actual  frequency  resolution  equals  the  sampling  rate  divided  by 
the  number  of  points  transformed.  Sampling  rate  and  number  of 
points  were  5  Hz  and  8192  points  for  three  8x26  buoy  recordings 
and  one  5x11  buoy,  giving  a  frequency  resolution  of  0.00061  Hz. 
Sampling  rate  and  number  of  points  for  the  remaining  5x11  buoy 
recording  were  2  Hz  and  4096,  giving  a  frequency  resolution  of 
0.00049  Hz. 


correlated  over  a  time  interval  of  about  16  seconds  for  the  8x26 
buoy  and  about  10  seconds  for  the  5x11  buoy.  Buoy  motion  is  not 
completely  random  and  correlation  with  the  flash  characteristic 
will  affect  detection  and  recognition  probabilities. 

To  illustrate  this  consider  the  ideal  case  of  a  buoy  rolling  with 
a  single  frequency  f  and  constant  amplitude.  Two  possible 
outcomes  are:  (1)  the  "light  flashes  when  buoy  is  vertical  and  is 
detected,  and  (2)  the  light  flashes  when  buoy  is  at  maximum  roll 
angle  and  is  not  detected.  A  continuum  of  outcomes  lie  between 
these  two,  where  detection  may  or  may  not  occur,  depending  on 
lens  vertical  divergence  and  distance.  If  the  light  flashes  with 
frequency  f-  then  outcome  1  (simultaneous  flash  and  vertical 
buoy)  occurs^  every 

1 

T  =  - -  (4) 


seconds.  In  this  special  case,  if  observation  time  is  greater 
than  T  as  given  in  equation  (4),  detection  probability  is  100%. 
Detection  probability  decreases  as  T  becomes  greater  than 
observation  time.  The  roll  period  for  the  8x26  with  10  degrees 
rms  roll  is  about  5.7  seconds.  Thus  we  expect  a  decrease  in  POD 
as  flash  periods  approach  5.7  seconds. 


To  test  this,  90%  POD  distances  (for  12-second  observation  time) 
were  determined  for  five  additional  nonstandard  flash 
characteristics  on  the  8x26  buoy  with  10-degree  rms  roll.  The 
characteristics  had  flash  periods  from  4.5  seconds  to  7  seconds, 
but  all  had  flash  pulse  durations  of  0.6  seconds  to  make  a  direct 
comparison  with  the  FL  6  characteristic.  Results,  tabulated  in 
Table  8,  show  a  minimum  90%  POD  distance  when  the  flash 
characteristic  period  matches  the  buoy  roll  period. 

Similar  results  occur  with  the  5x11  buoy  which  has  a  roll  period 
of  about  4.6  seconds.  It  produced  higher  PODs  for  the  FL6 
characteristic  than  for  the  FL4  for  observation  times  greater 
than  12  seconds. 

Equation  (4)  holds  as  long  as  buoy  roll  is  correlated  to  a 
constant  frequency.  As  buoy  roll  becomes  more  random 
(autocorrelation  time  decreases)  detection  probability  increases. 
This  accounts  for  the  slightly  higher  average  POD  distance  for 
the  5x11  buoy  than  on  the  8x26  buoy  with  8  degrees  rms  roll 
(Table  3) . 

On  the  other  hand,  correlation  between  buoy  roll  and  flash 
characteristic  improves  recognition  probability.  Recognition 
requires  (according  to  definition)  detection  of  at  least  two 
consecutive  flashes.  A  rolling  8x26  buoy  crosses  vertical 
approximately  every  2.9  seconds.  A  QK  FL  characteristic  flashing 
every  second  is  less  likely  to  have  two  consecutive  flashes  occur 
at  vertical  than  a  FL  2.5  or  FL  5.5.  We  should  expect  probabili- 


TABLE  8 


90%  POD  and  POR  Distances 
Versus  Flash  Period  Correlation 


Flash  Charac.  12-sec/90%  POD  Distance  14-sec/90%  POR  Distance 


FL 

7 

(0.6) 

1.8  n. miles 

1.06 

n. miles 

FL 

6.5 

(0.6) 

1.9 

1.11 

FL 

6 

(0.6) 

1.8 

1.11 

FL 

5.5 

(0.6) 

1.5  * 

1.18 

* 

FL 

5 

(0.6) 

1.9 

1.14 

FL 

4.5 

(0.6) 

2.6 

1.33 

FL 

2.5 

(0.3) 

2.2 

2.00 

* 

QK 

FL 

5.1 

1.42 

*  Maximum  or  minimum 


ty  of  recognition  to  peak  as  flash  period  approaches  an  integer 
multiple  of  one-half  the  buoy  roll  period.  This  is  confirmed  by 
the  90%  POR  distances  (for  14-second  obseirvation  time)  listed  in 
Table  8. 


7.  DISCUSSION 

7.1  Validity  of  Assumptions 

Analysis  of  POR  assumes  recognition  occurs  with  two  complete  and 
consecutive  flash  periods.  Recognition  is  more  complex  than  this 
definition  of  POR  implies.  Color,  familiarity  with  the  waterway, 
and  intelligent  elimination  of  possibilities  also  affect 
recognizing  a  flash  characteristic.  This  definition  is  applied 
uniformly  to  all  characteristics,  although  it  allows  confusion 
between  the  QK  FL  and  I  QK  FL  characteristics.  Requiring  seven 
consecutive  periods  of  QK  FL  would  distinguish  it  from  the  I  QK 
FL,  and  would  also  significantly  lower  its  POR  values.  On  the 
other  hand,  if  I  QK  FL  and  QK  FL  characteristics  were  never 
placed  together  in  a  channel,  the  mariner  would  not  have  to 
distinguish  between  the  two. 

The  definition  of  POR  was  chosen  to  model  a  task  more  difficult 
than  detection,  but  one  which  was  simply  calculated  by  a  computer 
in  a  reasonable  amount  of  time.  POR  values  may  not  accurately 
model  recognition,  since  it  is  such  a  complex  process,  but  they 
do  place  a  lower  limit  on  detection. 


Analyses  assumed  the  mariner  knows  where  to  look  for  the  buoy. 
Thus  every  flash  pulse  which  produces  at  least  threshold 
illuminance  at  the  observer's  position  is  assumed  detected.  If 
the  mariner  does  not  know  where  to  look  for  the  buoy,  PODs  and 
PORs  would  be  lower. 

Ignoring  lamp  nigrescence  time  in  the  calculations  was  also 
driven  by  computation  time  and  algorithm  complexity.  Figure  7 
indicates  that  highest  rotational  velocities  may  produce 
intensity  rise  and  fall  times  on  the  order  of  lamp  nigrescence 
time.  Much  of  the  intensity  modulation  is  slower  than  this 
however.  Including  nigrescence  time  in  the  calculations  would 
result  in  slightly  lower  probabilities.  Larger  reductions  would 
occur  with  shorter  pulse  duration  characteristics  such  as  the  QK 
FL  and  FL  2.5. 

The  possibility  of  false  flashing  within  a  single  pulse  due  to 
buoy  rolling  was  ignored.  Pulse  duration  was  used  to  define  the 
limits  of  integration.  This  assumption  results  in  slightly 
higher  EFI  values  than  those  using  the  Blondel-Rey-Douglas 
method.  However,  both  methods  are  still  only  approximations  to  a 
psychophysical  process.  Using  the  Blondel-Rey-Douglas  method 
would  significantly  increase  an  already  very  long  computational 
process. 

No  assumption  is  made  or  implied  that  the  data  collected 
represent  average  conditions  on  U.S.  waterways.  The  sampling  is 
far  too  limited  to  represent  even  local  conditions.  However, 
Milne,  Hsu,  and  Smith  (1975)  show  that  10-degree  rolls  on  8x26 
buoys  are  not  uncommon.  They  recorded  up  to  3 2 -degree  rolls  on 
an  8x26  in  a  sea  state  with  wave  heights  of  1.2  feet  rms. 

7.2  Scope  of  Research 

The  purpose  of  this  study  was  to  determine  if  a  POD  problem 
existed  due  to  buoy  motion  and  to  outline  a  method  to  quantify 
the  problem.  Although  a  greater  number  of  buoy  roll  recordings 
were  expected,  the  goals  of  the  study  have  been  accomplished 
without  them. 

This  study  does  not  determine  (nor  was  it  intended  to  determine) 
which  is  the  best  flash  characteristic,  best  lens  divergence,  and 
best  buoy  on  U.S.  waterways.  These  questions  would  be  answered 
by  a  parametric  study,  which  this  was  not.  This  study  only 
provides  the  tools  and  outlines  the  methodology  to  perform  a 
parametric  study  of  lighted  buoy  systems. 


8.  CONCLDSIONS 


Nominal  range  is  not  a  realistic  measure  of  detection  range  of  a 
lighted  buoy.  Flash  characteristic,  sea  state,  lantern 
divergence,  and  observation  time  influence  detection  and 
recognition  range  of  a  floating  aid  to  navigation.  For  the 
relatively  sheltered  conditions  studied  here,  some  startling 
facts  were  discovered.  Among  them  are: 

1.  Detection  ranges  (90%  Probability  of  Detection)  are  as 
little  as  30%  of  the  published  Nominal  Range. 

2.  The  most  commonly  used  characteristic,  FL  4(.4),  on  an 
8x26  buoy  has  an  effective  range  of  only  1/3  of  the 
published  nominal  range. 

3.  Doubling  lantern  vertical  divergence  will  nearly  double 
effective  range  of  floating  aids  to  navigation  with  10%  duty 
cycles. 

4.  Not  only  does  amplitude  of  buoy  roll  degrade  signal 
effectiveness,  but  so  does  correlation  between  period  of 
buoy  roll  and  flash  rhythm. 

The  following  general  observations  and  trends  are  derived  from 
the  results: 

1.  Significant  nominal  range  degradation  results  with  less 
than  10  degrees  rms  buoy  roll  (section  6.4) 

2.  Increasing  lens  vertical  divergence  can  significantly 
increase  POD  of  buoy  lights  (section  6.3) 

3.  Changing  a  flash  characteristic  can  significantly 
increase  POD  of  a  buoy  light  (sections  6.2,  6.5) 

4.  POD  and  POR  decrease  with  increasing  observation  distance 
(section  6.1) 

5.  POD  and  POR  increase  with  increasing  observation  time 
(section  6.1) 

A  unique  method  of  determining  probabilities  of  detecting  and 
recognizing  flashing  lights  on  rolling  buoys  has  been 
demonstrated.  The  method  takes  into  account  observer  distance, 
observation  time,  lens  vertical  divergence,  flash  characteristic, 
and  buoy  roll  spectra.  All  five  variables  affect  POD  and  POR. 


9.0  RECOMMENDATIONS 


A  better  measure  of  buoy  light  performance  than  nominal  range 
should  be  investigated.  The  90%  POD  distance  (for  12-second 
observation  time)  is  proposed  and  used  in  this  study.  Before 
this  measure  can  be  adopted,  further  research  on  probabilities  of 
detecting  and  recognizing  flashing  lights  on  rolling  buoys  is 
needed.  Recommendations  for  further  research  are  as  follows. 

Investigate  90%  POD  distance  for  use  as  a  measure  of  buoy 
light  performance.  Ninety  percent  (90%)  POD  allows  engineers 
to  optimize  system  parameters.  The  90%  POD  distance  (or  any 
other  criterion  POD  distance)  could  be  used  without  specific 
knowledge  of  local  wave  conditions  by  adopting  standard  or 
average  sea  spectra  for  the  United  States  or  several 
geographical  regions. 

Collect  additional  data  to  extend  the  results  of  this  study 
to  cover  a  more  diverse  range  of  sea  conditions.  The 
assumption  that  a  single  axis  of  measurement  is  adequate 
should  be  tested  by  measuring  roll  angles  in  two  mutually 
perpendicular  directions.  Video  recording  may  prove  the 
best  method  for  obtaining  these  new  data. 

Examine  the  effect  of  increased  vertical  divergence  on  POD. 
Explore  methods  of  accomplishing  this  without  designing  a 
new  lens  (i.e.  frosted  lamps,  longer  lamp  filaments,  new 
lamps) . 

Formulate  a  method  to  predict  POD  and  POR  as  a  function  of 
five  variables:  observation  time,  observation  distance,  lens 
divergence,  flash  characteristic,  and  roll  spectra. 

Investigate  resistance  of  various  buoy  designs  to  motion  in 
a  variety  of  sea  states. 
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APPENDIX  A 


SPECTRAL  ANALYSIS  OF  BUOY  ROLL 


I.  LINEARITY  AND  TRANSFER  FUNCTION 

If  a  buoy  is  assumed  to  be  a  linear  system  then  in  the  time 
domain  the  buoy's  output  response  y(t)  equals  the  convolution  of 
the  buoy's  impulse  response  h(t)  and  a  wave  input  function  x(t) . 
Convolution  in  the  time  domain  is  equivalent  to  multiplication  in 
the  frequency  domain.  Thus  a  buoy's  roll  spectra  Y(f)  equals  the 
product  of  the  complex  transfer  function  H(f)  and  the  sea  spectra 
X(f)  .  The  square  of  the  absolute  value  of  H(f)  is  frequently 
called  the  response  amplitude  operator  (RAO) ,  (although  some  also 
define  |H(f) |  as  the  RAO  e.g.  Price,  Milne,  et.  al.).  The  power 
spectral  density  (PSD)  of  buoy  roll  is  the  square  of  the  absolute 
value  of  the  roll  spectra.  Defining  F{  h(t)  }  to  mean  the 
Fourier  transform  of  h(t)  these  relationships  are  given  by: 


H(f) 

=  F(  h(t)  ) 

(Al) 

Y(f) 

=  F(  y(t)  } 

(A2) 

X(f) 

=  F(  X(t)  } 

(A3) 

Y(f) 

=  X(f)  H(f) 

(A4) 

PSD 

=  |Y(f)|^ 

(AS) 

Linearity  requires  a  buoy  to  satisfy  two  criteria:  1)  any 
frequency  component  of  buoy  roll  must  be  a  linear  function  of 
wave  amplitude  of  the  same  frequency,  and  2)  the  buoy's  response 
to  any  individual  wave  frequency  component  is  independent  of  it's 
response  to  any  other  wave  frequency  component.  Price  (1976) 
demonstrates  this  is  practically  true  for  at  least  some  Coast 
Guard  buoys.  Bhattacharyya  (1978)  states  that  linearity  fcr 
floating  vessel  motion  is  a  valid  assumption  if  the  amplitudes 
are  moderate.  Although  responses  are  actually  nonlinear, 
nonlinearities  can  be  ignored  in  practice. 

It  is  difficult  to  provide  an  impulse  input  to  a  buoy  and  measure 
its  impulse  response  directly.  However  a  step  input  can  be 
simulated  by  pulling  the  buoy  over  and  quickly  releasing  it.  The 
resulting  time  series  response  may  be  normalized  by  the  initial 
starting  angle  to  simulate  a  unit  step  response.  This  response 
is  actually  a  negative  step  response  as  the  input  force  (or  buoy 
angle)  is  unity  prior  to  time  zero  and  zero  afterwards.  The 
relationship  between  the  impulse  response,  h(t),  and  a  negative 
step  response  may  be  determined  by  starting  with  the  convolution 
integral . 


{ 

\ 


If  the  input  is  a  negative  unit  step  function: 


x(t)  = 


(A'  ) 


then  the  output  y(t)  is  the  negative  unit  step  response.  In  tie 
convolution  integral  above,  the  input  function  is  reversed.  This 
x(T-t)  is  a  positive  unit  step  beginning  at  T: 


x(T-t)  = 


1  if  t  >  T 
0  if  t  <  T 


(Af ) 


Equation  A8  implies  that  integral  A6  is  nonzero  only  for  tine 
greater  than  T.  Thi^  the  convolution  integral  reduces  to 

y(T)  =  y  h(t)  dt  =  H(oo)  -  H(T)  =  0  -  H(T)  (AS  ) 

T 

Taking  the  derivative  of  equation  A9  with  respect  to  time  shows 
the  negative  unit  step  response  is  related  to  the  impulse 
response  by: 

d  y(t) 

-  =  -  h(t)  (AlC) 

dt 

The  derivative  theorem  for  Fourier  transforms  states  that  tl.e 
transform  of  the  derivative  of  y(t)  equals  the  product  of  Y(:) 
and  the  imaginary  radian  frequency.  Multiplication  by  f requem y 
enhances  higher  frequencies  and  attenuates  lower  frequencies . 
Multiplication  by  the  square  root  of  negative  one  implies  a  9{  - 
degree  phase  shift.  Thus  the  transfer  function  H(f)  may  le 
determined  from  the  Fourier  transform  of  the  negative  step 
response  Y(f)  according  to: 

H(f)  =  -  j  27r  f  Y(f)  (All) 


/-.A 


II.  POWER  SPECTFRAL  DENSITY  OF  BUOY  ROLL 


According  to  equations  A1  through  A5 ,  power  spectral  densit  y 
(PSD)  curves  can  be  predicted  from  transfer  functions  and  s*  a 
spectra.  The  usefulness  of  PSD  can  be  seen  from  Parseval  s 
theorem  which  states  that  the  energy  in  a  waveform  x(t)  in  the 
time  domain  equals  the  energy  of  its  transform  X(f)  in  the 
frequency  domain.  For  discrete  or  sampled  functions  Parseval 's 
theorem  is  given  by: 

N-1  N-1 

^  [  x(k)  2  t  X(n)  ]^  /  N  (Ai;  ) 

k=0  n=0 


A- 2 


if  the  discrete  Fourier  transform  is  defined  as: 


X(n) 


N-1 

^xck, 

k=0 


exp(-j2Tr  nk/N) 


(A13) 


and  its  inverse  as: 

N-1 

x(k)  =  n"^  X(n)  exp(  jZ^rnk/N)  (A14) 

n=0 

The  integral  of  the  PSD  function  from  negative  cutoff  frequency 
to  positive  cutoff  frequency  equals  the  variance  (mean  square)  of 
the  roll  amplitude.  (Cutoff  frequency  is  one  half  sampling  rate. 
Only  positive  frequencies  were  plotted  since  the  negative  half  is 
a  mirror  image  of  the  positive  half  of  the  PSD.)  Thus  PSDs 
indicate  the  energy  content  of  various  frequency  components  of 
roll.  The  range  or  bandwidth  of  predominant  frequencies  gives  an 
estimate  of  output  frequency  correlation  time. 

Research  showed  that  roll  PSD  is  one  of  the  variables  affecting 
detection  and  recognition  probabilities,  due  to  correlation 
effects  with  the  flash  characteristic.  PSD  can  be  reduced  to  a 
single  number  via  the  cross  correlation  coefficient.  Additional 
research  is  needed  to  determine  the  relationship  between  the 
correlation  coefficient  and  POD  /  POR.  Since  wave  PSDs  have  been 
measured  for  numerous  oceans  and  waterways  throughout  the  world, 
spectral  analysis  may  prove  to  be  a  useful  tool  for  determining 
optimal  flash  characteristics  and  optimal  buoy  designs. 


[  BLANK  ] 


Program  RECORD_BUOY_MOTION 


PROGRAM  RECORD_BUOY_MOTION  (INPUT.  OUTPUT)! 

{Records  buoy  roll  angle  as  a  function  tine  and  stores  data  on  silo  mass 
storage  device.  Plots  angle  in  real  tine  on  CRT.  by  Dan  Broun,  APR  87> 
IMPORT 

0GL_LIB,  {get  graphics  routines) 

HPIB_0. 

GENERAL_0 . 

GENERAL.! , 

GENERAL.Z  . 

SYSGLOBALS. 

SYSDEUS  . 

lOCOMASM. 

lODECLARATIONS) 

CONST 

HPIB  “  7i  {  hpib  interface  select  code  > 

GPIO  “  12)  {  gpio  interface  select  code  > 

HP5334A  *  704;  {  hp  univ.  counter  isc  > 

K48S  “  722s  {  auto,  piconeter  ics  > 

MIDNIGHT  •  8G40000S  {  centi  seconds  per  day  > 

ESC  »  CHR<27);  {constants  for  silo  connands) 

CrtAddr  *  3s  {  address  of  internal  crt  > 

ControlUord  "Os  {  device  controls  0  for  crt  > 

OPEN.LENS  ■  3s  {  i/o  controls  open  lens,  notor  on  ) 

CLOSE.LENS  ■4s  {  i/o  control:  close  lens,  notor  off  > 

MAXANGLE  ■  45s  {nax  degrees  shown  on  CRT  plot  } 

WAUECAL  ■  3379600s  {wave  rider  calibration,  ft/anp  > 

TYPE 

MSG.TYPE  -  STRING! 255 Is 

DIR.TYPE  ■  <  RIGHT,  LEFT  )s  {enunerated  direction  type) 

UAR 

STARTTIME  :  INTEGERS  {  daily  start  tine  for  data  collection) 

STOPTIME  :  INTEGERS  {  daily  stop  tine  for  data  collection) 

NEXTONTIME  :  INTEGERS  {  next  start  tine  for  data  collection  ) 

ONDURATION  :  INTEGERS  {  tine  stepper  notor  is  on  ) 

OFFOURATION  :  INTEGERS  {  rest  tine  between  blocks  of  data  ) 

FINALSTOPDATE, 

NEXTSTARTDATE  , 

NEXTSTOPDATE, 

TODAYSDATE  :  OATERECs 

SAMPLERATE  :  REALs  {  rate  buoy  angle  is  noasured  ) 

SIGNALFREO  :  REALs  {  signal  frequency  fron  PMT  detector  ) 

OFFSET  :  REALs  {  uc_5334a  counter  tine  when  buoy  vertical 

I  :  INTEGERS 

ANS  :  CHARs 

TIME  :  TIMERECs 

PTR  :  TEXTS 

F  :  TEXTS 

BUOYTYPE  :  STRING  CtSIs 


. . ***** . . 

<  ♦  *  ) 

( *  PROCEDURE  GET.OATE  *  ) 

(  *  •  ) 

(*  Displays  pronpt  ,  reads  date  input  and  checks  *) 

<♦  fornat  of  Input.  Returns  date  to  SET.OATES.  *> 


(  ••**«««*»*****«***«*«««*«««««*««*«**«**»***4**4««*«***  ) 

PROCEDURE  GET_DATE  ( UAR  DATE  :  DATERECi  MSG  :  MSG_TYPE>i 


UAR 

J 

INTEGER; 

{ 

loop 

counter 

> 

DATEUALID 

BOOLEAN; 

< 

date 

status 

> 

DSTR 

STRINGE40]; 

DAY 

1  .  .31 ; 

MONTH 

1  .  .12; 

YEAR 

86.. 88;  -(PROGRAM  WILL  BOMB 

IN  891 } 

DELIMIT 

CHAR; 

BEGIN 

REPEAT 

PAGE; 

DATEUALID  := 

TRUE; 

FOR  J  :=  1  TO  10  DO  WRITELNs 
WRITELN  (MSG)i 


WRITELN; 

WRITE(’The  current  date  is:  '); 

WRITELN<OATE.MONTH:2,V'  .DATE .  DAY :  2  ,  ’  /  ’  .DATE .  YEAR :  2  )  ; 
WRITELN; 

WRITELN( ’Enter  new  date  in  format:  mn/dd/yy’); 

REAOLN(  DSTR  ); 

IF  STRLEN(  DSTR  )  >  0  THEN 
BEGIN 
TRY 

STRREAD<  DSTR  .1.1  .MONTH  .DELIMIT  .DAY  .DELIMIT  .YEAR  ) 
DATE. MONTH  :»  MONTH; 

DATE. DAY  DAY; 

DATE. YEAR  :=  YEAR; 

RECOUER 

BEGIN 

PAGE; 

WRITELN<  ' Unrecognized  date  format.  Try  again.’) 
WRITELN< ’Example:  07/30/86’); 

DATEUALID  :-  FALSE; 

END; 

END; 

UNTIL  DATEUALIDi 

END;  (*  PROCEDURE  GET_DATE  *) 

{  44444  4444444444444444  4444  44444444444444444444  444444444  ) 


(4  4  ) 

<  *  PROCEDURE  SET_DATES  ♦  ) 

<  *  *  ) 

<*  Sets  dates  (todaysdate.  startdate.  stopdate).  *) 

(*  Promts  for  date  input.  Displays  current  value.  *) 

(*  I  f  OK  then  <'Roturn>  accepts  current  value  for  the  *) 
(*  date.  Each  date  is  obtained  via  procedure  ♦) 

<*  GET_DATE.  *) 

(4  4  ) 

(*  units  called  :  SYSDATE  (  in  SYSDEUS  )  ♦) 

(*  SETSYSDATE  (  in  SYSDEUS  )  *) 

(  *  Procedure  GET_DATE  *  ) 

<4  4  ) 


PROCEDURE  SET_DATES(  UAR  TODAYSDATE ,  STARTDATE .  STOPDATE 
UAR 

MESSAGE  :  STRINGCSGli 
BEGIN 

MESSAGE  ’ENTER  TODAYS  DATE’? 

GET_DATE(  TODAYSDATE .  MESSAGE  )? 

SETSYSDATE  < TODAYSDATE  )  i  {  set  system  date  > 

MESSAGE  ’ENTER  STARTING  DATE’; 

GET_DATE(  STARTDATE,  MESSAGE  ); 

MESSAGE  :=  ’ENTER  PROGRAM  STOP  DATE’; 

GET_DATE(  STOPDATE.  MESSAGE  ); 

END;  (*  PROCEDURE  SET_DATE  *) 


DATEREC  ); 


( »•«««*«***»«**«««*««***•«•«•««*«•«*•••««*«•««»*«*««**•  ) 

(  *  *  ) 

(*  PROCEDURE  CONUERT_TIMES  *) 

(  *  «  ) 

(*  Converts  time  in  centiseconds  to  a  time  record  if  •) 
direction  is  "right".  If  direction  is  left  then  *) 

(*  converts  timerec  to  centiseconds.  *) 

(  *  *  ) 

(  **«*■****•»*■•♦*♦***********#*■»********«*##*****♦****♦***  ) 

PROCEDURE  CONUERT_TIMES(  UAR  TIMEC5EC  :  INTEGER; 

DIRECTION  :  DIR_TYPE; 

UAR  TIMERECORD  :  TIMEREC  ); 

UAR  CSEC  :  INTEGER; 

BEGIN 

CASE  DIRECTION  OF 

RIGHT  :  BEGIN  <  convert  centiseconds  to  timerec  > 

CSEC  TIMECSEC; 

TIMERECORD. HOUR  CSEC  DIU  3B0000; 

CSEC  CSEC  MOD  360000; 

TIMERECORD. MINUTE  CSEC  DIU  6000; 

TIMERECORD. CENT  I  SECOND  CSEC  MOD  6000; 

END; 

LEFT  :  BEGIN  Iconvert  timerec  to  centiseconds  > 

TIMECSEC  :=  3G0000*TIMERECORD . HOUR  +  6000*TIMERECORD . MINUTE  + 
T I MERECORO . CENT  I  SECOND ; 

END; 

END;  <CASE> 

END;  <  CONUERT  TIMES  > 


(  ***** . . . . . * . *****  ) 

<  *  *  ) 

(*  PROCEDURE  GET_TIME  *) 

(  *  *  ) 

<♦  Displays  prompt  and  takes  input  for  time.  ♦) 

(*  If  <CR>  then  accepts  displayed  time  for  input.  *) 
(  *  •  ) 


( *«**«•«*«***«**«***•«•**«*«**«•**••••*««*«««*****••**•  ) 

PROCEDURE  GET_TIME  (UAR  CENTISECONDS  :  INTEGER; 

MSG  :  MSG_TYPE); 

UAR 

J  :  INTEGER;  <  loop  counter  > 

TIMEUALID  :  BOOLEAN; 

TSTR  :  STRINGCBI;  <  input  buffer  > 


{  separator 


> 


DELIMIT  :  CHAR; 

TIME  :  TIMERECj 

BEGIN 
REPEAT 

TIMEVALID  :=  TRUE? 


PAGE; 

-  FOR  J  :=  1  TO  10  DO  WRITELN;  <  set  up  pronpt  display  and  read  screen  > 
WRITELN  (MSG);  {  write  what  tine  is  for  > 

WRITELN; 

WRITE  (’It  IS  now  set  to:  ’); 

CONUERT_TIMES(  CENTISECONDS .RIGHT .TIME  ); 

WRITELN  <TIME.H0UR:2  ’  .TIME.MINUTE:2  ’  .TIME.CENTISECOND  DIU  100:2); 
WRITELN; 

WRITE  ( ’Enter  the  new  tine  in  the  forn:  hh:nn:3s  ’  ); 

READLN  (TSTR); 

IF  STRLEN  (TSTR)  >  0  THEN 
BEGIN 
TRY 

STRREAD( TSTR  .  1  .J .TIME. HOUR  .DELIMIT  .TIME. MINUTE  .DELIMIT  . 

TIME.CENTISECOND  ); 

TIME.CENTISECOND  :=  TIME.CENTISECOND  *  100; 

CONUERT_TIMES<  CENTISECONDS  .LEFT  .TIME  ); 

RECOUER 

BEGIN 

PAGE; 

WRITELN  ( ’Urecognized  tine  fornat.  Try  again.’); 

WRITELN( ’For  exanple:  12:34:Be’); 

TIMEUALID  :-  FALSE; 

END; 

END; 

UNTIL  TIMEUALIO? 

END; 


(  *  ♦  ) 

(  *  PROCEDURE  SET_TIMES  *  ) 

(  ♦  *  ) 

(*  Sets  the  various  tines  used  by  progran.  Pronpts  *) 
(*  for  input  by  user.  Tines  are  all  in  centiseconds  ♦) 
(  *  *  ) 

(•  units  called  :  SETSYSTIME  (  in  SYSOEUS  )  *) 

<*  Procedure  GET_TIME  *) 

(  *  ♦  ) 


( **♦**♦***♦***#**«****♦««•*#****«**«*#*«***♦*♦♦****♦***  > 

PROCEDURE  SET_TIMES(  UAR  STARTTIME. 

STOPTIME  . 

ONDURATION. 

OFFDURATION  :  INTEGER  ); 


UAR 

I 

MESSAGE 

CURRENTTIME 

TIME 

BEGIN 


INTEGER;  <  LOOP  COUNTER  > 

MSG_TYPEi 

INTEGER; 

TIMEREC; 


MESSAGE  :=  ’ENTER  CURRENT  TIME  OF  DAY’; 
CURRENTTIME  :=  SYSCLOCK; 


GET_TIME  (  CURRENTTIME.  MESSAGE); 

CONUERT_TIMES(  CURRENTTIME  .RIGHT  .TIME  ); 

SETSYSTIME(TIME  ); 

MESSAGE  :=  ’ENTER  DAILY  START  TIME  FOR  DATA  COLLECTION  PHASE’; 
GET_TIME  (STARTTIME.  MESSAGE); 

MESSAGE  ’ENTER  DAILY  STOP  TIME  FOR  DATA  COLLECTION  PHASE’; 
GET_TIME  (STOPTIME.  MESSAGE); 

MESSAGE  :«  ’ENTER  TIME  LENGTH  OF  A  DATA  BLOCK’; 

GET_TIME  (ONDURATION.  MESSAGE); 

MESSAGE  :=  ’ENTER  TIME  LENGTH  BETWEEN  BLOCKS  OF  DATA  (REST  TIME)’ 
GET_TIME  (OFFDURATION.  MESSAGE); 

END;  (*  PROCEDURE  SET_TIMES  *) 

( «****««****«•*«*•**«**•****«••**••«•«*«*»***««•*«•«***  ) 


(  ♦  *  ) 

(  *  PROCEDURE  GET_NUMBER  *  ) 
(  *  *  ) 
(*  Seta  up  diaplay  and  prompta  operator  for  input  *) 
(*  of  a  real  number,  Passea  input  back  to  *) 
(*  SET_DATARATES.  *) 
(  *  ♦  ) 


( **•*•«***»««*«*•*«*«««**«*«««•*«•«•««*«««»«««•««•**«•*  ) 

PROCEDURE  GET_NUMBER  ( UAR  NUM  :  REAL;  MESSAGE  :  MSG_TYPE); 

UAR 


NSTR 

:  STRINGCBI;  <  input  buffer  > 

I 

:  INTEGER;  <  counter  > 

UALIDNUMBER 

:  BOOLEAN; 

BEGIN 

REPEAT 

PAGE; 

UALIDNUMBER 

TRUE; 

FOR  I:=l  TO  10  DO  WRITELN; 

WRITELN(  MESSAGE  ); 

WRITELN(’It  ia  currently  aet  to  ’.NUM:G:4); 

READLN(  NSTR  ); 

IF  STRLEN  (NSTR)  >  0  THEN 
BEGIN 
TRY 

STRREAD  (NSTR.I  ,I  .NUM); 

RECOVER 

BEGIN 

PAGE; 

WRITELN  ('ERROR  :  Invalid  data  entry.  Try  again.’); 
VALIDNUMBER  FALSE; 

END; 

END; 

UNTIL  UALIDNUMBER; 

END;  (*  PROCEDURE  GET_NUMBER  *) 

(  »*♦*****♦**#*•»*•»********•»*#***###***#****«******♦♦♦♦**  ) 


(  ♦  *  ) 

(  *  PROCEDURE  SET_DATARATE5  *  ) 
(  #  ♦  ) 
(*  Prompta  operator  for  buoy  angle  aampling  rate,  *) 
(•  counter  time  offset,  and  PMT  signal  frequency.  *) 


(*  Prompts  are  repeated  until  valid  response.  *) 
(*  <CR>  will  take  the  current  value  being  displayed  *) 
(  *  *  ) 
<*  units  called  :  SET_NUMBER  *) 
(  *  *  ) 
{•  limits  :  off  set  constant  >=  0  ♦) 
(  *  frequency  constant  >0  ♦  ) 
(  *  *  ) 


( *****•*«********«**«*«•»****«*««•«•*«•*******«•«*••«««  ) 

PROCEDURE  SET_DATARATES  ( UAR  SAMPLERATE  :  REAL: 

UAR  OFFSET,  SIGNALFREQ  :  REAL  ); 

UAR 

I  :  INTEGER:  {  loop  counter  } 

NUM  :  REAL: 

MESSAGE  :  MSG_TYPE; 

CRLF  :  STRINGI2]; 

BEGIN 

CRLF  :=  CHR(  t  3  )+CHR(  I  0  ) ;  {carnage  return  line  feed> 

MESSAGE  :=  'INPUT  THE  BUOY  ANGLE  SAMPLING  RATE  IN  HERZ’+CRLF+ 

'3  Hz  or  leas  ->  both  buoy  and  uaverider  signals  stored ' +CRLF+ 
'greater  than  3  Hz  to  8  Hz  ->  only  buoy  signal  is  stored'+CRLF+ 

' .  default  IS  8  Hz': 

GET_NUMBER<  SAMPLERATE  .MESSAGE  )  i 

MESSAGE  :=  'INPUT  COUNTER  TIME  (IN  SECONDS)  WHEN  BUOY  IS  UERTICAL’: 
GET_NUMBER(  OFFSET,  MESSAGE  ): 

MEsiAGE  :=  'INPUT  PMT  SIGNAL  FREQUENCY  (IN  HERTZ)': 

GET_NUMBER(  SIGNALFREQ,  MESSAGE  >: 

END: 


( *»**•»«•««*••*•**«•**««»«»*•**•«*•**««««»*«*«»»*••*«*• ) 

(  ♦  *  ) 

(*  PROCEDURE  SHOW_SETTINGS  *) 

(  ■*  *  ) 

(♦  Displays  all  input  dates,  tines,  and  constants,  *) 

(*  and  prompts  user  to  accept  or  reject  values.  *) 

<*  Must  respond  with  a  Y  or  a  N,  This  response  is  ♦) 

<*  returned  to  the  calling  unit.  *) 

(  ♦  *  ) 

( *«••«****««*•*•**«••*««•*•«*•*•»•*•««««•****»«*«••*«•« ) 

PROCEDURE  SHOW_SETTINGS  (UAR  ANS  :  CHAR; 

TODAYSDATE.NEXTSTARTDATE  .FINALSTOPDATE  :  DATEREC; 
STARTTIME, STOPTIME. 

ONDURATION,  OFFDURATION  :  INTEGER: 

SAMPLERATE.  OFFSET.  SIGNALFREQ  :  REAL): 

UAR 

I  :  INTEGER:  {  loop  counter  > 

CSEC  :  INTEGER 

TIME  :  TIMEREC 

BEGIN 
PAGE: 

FOR  I  :=  1  TO  7  DO  WRITELN; 

WRITE(  'TODAYS  DATE  IS  .  '  ): 

WRITELN( TODAYSDATE. MONTH: 1  ,'/’  .TODAYSDATE.DAY: 1  ,'/'  .TODAYSDATE . YEAR : 1  ): 

URITE(  'CURRENT  TIME  OF  DAY  IS  .  '): 

SYSTIME(  TIME  ); 

WRITELN( TIME.H0UR:2  , ' : '  .TIME .MINUTE:2  . ' : ’  .TIME. CENT  I  SECOND  DIU  100:2  ): 


WRITE(  'DATA  COLLECTION  START  DATE  IS  .  ’)i 

WRITELN(NEXTSTARTDATE.MONTH:2  .NEXTSTARTDATE . DAY : 2  , 

NEXTSTARTDATE.YEAR:2  ): 

UIRITEC  'DAILY  STARTTINS  TIME  IS  .  '); 

C0NUERT_TIME5(  STARTTIME  .RIGHT  .TIME  ); 

WRITELN( TIME.HOUR:2  : ’ .TIME.MINUTE:2  ,  ’ : ' .TIME . CENTISECOND  DIU  100:2  )i 

WRITEI’DATA  COLLECTION  STOP  DATE  IS  .  '); 

WRITELN(FINALSTOPDATE.MONTH:2 . .FINALSTOPDATE . DAY : 2  . ' / '  . 
FINALST0PDATE.YEAR:2  >; 

WRITE( 'DAILY  STOP  TIME  IS  .  ')i 

CONOERT_TIMES(  STOPTIME  .RIGHT  .TIME  )» 

WRITELN(TIME.H0UR:2  . ' : '  .TIME . MINUTE : 2  : ’  .TIME . CENTISECOND  DIU  100:2  )i 

WRITE( 'LENGTH  OF  TIME  FOR  CONTINUOUS  BLOCK  OF  DATA  IS  . . , 


CONUERT_TIMES(  ONDURATION .RIGHT  .TIME  )j 

WRITELN<TIME.H0UR:2  ’  .TIME.MINUTE:2 . ’ : ’  .TIME. CENTISECOND  DIU  100:2); 

WRITE( 'OFF  TIME  BETUEEN  DATA  BLOCKS  IS  . 

CONUERT_TIMES(  OFFOURATION  .RIGHT  .TIME  ); 

WRITELN( TIME.H0UR:2  . ' : ’  .TIME . MINUTE : 2  .’ : '  .TIME . CENTISECOND  DIU  100:2  >; 

WRITE<’DATA  SAMPLING  RATE  IS  .  '>; 

WRITELN<  SAMPLERATE;2;2  >; 

WRITE<  'COUNTER  TIME  WITH  BUOY  UERTICAL  IS  .  '); 

WRITELN<0FFSET:4:3); 

WRITE<'PMT  SIGNAL  FREQUENCY  IS  . . .  '); 

URITELN(SIGNALFRE0:4:3); 

WRITELN; 

URITELN; 

WRITE  ('Do  you  want  to  change  any  values?  Y/N  ’>; 


REPEAT 

READ  (ANS)', 

UNTIL  (ANS  -  'N')  OR  (ANS  «  ’Y’)i 
END;  (*  PROCEDURE  SHOU_SETTINGS  *) 

(  *♦****«♦♦*♦♦**♦*♦*»*«♦*****♦♦*♦****•***■»■»*•»**♦****•♦♦♦**  ) 


(  *  *  ) 

(*  FUNCTION  NEXTDAY  *> 

(  •  *  ) 

(*  Increnents  the  date  passed  to  it  by  one  day.  ♦> 

(  *  *  ) 


( ***«««*«•«•**«**««•«•«•«•«•*•«*•*»«•••«««*«•«*•**•••**  ) 

FUNCTION  NEXTDAY  (DATE  :  DATEREC )  :  DATEREC; 

UAR 

DAYSINMONTH  :  28,  .31  ; 

BEGIN 

CASE  DATE. MONTH  OF 

2  :  DAYSINMONTH  28; 

4.8.9.11  :  DAYSINMONTH  30; 

OTHERWISE  DAYSINMONTH  31; 

END; 

DATE. DAY  (DATE. DAY  MOD  DAYSINMONTH)  +1;  <  next  day  > 

IF  DATE. DAY  -  1  THEN 

BEGIN  {  if  next  month  > 

DATE. MONTH  :-  (DATE. MONTH  MOD  12)  +  1; 

IF  DATE. MONTH  -  1  THEN  {  if  next  year  > 

-  DATE. YEAR  :-  DATE. YEAR  +  1; 

END; 


NEXTDAY  :=  DATEi  <  return  new  date  > 

ENDi  (*  FUNCTION  NEXTDAY  *) 


I  ) 


‘0 

Ui 


tiKd 


( *«**«»*****•**••««*«••**•*•*•«**•*«««*•*«*««**««**•«•« ) 

<  ♦  ♦ ) 

(*  FUNCTION  TIMEDATE  ♦> 

(  *  *  ) 

(*  Deternines  if  system  time  and  date  are  equal  to  *> 

(*  or  after  the  time  and  date  passed  into  function.  *) 

(*  If  so.  function  returns  true.  Otherwise  false  is  *) 

( *  returned .  *  > 

(  *  #  > 

<  ««***•«*«««•«««*«««••«*««*••»«*•«»•**«»«««*«•«**««««•«  ) 
FUNCTION  TIMEDATE  (TIME  :  INTEGERi  DATE  :  DATEREC )  :  BOOLEAN; 
UAR 

SDATE  :  DATEREC; 

SDAYS.  DAYS  :  INTEGER; 

BEGIN 

TIMEDATE  :=  FALSE; 

SYSDATE(  SDATE  ); 

SDAYS  :=  SDATE. YEAR*372  +  SDATE . M0NTH*3 1  +  SDATE. DAY; 

DAYS  DATE. YEAR*372  +  DATE.M0NTH*31  +  DATE. DAY; 

IF  SDAYS  >  DAYS  THEN 
TIMEDATE  TRUE 

ELSE 

BEGIN 

IF  SDAYS  =  DAYS  THEN 

IF  SYSCLOCK  >  TIME  THEN  TIMEDATE  :=  TRUE; 

END; 

END;  (*  FUNCTION  TIMEDATE  *) 


( **«««•«««••*•••«**««««**«*•*«•••«**«*•«•««*«*«»«««*«*» ) 

(  *  *  ) 

( *  PROCEDURE  PAUSE  ♦ ) 

<  *  *  ) 

(♦  Pauses  program  and  displays  current  time  and  ne)<t  *) 

(  ♦  starting  time.  *  ) 

(  ♦  ♦  ) 

(  ***•«««««•****«*«««••«**«««•«••'(  }••«««••«•«•**«*«*««•»  ) 
PROCEDURE  PAUSE  (NEXTONTIME  :  INTEGER;  NEXTSTARTDATE  :  DATEREC) 
UAR 


INTEGER; 

CURRENT  :  TIMEREC; 

TIME  :  TIMEREC; 

BEGIN 

PAGE; 

FOR  I  1  TO  10  DO 
WRITELN  ; 

WRITELN  (’WAITING  FOR  NEXT  DATA  COLLECTION  PHASE’); 
CONUERT_TIMES(  NEXTONTIME .RIGHT .TIME  ); 

URITELN( ’START  AT  ’  .TIME. HOUR:2  .TIME. MINUTE:2  . 

TIME.CENTISECOND  DIU  100  :2)| 

REPEAT  {  wait  for  designated  time  and  date  > 

SYSTIME  (CURRENT); 

WRITELN(’IT  IS  NOW  ’  .CURRENT . HOUR : 2  .CURRENT . MINUTE : 2  .  ’ 
CURRENT. CENTISECOND  DIU  100  :2); 
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iV*’-  •'  *1%  ''-S'  i’rV.V’VV<”'S'.  »*  .LN* 


UIRITE<  CHR(  31  ))i  <moves  cursor  up  one  position) 

FOR  I:=l  TO  100000  DO  BEGIN  END; 

UNTIL  TIMEDATEt  NEXTONTIME.  NEXTSTARTDATE  >i  {wait  till  tine) 
END!  (*  PROCEDURE  PAUSE  *) 


i 


'ft 

m 


( ********•**«**«*«*«•*«*««•»«••••*«•*««•*••*«««««*•*«•• ) 

( *  * ) 

( *  PROCEDURE  5T0REHEADER  * ) 

(  *  *  ) 

(*  Writes  a  header  to  disk  file.  Note  that  for  *) 

(♦  filename  to  be  a  unique  name,  sum  of  onduration  ♦) 

(*  and  off  duration  should  be  greater  than  one  hour.  ♦> 

(  •  ♦  ) 

( *•*«********«*•**•**««*••«•«••«•*«••*•«•*«•*«***«*«««*  ) 

PROCEDURE  STOREHEADER  < UAR  F  :  TEXT? 

NEXTSTARTDATE  :  DATEREC; 

NEXTONTIME, 

ONDURATION  :  INTEGER; 

SAMPLERATE, 

ROLLAMPLITUDE. 

ROLLPERIOD  :  REAL)  i 

UAR 

TIME  :  TIMERECi 
FILENAME  :  STRING  C20]i 
T  :  INTEGER! 

BEGIN 

FILENAME  :=  ’#4:B’! 

STRWRITE  < FILENAME, B,T, NEXTSTARTDATE. MONTH  :0)! 

FILENAME  FILENAME  +  i 

STRWRITE  ( FILENAME. T+1  ,T. NEXTSTARTDATE. DAY  :0>! 

FILENAME  FILENAME  +  ’H’l 

STRWRITE  (FILENAME  .T+1  .T  .NEXTONTIME  DIU  3G0000  :0)! 

FILENAME  :«  FILENAME  +  ’.TEXT’; 

REWRITE  (F,  FILENAME)! 

WRITE  (F.’DATE  iMM/DD/YY)  »  ’  .NEXTSTARTDATE . MONTH : 0 ) i 
WRITE  <F.’/’  , NEXTSTARTDATE. DAY:0)| 

WRITELN  <F.’/’  , NEXTSTARTDATE. YEAR:0)! 

WRITELN<PTR  , 'START  DATE  :  ’  .NEXTSTARTDATE. MONTH : 2 .NEXTSTARTDATE .DAY : 2  , 
. NEXTSTARTDATE. YEAR:2)! 

WRITELN  (F, 'START  TIME  (csec)  =■  ’.NEXTONTIME  :0)i 
CONUERT_TIMES( NEXTONTIME .RIGHT .TIME ) j 

WRITELN(PTR, 'START  TIME  :  ’  .TIME. HOUR:2  .TIME. MINUTE:2  , 

TIME.CENTISECOND  DIU  100  :2>! 

WRITELN  (F.'ON  DURATION  (csec)  -  ’  .ONDURATION: 0  )  i 

WRITELN(PTR  ,'DATA  BLOCK  LENGTH  :  '.ONDURATION  DIU  B000  :4.'  MINUTES')! 
WRITELN  (F, 'SAMPLE  RATE  •  '  .SAMPLERATE  ) i  {  rate  buoy  angle  is  recorded  > 
WRITELN(PTR,' SAMPLE  RATE  :  ’  .SAMPLERATE: 4 , ’  Hz'); 

WRITELN  (F.'ROLL  AMPLITUDE  -  ’.ROLLAMPLITUDE)! 

WRITELN  (F.'ROLL  PERIOD  -  ’.ROLLPERIOD)! 

ENDi  {  STOREHEADER  > 


a 


( * . * . * . . . . . . . 

( *  * ) 

( *  PROCEDURE  PLOTANGLE  *  ) 

(  *  *  ) 

(*  Plots  a  line  in  window  of  CRT.  If  right  hand  side  ♦) 


li 
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( *  of  CRT  is  net  then  screen  is  cleared  and  pen  noved  ♦  ) 

(*  to  left  nargin  of  screen.  ♦) 

(  •  *  ) 

( *•**•«**•****«***•*«*««*»••»*«**«***«*»•«•«««•*•*«••*** ) 

PROCEDURE  PLOTANGLE  (ANGLE  :  REAL;  UAR  TIME  :  INTEGER;  DTIME 
BEGIN 

IF  TIME  >=  3000  THEN 
BEGIN 

TIME  •-  0;  <  set  for  left  margin  > 

MOUE  < TIME  .ANGLE  ) ;  {  move  pen  to  left  nargin  > 

CLEAR_DISPLAY ;  {  clear  graphics  display  > 

END 

ELSE 

BEGIN  {  plot  line  > 

TIME  :«  TIME  +  DTIME; 

LINE  ( TIME  .ANGLE  );  {  draw  line  > 

END; 

END;  (*  PROCEDURE  PLOTANGLE 

( ***««*•*»•*«*•«««•***«*«*««*««•*«•«»*««««•*•««•««*««*«*  ) 

(  *  *  ) 

(♦  GRAPHICS  PLOT_AXES  *) 

(  *  *  > 

(♦  Initiates  graphics,  sets  viewport.  Draws  axes  *> 

(*  for  buoy  angle  versus  tine  plot.  Axes  is  drawn  *> 

(*  on  alpha  screen,  plotting  is  done  on  graphics  ♦> 

(*  screen.  This  saves  tine  in  clearing  screen.  •) 

(  •  «  > 

( *«•*••**•••««*«•**«**««•*«•«•«««*••••««»•»»»*»»*»»**»** ) 
PROCEDURE  PL0T_AXES(  MAXANGLE  :  REAL  ); 

UAR 

I  :  INTEGER; 

DY  :  REAL; 

BEGIN 

PAGE;  <  clear  alpha  screen  > 

WRITELN; 

SET_ASPECT  (511  ,389);  <  set  graphics  screen  > 

SET_UIEWPORT  ( 0 . 1 05 .0. 930 .0 . 1 64 .0.692 ); 

SET_WINDOW  (0  .3000  , -MAXANGLE  .MAXANGLE  >; 

DY  MAXANGLE  /  4.5; 

FOR  I  4  DOUNTO  1  DO  BEGIN  {  display  upper  y  axes  > 
WRITELN( '  -I’); 

WRITELN  (  I*DY:3:0,’  — I  ’  ); 

END; 

WRITELN( ’  -I’); 

WRITE  (0:3,'--!’);  {  display  zero  degree  line  > 

FOR  I  I  TO  71  DO 
WRITE  (’-'); 

WRITELN; 

FOR  I  :=  1  TO  4  DO  BEGIN  <  display  lower  y  axes  > 

WRITELN( ’  -I’); 

WRITELN  (-I*DY:3:0,'--I ’ ); 

END; 

WRITE!'  -1’);  {  display  x  axes  > 

FOR  I  :»  1  TO  71  DO 
WRITE  (’_’); 


INTEGER  ); 


.  ‘‘r  *  ‘  l. 
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<  ticks  on  X  axes  > 


URITELN; 

FOR  I  I  TO  n  DO 
WRITE  (  ’  ' ’  )i 

WRITELNi 

FOR  I  :=  0  TO  10  DO  {  number  x  axes  > 

WRITE  (1*3:7  )| 

ENDi  (•  GRAPHICS  PROCEDURE  PLOT  AXES  *) 


(  *  ♦  ) 

<  *  PROCEDURE  TAKEDATA  *  ) 
(  *  *  ) 
(*  Records  time  interval  on  HPB334A  and  converts  *) 
(*  to  buoy  angle.  Readings  taken  at  regular  *) 
(*  intervals  specified  by  sample  frequency.  Sends  ♦) 
(*  readings  to  silo.  If  sample  freq.  3H2  or  less  *) 
(*  then  also  records  waverider  receiver  output  volts  *) 
(*  Plots  buoy  angle  on  CRT  in  real  time.  ♦> 
(  *  »  ) 


( •«««•«*«••*•«•*«•••**««»***«•««««••«•«•«••»•••*••****«  ) 

PROCEDURE  TAKEDATA  (UAR  F  :  TEXT; 

ONDURATION  :  INTEGER; 

SAMPLERATE  . 

SIGFREQ, 

OFFSET  :  REAL); 


{Note:  for  this 
between  signal 
motor  shaft  to 
UAR 


to  work  properly,  there  must 
and  reference  wave  when  buoy 
adjust  phase  shift.) 


exist  a  SO-degree  phase 
is  vertical.  Rotate  gear 


shi  ft 
on 


B 

W 

I.  T 

DATA 

BUOYANGLE 

BU0YANG2 

WAUEHEIGHT 

SAMPTIME 

NUMREADINGS 

INCLUDEWAUE 

NEWOATA 

OUERFLOU 

ERRORRETURN 

TIME 

TDATA 

0LDI5RH00K 


FILE  OF  REAL; 
TEXT; 

INTEGER; 
STRINGI21 ]; 
REAL; 

REAL; 

REAL; 

INTEGER; 

INTEGER; 

BOOLEAN; 

BOOLEAN; 

BOOLEAN; 

INTEGER; 

INTEGER; 

TIMERDATA; 

KBDHOOKTYPE; 


{  counter  > 

{  input  buffer  } 


{  time  in  centisecs  between  readings  > 

{  counter  of  readings  taken  } 

{  true  if  waverider  signal  is  recorded  > 


{variable  for  initialization  outcome) 
{time  in  centiseconds  from  plot  origin) 


{•»•**•*««***••«**••••••*••*•••••••••••••****««**««««*•»•«*•«•*••••*•••««• 

PROCEDURE  GETREADING  (UAR  STATBYTE .  DATABYTE  :  BYTE;  UAR  DOIT  :  BOOLEAN); 
{this  procedure  uses  global  data  from  TAKEDATA) 

BEGIN 

IF  OUERFLOW  THEN 
BEGIN 

URITELNC ’DIGITIZING  FREQUENCY  TOO  HIGH’); 

lOCONTROL  ( GPIO ,3 ,CLOSE_LENS ) ;  {  close  lens  and  shut  off  engine  ) 

CL0SE(F>;  {close  without  saving) 

END; 


OVERFLOW  :=  TRUE;  {if  this  is  not  reset  to  false,  proc  never  finished} 
READSTRINO^  HP5354A ,  DATA  )) 

STRREAD(  DATA  .2  .T  .BUOYANGLE  )  •, 

READSTRING(  HP5334A .  DATA  ); 

STRREAD<  DATA  .2  .T  .BUOYANGZ 

BUOYANGLE  :=  (BUOYANGLE  +  BUOYANGZ)  /  2i  {average  out  some  of  noise) 

BUOYANGLE  :=  (BUOYANGLE  -  OFFSET)  *  SIGFREQ  *  180; 

WRITE(  B  ,  BUOYANGLE  ); 

IF  INCLUDEWAUE  THEN 
BEGIN 

READSTRING(  K48S .  DATA  ); 

WRITELN  (W,  DATA); 

END; 

NEWDATA  :=  TRUE; 

OVERFLOW  :=  FALSE; 

END; 

{ft*************************************************************************} 

BEGIN  {TAKEDATA} 

{REWRITE  (B  .  ’RAM:BDATA’  ); 

REWRITE  <W.  ’RAM:WDATA’  );>  REWRITE  <B);  REWRITE  (W); 
lOCONTROL  ( GPIO  ,3  ,0PEN_LENS  ) ;  {  open  lena  ,  start  engine  > 

FOR  I:  =  t  TO  200  DO  READSTRINGf  HP5334A .  DATA  );  {  waits  about  20  sec  > 

GRAPHICS_INIT; 

DISPLAY_INIT( CRTADDR .CONTROLWORD  .ERRORRETURN  ) ; 

IF  ERRORRETURN  <>  0  THEN  HALT;  {  if  error  halt  > 

PL0T_AXES( MAXANGLE  / ;  <  plot  graph  axes  on  alpha  screen  > 

TIME  :=  3000; 

NUMREADIN6S  :=■  R0UND(  (  ONOURATION  DIV  1  00  )  *  SAMPLERATE  ) ; 

SAMPTIME  ROUND  (100  /  SAMPLERATE);  {time  in  csec  between  samples} 

IF  SAMPLERATE  <>  3  THEN  INCLUDEWAVE  TRUE 
ELSE  INCLUDEWAVE  :=  FALSE; 

PLOTANGLE  ( 0 . 0  .TIME  .SAMPTIME  ) ;  <  place  pen  > 

OVERFLOW  :=  FALSE; 

NEWDATA  :=  FALSE; 

BUOYANGLE  :=  0; 

TRY 

OLDISRHOOK  :»  TIMERISRHOOK ;  {save  old  ISR  to  later  restore) 
TIMERISRHOOK  GETREADING;  {assign  new  ISR) 

TDATA, COUNT  SAMPTIME; 

CALL(  TIMERIOHOOK.  CYCLICT.  SETT,  TDATA  );  {sot  cyclic  timer) 

CALL(  MASKOPSHOOK.  TIMERMASK,  0  );  {enable  timer) 

I  :=  0; 

REPEAT 

IF  NEWDATA  THEN 
BEGIN 

I  I  +  1 ; 

NEWDATA  FALSE; 

PLOTANGLE  ( BUOYANGLE  .TIME  .SAMPTIME  ) ; 

END; 

UNTIL  I  >-  NUMREADINGS; 

ESCAPE  (0);  {generate  an  error) 

RECOVER 

BEGIN 


CALL!  MPiSKOPSHOOK  .  0,  TIMERMASK  );  {disable  tiner  isr) 

TDATA. COUNT  :=  0i 

CALK  TIMERIOHOOK,  CYCLICT,  SETT,  TDATA  );  {zero  timer) 
TINERISRHOOK  :=  OLDI5RHOOK;  {restore  old  ISR) 

END: 


CLEAR_D I  SPLAY ;  {  clearm  graphics  from  CRT  > 

GRAPHICS_TERM ;  {terminate  the  graphics  package) 

lOCONTROL  ( bPIO  .3  ,CLOSE_LENS  )  i  {  close  lens  and  shut  off  engine  ) 


RESET  (B);  {transfer  data  to  disk) 

RESET  (W); 

WHILE  NOT  EOF  (B)  DO 
BEGIN 

READ  (B.  BUOYANGLE): 

WRITELN(  F.  BUOYANGLE  ); 

IF  INCLUDEWAUE  THEN 
BEGIN 

READLN  (W,  DATA); 

STRREAD(  DATA  ,5 .T .UAUEHEIGHT  ); 

WAUEHEIGHT  WAUEHEIGHT  *  UAUECAL; 

WRITELN(  F  .’U’  .WAUEHEIGHT  ); 

END; 

END; 

CLOSE  (B  .  ’SAUE’  ); 

CLOSE  (W  .  ’SAUE’  ); 


END;  (*  PROCEDURE  TAKEDATA  *) 


u 


i  I 


( •*«*«***»••*•****«**********»*•«••*««»«««•*«««****««*• ) 

(  *  *  ) 

( *  MAIN  PROGRAM  * ) 

(  *  *  ) 

(  ♦*♦•»♦**■»******♦**♦•***♦***«**********■»*♦*«♦****♦**♦***  ) 

BEGIN 

WRITELN  ('HAUE  YOU  INSTALLED  A  RAM  UOLUME  NAME  "RAM:"?  (Y)  or  (N)’); 
READ  (ANS): 

IF  (ANS  -  ’N’  )  OR  (ANS  =■  ’n’)  THEN  HALT; 

WRITELN; 

REWRITE(PTR  , ’PRINTER: ’  ); 

lOINITIALIZE;  {reset  all  interfaces) 

lOCONTROLf GPIO  ,3  ,CLOSE_LENS > ;  {insure  lens  is  closed  first  thing) 
SYSDATE  ( TODAYSDATE  ) ;  {  set  current  date  to  system  date  ) 

NEXTSTARTDATE  TODAYSDATE;  {  default  settings  ) 

FINALSTOPDATE  :=  TODAYSDATE; 

STARTTIME  :=  75G0000;  {21:00:00) 

STOPTIME  :=  8280000;  {23:00:00) 

ONDURATION  :=  180000;  {00:30:00) 

OFFDURATION  :=  360000;  {01:00:00) 

SAMPLERATE  :=  5; 

OFFSET  :-  1G.0E-3; 

SIGNALFREQ  31.5; 

REPEAT  {  get  user  provided  inputs  ) 

SET_0ATES(  TODAYSDATE,  NEXTSTARTDATE.  FINALSTOPDATE); 

SET_TIMES(  STARTTIME,  STOPTIME.  ONDURATION,  OFFDURATION); 
SET_DATARATES  (SAMPLERATE.  OFFSET,  SIGNALFREQ); 


WRITELN  (’INPUT  BUOY  TYPE  (8x26,  5x11,  etc.)'  )) 

READLN  (BUOYTYPE)i 

SHOW_SETTINGS  ( ANS  ,  TODAYSDATE  .  NEXTSTARTDATE  ,  F INALSTOPDATE  , 

STARTTIME,  STOPTIME, 

ONDURATION,  OFFDURATION. 

SAMPLERATE,  OFFSET,  SIGNALFREO); 

UNTIL  ANS  =  ’N’  ; 

WRITELN; 

WRITELN  (’****INSERT  DATA  OUTPUT  DISK  IN  RIGHT  HAND  DRIUE****'); 

READ  (ANS): 

NEXTONTIME  :=  STARTTIMEi 

IF  STARTTIME>STOPTIME  THEN  -(crosses  midnight) 

NEXTSTOPDATE  NEXTDAY(  NEXTSTARTDATE  ) 

ELSE 

NEXTSTOPDATE  :=  NEXTSTARTDATE; 

REPEAT  {  go  until  final  stop  date  and  time  > 

PAUSE  (  NEXTONTIME.  NEXTSTARTDATE  >;  {  wait  till  collection  phase  > 

STOREHEADER  ( F  .NEXTSTARTDATE  .NEXTONTIME  .ONDURATION  .SAMPLERATE  ,0,0); 
WRITELN  (F  .BUOYTYPE ,’  BUOY  ANGLE  US  TIME’); 

TAKEDATA  ( F  .ONDURATION  .SAMPLERATE  .SIGNALFREQ  .OFFSET  ) ; 

CLOSE  (F  .  ’SAUE’  ); 

WRITELN(PTR  , 'END  OF  DATA  ACQUISITION’); 

WRITELN(PTR  ); 

WRITELNIPTR); 

NEXTONTIME  (NEXTONTIME  +  ONDURATION  +  OFFDURATION)  MOD  MIDNIGHT; 
SYSDATE  (TODAYSDATE);  {update  todaysdate) 

IF  NEXTONTIME  <  SYSCLOCK  THEN  {crosses  midnight) 

NEXTSTARTDATE  NEXTDAY(  TODAYSDATE  ) 

ELSE  {do  not  cross  midnight  yet) 

NEXTSTARTDATE  TODAYSDATE; 

IF  TIMEDATE(  STOPTIME.  NEXTSTOPDATE  >  THEN  {end  of  nights  collection) 
BEGIN 

NEXTONTIME  STARTTIME; 

IF  STOPTIME  >  STARTTIME  THEN  {collection  does  not  cross  midnight) 
BEGIN 

NEXTSTARTDATE  :•  NEXTDAY(  NEXTSTOPDATE  ); 

NEXTSTOPDATE  NEXTSTARTDATE; 

END 

ELSE  {collection  does  cross  midnight) 

NEXTSTOPDATE  NEXTDAY(  NEXTSTARTDATE  ); 

END; 

UNTIL  TIMEDATE(  STOPTIME,  F INALSTOPDATE  ); 

PAGE;  {  clear  screen  ) 

FOR  I  :=  1  TO  10  DO  {  terminate  message  ) 

WRITELN; 

SYSTIME  (TIME); 

WRITE  (’Program  Terminated  at  ’  .TIME . HOUR : 2  ,  ’ : ’  ) ; 

WRITELN  ( TIME. MINUTE:2  .TIME. CENTISECOND  DIU  100:2); 

END.  (*  PROGRAM  BUOY_MOTION  *) 


Program  BUOY_PROBABILITIES 


PROSRAM  BUOY_PROBABILITIES  (INPUT,  OUTPUT); 

(Determines  the  probabilities  of  detection  and  recognition  of  flashing  lights 
on  rolling  buoys.  Input  data  is  in  user  specified  file  on  left  disk.  Data 
was  stored  using  either  PROGRAM  BU0Y_M0TI0N  or  PROGRAM  ROLLSIM  &  must  be  on 
left  disk.  Output  stored  on  right  disk  (i4:).  LANT.REAL  must  be  on  left  (43:) 
disk.  LANT.REAL  was  generated  with  TRANSLAMP . CODE .  Before  executing  program 
operator  (or  Boot  disk)  must  set  up  a  ram  volume  named  RAM:  . 

By  Dan  Broun.  Apr  1987.> 

IMPORT 
RND . 

lODECLARATIONS  , 

SYSDEUS  , 

DGL_LIB , 

DGL.INQ; 

CONST 

BUOYRATE 
NUMANGLES 
RAMFILE 
NUMDISTS 
NUMTIMES 
08STIMEINC 
RATEMULTIPLIER 
WINOOWLENGTH 
NUMUINOOWS 
ESC 

TYPE 

CHARACTERISTICS  =  ( FLS  ,  FL4  .  FL2S  .  QKFL  ,  GPFL5  .  GPFL6  .  IQKFL,  MOA  . 

EINTG,  0CC4)i  (the  ten  standard  flash 

characteristics  I'isted  on  page  6-52  of  AtoN  Technical> 
MSG_TYPE  -  STRINGC255]; 

0IR_TYPE  ■  (  RIGHT,  LEFT  );  (enumerated  direction  type> 

PROBRECORD  -  RECORD 

POD  :  ARRAY  [ I . .NUMDISTS  ,  1 .. NUMTIMES 1  OF  REAL; 

POR  :  ARRAY  M .. NUMDISTS  .  1 .. NUMTIMES  1  OF  REAL; 

EFI  :  ARRAY  C 1 . .NUMDISTS  1  OF  REAL; 

DISTANCEINC  :  REAL; 

TIMEINC  :  REAL; 

END; 

PERRECORD  -  RECORD 

DETECDIST  :  REAL; 

RECOGDIST  :  REAL; 

END; 

PERIODINFO  -  ARRAY! 1 . .NUMTIMES]  OF  PERRECORD;  (gives  detection 

and  recognition  distance 
of  each  flash  period  in 
a  window  of  data) 

HEADRECORD  -  RECORD 

DAY  :  1..3I; 

MONTH  ;  1 .  .  I  2  ; 

YEAR  :  86.. 90; 

STARTTIME  :  INTEGER;  (time  in  centiseconds  past  midnight) 


>  6;  (maximum  buoy  angle  sample  rate) 

•  455;  (number  inten  vs  angle  data  points  in  LANTARRA') 

“  ’RAM:DATA.REAL’ ;  (file  name  of  inten  vs  time) 

■  16;  (number  of  different  distances  from  light) 

-  15;  (number  of  different  observation  times) 

-  2;  (integer  no.  seconds  per  time  increment) 

-  7;  (no.  sample  points  multiplied  for  better  resolution) 

-  OBSTIMEINC  *  NUMTIMES  *  RATEMULTIPLIER  •  BUOYRATE; 

■  500;  (number  data  windows,  NUMTIMES  wide,  analyzed) 

-  CHR(27); 
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ONDURATION  :  INTEGER)  {data  block  length  in  cent i seconds > 
SAMPLERATE  :  REAL)  {rate  in  Herzl 

ROLLAMPLITUDE  :  REALi  {amplitude  in  degrees  off  vertical} 
MEANROLLANGLE  :  REAL)  {expected  value  of  tilt  angle} 

ROLLPERIOD  :  REAL)  {period  in  seconds} 

END) 

LANTRECORD  «  RECORD 

ANGLE  :  REAL) 

INTENSITY  :  REAL) 

END) 

LANTARRAY  =  ARRAY  C 0 . . NUMANGLES I  OF  LANTRECORD)  {array  of  intensity  vs 

angle  from  user  specified 
lantern  file} 

FILENAMETYPE  =  STRING  1201) 

FILENAMEARRAY  -  ARRAY  Cl.. 10]  OF  FILENAMETYPE)  {array  of  disk  locations 

of  data  sets  to  be  analyzed} 

DATAUINDOWS  -  RECORD 

INTENSITY  :  ARRAY  [ 0 . . WINDOWLENGTH ]  OF  REAL) 

SAMPLERATE  :  REAL) 

MAXEFI  :  ARRAY  1 1 . .NUMTIMESl  OF  REAL)  {detection  EFI} 
MINEFI  :  ARRAY  1 1 . . NUMTIMES ]  OF  REAL)  {recognition  EFI} 
END) 

PROBFILE  -  FILE  OF  PROBRECORD) 

FLASHARRAY  »  ARRAY  [0..100  *  RATEMULTIPL lER ]  OF  0..1) 


DATAUINDOWS 
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PROBABILITY 

HEADEP 

FILENAMES 

LANTERNFILE 

MINUTES 

FLASHCHAR 

OUTDISK 

{follow 

FILEPOS 

SEED 

RANGE 

NUMREA0ING5 

{  f  ol  low 

UINOOU 

FLASH 

FLASHLEN 

W 

STARTPOS 

TIME 


:  INTEGER) 

:  PROBRECORD) 

:  HEADRECORD) 

:  FILENAMEARRAY) 

:  FILENAMETYPE)  {file  name  of  inten  vs  angle} 

:  INTEGER)  {number  of  minutes  of  data  set  to  analyze} 

:  CHARACTERISTICS)  {flash  characteristic  being  analyzed} 

:  PROBFILE) 

ing  variables  refer  to  intensity  record  on  ram  disk} 

:  INTEGER)  {file  position  (in  ram  file)  of  data  window} 

;  INTEGER)  {dummy  variable  for  random  number  generator} 

INTEGER)  {range  of  random  numbers  for  file  position} 

:  INTEGER)  {number  of  data  points  stored  on  ram  file} 
ing  variables  refer  to  WINDOW. INTENSITY  array} 

:  DATAUINDOWS)  {15  periods  of  intensity  vs  time  data} 

:  FLASHARRAY)  {normalized  flash  period} 

:  INTEGER)  {number  of  data  points  in  period  in  FLASH} 

:  1 . .NUMWINDOWS)  {counter  for  number  of  windows  analyzed} 
;  INTEGER)  {start  position  in  FLASH  for  normalizing  win.} 
:  I.. NUMTIMES)  {counter  to  zero  arrays) 

:  TEXT) 


PROCEDURE  PL0T_WIND0W  < UAR  W  :  DATAUINDOWS)) 

{Plots  data  array  pointed  to  by  P.  Used  for  testing.} 

UAR 

I  :  INTEGER)  {loop  counter} 

CHARWIDTH,  CHARHEIGHT  :  REAL) 

XMIN,  XMAX,  YMIN,  YMAX  ,  X.  Y  :  REAL)  {x  and  y  plot  values) 
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FIRSTPOINT  :  BOOLEAN? 

MAXMAG  :  REAL;  <maximun  magnitude  in  data) 

{•***««**'»***«*«*****««««*««*««*«#*«***4«******«««**««****«*«***«*»*********> 

PROCEDURE  INITIALIZE_PLOT; 

CONST 

CRTADDR  -  3;  <3-CRT.  705=PLOTTER> 

PLTADDR  -  705; 


CONTROLWORD  *  0; 
GET  ASPECT=254; 


ERROR  :  INTEGER; 

DUMMY  :  INTEGER; 

RATIO_LIST  :  ARRAY! 1.. 21  OF  REAL; 

BEGIN 

GRAPHICS_INIT; 

DISPLAY_INIT( CRTADDR .CONTROLWORD .ERROR ) ; 

INQ_WS( GET_ASPECT .0.0.2 .DUMMY .DUMMY .RATIO_LIST .ERROR ) ; 

IF  ERROR-0  THEN  SET_ASPECT< I . 0 .RATIO_LISTC 2 ] ) : 

■CSET_UIEWPORT(0. 0.1  .00.0.0.0.74);}  {Sets  plotting  area  of  CRT) 

SET_UIEWPORT(0.07 .1 .00  .0.04  .0.7);  {Sets  plotting  area  of  CRT) 

END: 

{t*********************************************#****************************} 

PROCEDURE  SET_USER_SCALE  ( UAR  W  :  DATAWINDOWSi 

UAR  XMIN. 

XMAX  . 

YMIN. 

YMAX  , 

CHARWIDTH, 

CHARHEIGHT  :  REAL > ; 


VALUE.  MAXUALUE  :  REAL; 

I  :  INTEGER; 

BEGIN 

MAXUALUE  0; 

WITH  W  DO 

FOR  I:=0  TO  WINDOWLENGTH  DO  {find  maximum  value  in  data  array) 

BEGIN 

VALUE  INTENSITY  CIl; 

IF  VALUE>MAXVALUE  THEN  MAXUALUE : -VALUE ; 

END; 

XMIN:-0.0; 

XMAX:^  30;  {number  of  seconds) 

YMIN:-  0; 

YMAX:-  MAXUALUE; 

SET_WINDOW(XMIN. XMAX. YMIN, YMAX);  {Sets  user  coord,  scale) 

CHARWIDTH  :-  < XMAX  -  XMIN)  /  50; 

CHARHEIGHT  :=  (YMAX  -  YMIN)  /  25; 

SET_CHAR_S I ZE( CHARWIDTH .CHARHEIGHT ) ;  {sets  constant  relative  size) 

END; 

<*****♦**♦*♦*♦**»#•#*♦*#**#♦*********♦*»*«**♦******#*»*■»#*********»**♦*****»} 

PROCEDURE  LABEL_PLOT  <  XMIN, 

XMAX  , 

YMIN  , 

YMAX  , 

CHARWIDTH. 

CHARHEIGHT  :  REAL  ); 
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{vertical  axis  tic  marks) 


TITLE  :  STRIN6I80]-, 

X ,  DX ,  Y .  DY  :  REAL; 

I  ,  J ,  K  :  INTEGER; 

BEGIN 

MOL'E(XMAX  .0  ); 

LINE(0,0); 

X  :=  0; 

DX  :=  1.0;  {Label  x  axis) 

WHILE  X  <  XMAX  DO  {major  tick  marks,  horizontal  axis) 

BEGIN 

X  :=  X  +  DX; 

M00E(X  ,0  ); 

LINE(X  ,-CHARHEIGHT  >; 

TITLE:*’’;  {label  major  ticks) 

STRWRITE(TITLE,1  ,K  .X:0:0); 

MOOE(X-STRLEN( TITLE  )*CHARWIDTH+CHARWIDTH/2 .-1 . 2*CHARHEI6HT ) ; 

GTEXT( TITLE) 

END; 

MOOE<0,YMAX  ); 

LINE(0,YMIN); 

LINE(-CHARWIDTH,YMIN); 

Y  :=  YMIN; 

DY  :*  ( YMAX  -  YMIN)/10;  {vertical  axis  tic  marks) 

FOR  I: -I  TO  10  DO 
BEGIN 

Y  :*  Y  +  DY; 

MOOE(0 .Y  >; 

LINE(-0.5  *  CHARWIOTH.Y); 

END; 

TITLE:*’ ’ ; 

STRWRITE(TITLE,1  ,K,YriAX:4:t  ); 

M0VE( -STRLEN( TITLE  )*CHARWIDTH ,YMAX ) ; 

GTEXT<  TITLE  ); 

TITLE  :=  ’Data  Window’; 

MOOE  (  <XMAX  -  XMIN  -  STRLEN< TITLE  )*CHARWIDTH  )/3  +  XMIN,  YMAX  ); 

GTEXT( TITLE  ); 

END; 

<*«••««•«««<•««••*•»«««*«««*««««•««««»««•«*•«•«••«««•««««••«•«•««•«•«•**•*••} 

BEGIN  {PL0T_DATA) 

INITIALIZE_PLOT; 

SET_USER_SCALE  ( W  ,XMIN  ,XMAX  .YMIN  ,YMAX  .CHARWIDTH  .CHARHEIGHT  > ; 

LABEL_PLOT  <  XMIN  ,XMAX  .YMIN .YMAX .CHARWIDTH .CHARHEIGHT ) ; 

MOOE(0  ,  W. INTENSITY  [0] ); 

I  :=  0; 

WHILE  X  <  30.0  DO 
BEGIN 

I  :=  I  +  1 ; 

X  :-  I  /  W.SAMPLERATE; 

LINE( X  ,W. INTENSITY  Cl]); 

END; 

CALL<  DUMPGRAPHICSHOOK ) ; 

CLEAR_DISPLAY; 

GRAPHICS_TERM; 

WRITELN  (PRIN,  ’MAXEFI  MINEFI’); 

FOR  I: -I  TO  NUMTIMES  DO 
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URITELN  (PRIN,  U.MAXEFI  [11  :0:2.  W.MINEFI  [I]  :40;2)i 
PAGE  (PRIN): 

ENDi 

{*«*««*****»********««««*««««»*««*««»«**«*«««**««***ii'««**«*****«**«**«***«*«*} 

{♦I***************************************************************************} 

{•«**»*««««*««**«*««*«*«*«**««««***««««**««*«*******«««*««*****4**«*«««*«****) 

PROCEDURE  INPUT_USER_DATA  ( UAR  FILENAMES  :  FILENAMEARRAY ! 

UAR  LANTFILE  :  F ILENAMETYPE  ) ; 

{FILENAMES  specifies  data  filenames  to  be  analyzed,  LANTFILE  specifies  the 
defined  lantern  type  (vertical  intensity  profile)  to  use.} 

UAR 

NUMSETS  :  1  .  . 10) 

I  .  J  :  INTEGER) 

C  :  CHAR) 

BEGIN 

FOR  I:=l  TO  10  00  WRITELN) 

WRITELN  (’HAUE  YOU  INSTALLED  A  RAM  VOLUME  NAMED,  "RAM"  ^  (Y>  or  (N>’)j 
READ  (C); 

IF  (C  =  ’N’  )  OR  (C=’n’)  THEN  HALT; 

PAGE) 

FOR  I:=1  TO  10  DO  WRITELN; 

URITELN<  ’ INPUT  NAME  OF  LANTERN  FILE  (excluding  volume  no.)’); 

READLN  (LANTFILE)) 

LANTFILE  :»  ’#3:’  +  LANTFILE; 

PAGE) 

FOR  I :=1  TO  10  DO  WRITELN; 

WRITELN( ’ INPUT  NUMBER  OF  SETS  OF  DATA  TO  ANALYZED  (1..10)’)) 

READLN  (NUMSETS); 

FOR  J:»  1  TO  NUMSETS  DO 
BEGIN 
PAGE) 

FOR  I:-1  TO  10  DO  WRITELN) 

WRITELN( ’ INPUT  DATA  FILE  NAME  (excluding  vol .  no.)  FOR  RUN  NUMBER  ’,J)) 
READLN(  FILENAMES  [JI  )) 

FILENAMES  [J]  ’#3:’  +  FILENAMES  IJl) 

END) 

FOR  J:-  NUMSETS  +  1  TO  10  DO 

FILENAMES  [J]  :*  ’NONE')  {marks  end  of  sets} 

PAGE; 

FOR  I :-1  TO  10  00  WRITELN; 

WRITELN  (’  ’:22,  ’PLACE  DATA  INPUT  DISK  IN  LEFT  HAND  DRIVE _ (#3)’); 

WRITELN  (’  ’:22,  ’THIS  DISK  MUST  ALSO  HAVE  LANTERN  FILE  ON  IT’); 

WRITELN  (’  ’:22,  ’PLACE  OUTPUT  DISC  IN  RIGHT  HAND  DRIVE _ (#4)’); 

WRITELN) 

WRITE  (’  ’:22.  ’TYPE  ANY  CHARACTER  WHEN  READY...  ’); 

READ  (C); 

WRITELN) 

END) 

<*»♦♦***♦♦#*** ************** ♦**#*♦#***♦*♦*♦*****«*»************ *♦**»*»*♦*****} 
{****************************************************************************} 
<****************************************************************************} 

PROCEDURE  INITPROB_RECORD  (VAR  PROBABILITY  :  PROBRECORD ) ) 

{sets  all  probabilities  to  zero,  sets  the  observer  distance  and  effective 
intensity  fields} 

VAR 


DIST  :  INTEGER) 
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OBSTIME  :  INTEGER; 
BEGIN 

WITH  PROBABILITY  DO 
BEGIN 

DISTANCEINC  :=  0.5i 


TIMEINC  :=  OBSTI 
EFI  [ 1  ]  0.  IBS 

EFI  [2]  :=  0.905 
EFI  13]  :=  2.3G8 
EFI  [4]  4.894 

EFI  [51  :=  8.889 
EFI  CGI  :=  14.88 
EFI  [7]  :=  23.54 
EFI  C8]  :=  35.74 
EFI  [9]  :=  54.59 
EFI  1101  :»  75.4 
EFI  cm  :=  106. 
EFI  [121  :=  146. 
EFI  [13]  200. 

EFI  [14]  :=  270. 
EFI  [15]  360. 

EFI  [16]  :=  470. 
FOR  DIST:=1  TO  N 
FOR  OBSTIME: 
BEGIN 


OBSTIMEINCi 

0.195; 


0 . 905 ;  { 1 . 0  n 

2.368;  <1 .5  n 

4.894;  <2.0  n 

8.889;  <2.5  n 

14.881;  <3.0 

23.545;  <3.5 

35.749;  <4.0 

54.597;  <4.5 

75.484;  <5.0 

106.176;  <5. 

146.888;  <6. 

200.399;  <6. 

270.177;  <7. 

360.545;  <7.! 

470.871;  <8. 

TO  NUMDISTS  DO 
TIME:-1  TO  NUMTIMES  DO 


<0.5  n .  tni  les  > 
<1.0  n. miles > 
<1.5  n.mi les> 
<2.0  n. miles} 
<2.5  n. miles} 
<3.0  n. miles} 
<3.5  n. miles} 
<4.0  n . mi les } 
<4.5  n. miles} 
<5.0  n. miles} 
<5.5  n. miles} 


<6.0  n. mi les } 
<6.5  n. miles} 
<7.0  n. miles} 
<7.5  n. miles} 
<8.0  n. miles} 


POD  [DIST, OBSTIME] 
POR  [DIST, OBSTIME] 
END; 


{•*•«««««•««*««**«««*•«*«**«*«•«••««•**••*«*••«««««««*«««««*«**«»«««««««««*•»}. 
<»*«**«««*•*«*»**«•«**««•**«**«••*«««*«*••««*«« *'«**«**««****«4«*««*«««*««««««} 

{»***•«•••*****«•«***«»••«••••••««•«•*•*••*«••*«•«••**•*•*««*««««««««*««•«*•«> 

PROCEDURE  STOREDATA_IN_RAMFILE  (FILENAME  :  FILENAMETYPE ; 

LANTERNFILE  :  FILENAMETYPE; 

UAR  HEADER  :  HEADRECORD; 

UAR  READINGS  :  INTEGER); 

<reads  disk  data  input,  converts  buoy  versus  time  data  into  intensity 
versus  time  data  for  a  fixed  on  light.  FILENAME  is  data  location  on  disk. 
LANTERNFILE  is  specified  file  name  of  intensity  vs  angle  data  on  floppy  disk. 
If  data  is  actual  field  data  then  program  will  print  RMS  roll 
amplitude.  If  data  simulated  (with  ROLLSIM)  then  peak  roll  amplitude  and 
roll  period  are  printed  out.  Actual  vs  simulated  is  determined  by  whether  or 
not  HEADER . ROLLPERIOO  is  nonzero.  Returns  number  of  READINGS  stored  on  ram 
disk.  Multiplies  number  of  data  points  by  RATEMULTIPLIER  for  greater 
resolution  when  determining  effective  intensities.} 

UAR 

C  :  CHAR; 

T,  I  :  INTEGER; 

NUM  :  REAL; 

LANTERNDATA  :  LANTARRAYi 

DATASIMULATED  :  BOOLEAN;  (distinguish  between  actual  and  simulated} 
WAUEINCLUDED  :  BOOLEAN;  (true  if  wave  rider  data  included} 

RMSROLL  :  REAL;  (root  mean  square  roll  amplitude,  for  actual  data} 

MEANROLL  :  REAL; 
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TEMP  :  ARRAY  [ 0 . . RATEMULTIPLIER ]  OF  REALs  ^temporary  array  of  buoy  data) 
AN6LEINC  :  REAL; 

F  :  FILE  OF  REAL;  {ram  disk  of  intensity  vs  tirne> 

PRIN  :  TE.XTj  {printer) 

DSK  :  TEXT;  {floppy  disk  of  input  data) 

STR  :  STRING  140];  {variable  to  read  DSK) 

{*♦♦*♦♦******♦**♦♦♦*«**•»»**■»♦«*»****#♦*****♦«*♦■»**♦#*♦*****■»**********♦♦»♦*♦#} 

PROCEDURE  GET_LANTERN_DATA  (UAR  LANTERNDATA  :  LANTARRAY ; 

OAR  LANTERNFILE  :  FILENAMETYPE  ) ; 

{reads  intensity  versus  angle  data  from  lantern  file  stored  on  floppy  disk 
specified  in  CONST  block.  Data  used  to  convert  angle  vs  time  data  on  silo  to 
intensity  vs  time  which  is  stored  on  ram  disk.  If  different  lantern  file  is 
needed,  must  change  file  name  in  main  CONST  block.) 

CAR 

F  :  FILE  OF  REAL; 

I  :  INTEGER; 

BEGIN 

RESET  (F,  LANTERNFILE); 

LANTERNDATA  [01.  ANGLE  :=  -95.0;  {assures  that  up  to  90-deg  rolls  can) 

LANTERNDATA  1 0 1 . INTENSITY  :=  0.0;  {be  converted  to  intensity  values) 

I  :  =  1  ; 

UHILE  NOT  E0F<F)  DO 

WITH  LANTERNDATA  C I ]  DO 
BEGIN 

READ  (F,  ANGLE.  INTENSITY); 

I  :=I  +  1  ; 

END; 

CLOSE  (F); 

WHILE  I<=NUMANGLES  DO 

WITH  LANTERNDATA  [ 11  DO 
BEGIN 

ANGLE  9B; 

INTENSITY  :=  0; 

I  :=I  +  1  ; 

END; 

END; 

{*»»**»•»•»»»»••»*•••»»*»«•»»•««*•«••«•«•««***«*•*•••••««««««*«««***•«*•«•*•*) 

FUNCTION  ANSLE_TO„INTENSITY  (OAR  LANTERNDATA  :  LANTARRAY; 

BUOYANGLE  :  REAL)  :  REAL; 

{converts  buoy  angle  to  lantern  intensity.  max  angle  ever  seen  by  this 
function  is  +  or  -  90  degrees  due  to  data  acquisition  algorithm) 

UAR 

I  ;  INTEGER; 

I I  ,  AI  :  REAL; 

BEGIN 

I  :-l  ; 

WHILE  BUOYANGLE  >  LANTERNDATA  III. ANGLE  DO  I:-I+l; 

WITH  LANTERNDATA  [I-ll  DO 
BEGIN 

II  :=  INTENSITY; 

Al  :=  ANGLE; 

END; 

WITH  LANTERNDATA  [ 1 1  DO 

ANGLE_TO_INTENSITY  :  =»  1 1 +(  INTENS I TY- 1 1  )*(  BU0YANGLE-A1  )/(ANGLE-At  ); 

END; 

{♦*♦**♦****♦***♦**#**♦**♦***♦*#*♦*******»♦**♦**♦**♦#♦♦•****•****♦*♦*♦****♦♦♦*) 


BEGIN  <STOREDATA_IN_RAMFILE} 

SET_LANTERN_DATA  < LANTERNDATA  .  LANTERNFILE  )  i  <get  intensity  vs  angle  data) 
{reset  devices) 

RESET  (DSK  .  FILENAME  )} 

REURITE  (PRIN,  ’PRINTER:’)) 

OPEN  (F,  RAMFILE))  {open  ran  file  for  writing  to) 

CLOSE  <F,  ’PURGE’);  {purges  old  file  if  it  exists) 

REWRITE  (F.  RAMFILE);  {open  ran  file  for  writing  to) 

{read  file  header) 

READLN  (DSK,  STR); 

WRITELN  (PRIN.  STR); 

T:  =  l  ; 

REPEAT 

STRREAO  (STR.T.T.C); 

UNTIL  C=’  =  ’  ; 

STRREAO  (STR. T.T. HEADER. MONTH); 

STRREAD  (STR. T+1  .T, HEADER. DAY); 

STRREAD  (STR. T+1  ,T, HEADER. YEAR); 

WRITELN  (HEADER. MONTH.’/’ , HEADER. DAY  .’/ ’  .HEADER . YEAR  ) ; 

READLN  (DSK.  STR); 

WRITELN  (PRIN.  STR); 

T:-1  ; 

REPEAT 

STRREAD  (STR.T.T.C); 

UNTIL  C=’-’ ; 

STRREAD  ( STR  .T  ,T  . HEADER. STARTTIME  >; 

READLN  (DSK.  STR); 

WRITELN  (PRIN,  STR); 

T:  =  1  ; 

REPEAT 

STRREAO  (STR.T.T.C); 

UNTIL  C"=’  =  ’; 

STRREAD  (STR  ,T  .T  .HEADER . ONDURATION > ; 

WRITELN(’ON  DURATION  -  ’ .HEADER . ONDURATION  DIU  G000  :0,  ’  ninutes’); 
WRITELN; 

READLN  (DSK.  STR); 

WRITELN  (PRIN,  STR); 

T:  =  1  ; 

REPEAT 

STRREAD  (STR.T.T.C); 

UNTIL  C=’=’ ; 

STRREAD  ( STR  . T  . T  .HEADER . SAMPLERATE ) ; 

IF  HEADER. SAMPLERATE  >  3  THEN 
WAUEINCLUDED  FALSE 

ELSE 

WAUEINCLUDED  :■  TRUE;  {3  Hz  or  less  includes  wave  rider) 

READLN  (DSK  .  STR); 

WRITELN  (PRIN.  STR); 

T:«1  ; 

REPEAT 

STRREAO  (STR.T.T.C); 

UNTIL  C=’  =  ’  ; 

STRREAO  (STR  ,T  ,T  .HEADER .ROLLAMPLITUOE  ) ; 

READLN  (DSK.  STR); 

WRITELN  (PRIN,  STR); 

T:-l  ; 


REPEAT 

STRREAD  (STR.T.T.C); 

UNTIL  C=’  =  ’  i 

STRREAD  (STR  .T  .T  .HEADER. ROLLPERIOD  ); 

IF  HEADER. ROLLPERIOD  <>  0  THEN 
DATA5IMULATED  :=  TRUE 

ELSE 

DATASIMULATED  :=  FALSE) 

READLN  <DSK  ,  STR  ) ; 

URITELN  (PRIN,  STR)i 
PAGE  (PRIN)i 

{read  data  and  store  on  disk> 

READLN  (DSK,  NUM  ) ; 

WRITE  <F,  ANGLE_TO_INTENSITY  < LANTERNDATA .  NUM)  >; 

TEMP  [RATEMULTIPLIERI  :=  NUM; 

READINGS  :=  t; 

IF  NOT  DATASIMULATED  THEN 
BEGIN 

RMSROLL  NUM  *  NUM; 

MEANROLL  :=  ABS  ( NUM  ) ; 

END 

ELSE 

BEGIN 

RMSROLL  :»  0; 

MEANROLL  :=  0; 

END; 

IF  WAUEINCLUDED  THEN  {throw  away  wave  data} 

BEGIN 

READ  (DSK  .  C  >; 

READLN  (DSK,  NUM); 

END; 

WHILE  NOT  E0F(DSK  )  DO 
BEGIN 

READLN  (DSK  ,  NUM ) ; 

READINGS  READINGS  +  RATEMULTIPLIER ; 

IF  NOT  DATASIMULATED  THEN 
BEGIN 

RMSROLL  :•  RMSROLL  +  NUM»NUN;  {sun  the  squares} 

MEANROLL  :=  MEANROLL  +  ABS(  NUM  ); 

END; 

TEMP  [0]  TEMP  [RATEMULTIPLIER]; 

TEMP  [RATEMULTIPLIER]  NUM; 

ANGLEINC  (TEMP  [RATEMULTIPLIER]  -  TEMP  [0])  /  RATEMULTIPLIER 

FOR  I  1  TO  RATEMULTIPLIER  -  1  DO 

TEMP  [I]  I  •  ANGLEINC  +  TEMP  [0]; 

FOR  I  1  TO  RATEMULTIPLIER  DO 
BEGIN 

NUM  :=  ANGLE_TO_INTENSITY  (LANTERNDATA.  TEMP  [I]); 

WRITE  (F,  NUM); 

END; 

IF  WAUEINCLUDED  THEN  {throw  away  wave  data} 

BEGIN 

READ  (DSK  ,  C ); 

READLN  (DSK.  NUM); 


MM 


UD 


CLOSE  (DSK)i 
CLOSE  (PRIN); 

CLOSE  <F ,  ’SAOE’  ); 

IF  NOT  DATASIMULATED  THEN 
BEGIN 

RMSROLL  :=  SQRT(  RMSROLL  /  (READINGS  DIO  RATEMULTIPLIER  +  I  )  ); 
HEADER. ROLLAMPLITUDE  RMSROLLi  <R0LLPERI0D  still  zero  for  actual) 
HEADER. MEANROLLANGLE  :=  MEANROLL  /  (READINGS  DIO  RATEMULTIPLIER  +  1); 
ENDi 

URITELN( ’FINISHED  READING  DATA  FROM  DISK’); 


PROCEDURE  GET_DATA_WINDOW  (OAR  WIN  :  DATAWINDOWSt  (yindow  of  inton/tine) 

FILEPOS  :  INTEGER >1  {start  position  in  file) 
{Reads  a  yindow  of  data  from  ran  file  beginning  at  a  position  in  file 
indicated  by  FILEPOS.  Length  of  window  is  nax  observation  tine.) 

OAR 

F  :  FILE  OF  REAL; 

INTENOIF  :  REAL;  {inten  diff  between  two  actual  sanple  points) 
NUMPOINTS  !  INTEGER;  {nunber  of  data  points  to  read  into  window) 

I.  IL  :  INTEGER;  {indices) 

DI  :  0 .. RATEMULTIPLIER ;  {index  difference) 


INTENOIF  :  REAL;  {inten  diff  between  two  actual  sanple  points) 
NUMPOINTS  !  INTEGER;  {nunber  of  data  points  to  read  into  window) 

I.  IL  :  INTEGER;  {indices) 

DI  :  0 .. RATEMULTIPLIER ;  {index  difference) 

BEGIN 

OPEN  (F,  RAMFILE);  {point  to  start  of  data  record) 

NUMPOINTS  ROUND  (OBSTIMEINC  •  NUMTIMES  *  WIN . SAMPLERATE  ) ; 

READDIR  (F,  FILEPOS,  WIN. INTENSITY  [0]);  {positions  to  FILEPOS  and  reads) 


I  :  *  0 ; 

WHILE  KNUMPOINTS  DO 
BEGIN 

READ  (F,  WIN. INTENSITY  [I]); 
I  I  +  1 ; 

END; 

WHILE  HWINOOWLENGTH  DO 
BEGIN 

I  I  +  1  ; 

WIN. INTENSITY  C I ]  0; 

END; 

CLOSE  (F>; 


PROCEDURE  FLASH  PROPERTIES  (UAR  FLASH  :  FLASHARRAY; 


PROCEDURE  FLASH_PROPERTIES  (UAR  FLASH  :  FLASHARRAY; 

UAR  PERLEN  :  INTEGER; 

FLASHCHAR  :  CHARACTERISTICS; 

RATE  :  REAL  ) ; 

{Inputs  a  nornalized  flash  period  in  FLASH  to  be  used  for  finding  flashed 
intensity  vs  tine.  Sanplerate  in  FLASH  is  sane  as  WINDOW . INTENSITY ,  ie. 
buoy  roll  sanplerate  tines  ratenul t ip  1 ier .  Returns  nunber  of  points  per 
flash  period  in  PERLEN.) 

UAR 

I,  J.  TURNOFF,  TURNON  :  INTEGER; 
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BEGIN 

FLASH  [0]  :»  0i 
I  :=  1  i 

CASE  FLASHCHAR  OF 
FLB  :  BEGIN 

TURNOFF  :=  ROUND  (  0.B  *  RATE  )  +  1; 
WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  in  :=  1  ; 

I  :=  I  +  1  ; 

END) 

PERLEN  :=  ROUND! G  *  RATE); 

END; 

FL4  :  BEGIN 

TURNOFF  :=■  ROUND  (  0.4  *  RATE  )  +  1; 
WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  in  :=  t  ; 

I  :=  I  +  1 ; 

END; 

PERLEN  ;=  ROUND! 4  *  RATE); 

END; 

FL25  :  BEGIN 

TURNOFF  :=  ROUND  !  0.3  *  RATE  )  +  1; 
WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  C  n  :>  1 ; 

I  I  +  1  ; 

END; 

PERLEN  ROUND! 2. 5  ♦  RATE); 

END; 

QKFL  :  BEGIN 

TURNOFF  ;=  ROUND  !  0.3  ♦  RATE  )  +  1; 
WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  in  1  ; 

I  :  =  I  +  1  ; 

END; 

PERLEN  ROUND! RATE); 

END: 

GPFL5  :  BEGIN 

TURNOFF  ROUND  !0.4  *  RATE)  +  1; 
WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  in  t  ; 

I  I  +  1 ; 

END: 

TURNON  :=  ROUND  !RATE)  +  1; 

WHILE  I  <  TURNON  DO 
BEGIN 

FLASH  [ I ]  0; 

I  :=  I  +  1 ; 

END; 

TURNOFF  :=  ROUND  !1.4  •  RATE)  +  1; 
WHILE  I  <  TURNOFF  DO 
BEGIN 
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FLASH  [I]  :=  1 ; 

I  :=  I  +  1  i 
END; 

PERLEN  :=  R0UND<5  *  RATE); 

END; 

BEGIN 

TURNOFF  :=  ROUND  (RATE)  +  1; 

WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  [ I ]  :=  1 ; 

I  I  +  1  ; 

END: 

TURNON  :=  ROUND  (2  •  RATE)  +  1; 

WHILE  I  <  TURNON  DO 
BEGIN 

FLASH  in  :=  0; 

I  :*  I  +  1  ; 

END; 

TURNOFF  :=  ROUND  (3  *  RATE)  +  1; 

WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  [1]  :=  1 ; 

I  :«  I  +  1 ; 

END; 

PERLEN  :=  ROUND<G  *  RATE); 

END; 

BEGIN 

FOR  J  I  TO  6  DO 
BEGIN 

TURNOFF  ROUND(  (J  -  1  +  0.3)*RATE)  +  1; 
WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  [ I ]  :*  1  ; 

I  :»  I  +  1 ; 

END; 

TURNON  ROUND(  J  *  RATE  )  +  I; 

WHILE  I  <  TURNON  DO 
BEGIN 

FLASH  [I]  :«  0; 

I  :=  I  +  1 , 

END; 

END; 

PERLEN  ROUND(10  *  RATE); 

END; 

BEGIN 

TURNOFF  ROUND  (0.4  *  RATE)  +  1; 

WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  [I]  1 ; 

I  I  +  1  ; 

END; 

TURNON  ROUND  (RATE)  +  1; 

WHILE  I  <  TURNON  DO 
BEGIN 

FLASH  in  0; 

I  I  +  1  ; 
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ENDi 

TURNOFF  :=  ROUND  (3  *  RATE)  +  1; 
WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  II]  1 ; 

I  : «  I  +  1  5 
ENDi 

PERLEN  :=  ROUND! 8  *  RATE); 

END; 

EINTG  :  BEGIN 

TURNOFF  :=  ROUND  (  3.0  *  RATE  )  +  1; 
WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  Cl]  :»  1 5 
I  I  +  1  ; 

ENDi 

PERLEN  ROUND! S.0  *  RATE); 

END; 

0CC4  :  BEGIN 

TURNOFF  ROUND  !  3.0  *  RATE  )  +  ti 
WHILE  I  <  TURNOFF  DO 
BEGIN 

FLASH  in  1  ; 

I  I  +  1  i 

END; 

PERLEN  :=  ROUND! 4.0  *  RATE); 

END; 

END;  {CASE> 

WHILE  I  <=  100  *  RATEMULTIPLIER  DO 
BEGIN 

FLASH  in  0; 

I  I  +  1  i 
END; 


PROCEDURE  GET_FLASHED_INTENSITIES  ! OAR  WINDOW  :  DATAWINDOWS; 

UAR  FLASH  :  FLASHARRAY; 

PERLEN  :  INTEGER; 

K  :  INTEGER); 

{Converts  fixed  on  intensity  vs  tine  to  flashed  intensity  vs  tine  in  WINDOW. 
PERLEN  13  nunber  of  data  points  per  flash  period  in 

WINDOW  array.  Beginning  of  flash  period  is  indicated  by  K  upon  calling  this 
procedure.  Snail  nunber  added  to  intensity  assures  nonzero  during 
flash  pulse  but  does  not  noticably  change  effective  intensity.  This  assists 
locating  pulses  in  window  to  deternine  their  effective  intensities.! 

UAR 

I  :  INTEGER; 

BEGIN 

WINDOW. INTENSITY  [0]  0; 

I  :«  1  ; 

WITH  WINDOW  DO 

WHILE  I  <-  WINDOWLENGTH  DO 
BEGIN 

INTENSITY  III  ! INTENSITY  III  +  0.0000001)  *  FLASH  IK]; 
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I  :=  I  +  1  ; 

K  :=  ( K  +  1 )  MOD  PERLEN; 
END! 


ENDi 

•(*#♦***♦*********«< 

{*■*♦***♦**■*♦< 

{♦*♦*< 

PROCEDURE  GET_EUENT_EFI  < UAR  UINDOU 

FLASHLEN 


h  *  ♦  *  > 

)**«***«**«*«*«***«***«••*•«••*««*««} 
»•***«*««««*««•«•***«**«*«»«} 

DATAUINDOUSi 
INTEGER  )! 


{Determines  maximum  effective  intensity  for  all  observation  events. 

Stores  value  in  WINDOW. MAXEFI  [event],  the  EFI  for  detecting  that  event. 

Also  determines  minimum  EFI  of  all  flash  pulses  in  this  event  and 
stores  value  in  WINDOW. MINEFI  [event].  This  is  the  maximum  EFI  this  event 
produces  such  that  two  complete  flash  periods  are  recognized,  ie.  no  missing 
pulses  in  flash  period.  Uses  Schnidt-Clausen  method  to  calculate  effective 
intensity  <calls  EFFECTIUE_INTENSITY ) .  MINEFI  register  in  WINDOW  must  be 
reinitialized  because  flash  durations  longer  than  OBSTIMEINC  will  not 
record  an  efi  in  first  event . > 


UAR 

EUENT,  EUENTINC,  POINTS,  I,  K  :  INTEGER; 

EFI  :  REAL i  {equivalent  fixed  intensity,  effective  intensity! 

TIMEINC  :  REALi  {number  of  seconds  between  data  points  in  INTENT 
<*•«•*****»•**«•*«««*•««««*•«*«••«•*»«•»•«•**««•««*««**««•••««•«*•****«*«•»««} 
FUNCTION  EFFECTIUE_INTENSITY  (UAR  WINDOW  :  DATAWINDOWS; 

UAR  I  :  INTEGER; 

TIMEINC  !  REAL)  :  REAL; 

{INTEN  [I]  IS  first  nonzero  intensity  point  in  flash,  INTEN  [I-l]  is  zero. 

Integrates  until  INTEN  [I]  is  zero  again  and  STOP  is  reached.  This  function 
should  not  be  entered  unless  at  least  one  data  point  is  guaranteed  to  be 
nonzero,  otherwise  a  divide  zero  error  will  occur.  Uses  trapezoidal  formula 
to  estimate  definite  integral.! 

UAR 

AREA  :  REAL;  {integral  of  intensity! 

PEAK  :  REAL;  {peak  intensity  value  in  flash! 

BEGIN 

AREA  :=  0! 

PEAK  0; 

WITH  WINDOW  DO 

WHILE  INTENSITY  [ I ]  >  0  DO 
BEGIN 

AREA  AREA  +  INTENSITY  [I]; 

IF  INTENSITY  [I]  >  PEAK  THEN 
PEAK  INTENSITY  III; 

I  :=  I  +  1 1 
END; 

AREA  :=  AREA  *  TIMEINC; 

EFFECTIUE_INTENSITY  :=  AREA  /  <0.2  +  AREA/PEAK); 

END; 

{»*♦***♦♦*♦*♦*#*♦*♦***♦♦«•♦***«***#**********•************♦♦*****#****••*****} 

BEGIN 

WITH  WINDOW  DO 
BEGIN 

TIMEINC  1 /SAMPLERATE ;  {time  between  sample  points! 

POINTS  R0UND(  OBSTIMEINC  *  SAMPLERATE  );  {points  per  event  time  increm.! 
FOR  EUENT  : »  1  TO  NUMTIMES  DO 

BEGIN  {initialize  event  efi  registers! 


MINEFI  [EUENT]  :=  10000000-. 

MAXEFI  [EUENT]  :=  0; 

END  I 

EUENT  :«  1 ; 

I  :=■  0! 

WHILE  EUENT  <  NUMTIMES  DO  {analyze  all  pulses  in  flash  characteristic) 

BEGIN 

WHILE  INTENSITY  C I  1  =  0  DO 

I  :=  I  +  Ij  {move  to  first  non  zero  position) 

EFI  :=  EFFECTIUE_INTENSITY  (WINDOW.  I.  TIMEINC): 

EUENT  :=  (I  DIU  POINTS)  +  It  {determine  event  cell  number) 

IF  EUENT  <=  NUMTIMES  THEN 
BEGIN 

IF  EFI  >  MAXEFI  [EUENT!  THEN 

MAXEFI  [EUENT!  EFIt  {for  detecting  a  flash  in  increment) 

IF  EFI  <  MINEFI  [EUENT!  THEN 


MINEFI  [EUENT!  :=  EFIt  {detecting  all  flashes) 

END; 

END;  {WHILE  EUENT) 

{determine  detection  efi  vs  observation  time) 

FOR  EUENT  :=  2  TO  NUMTIMES  DO  {slide  maximum  up  to  right) 

IF  MAXEFI  [EUENT-n  >  MAXEFI  [EUENT!  THEN 
MAXEFI  [EUENT!  :=  MAXEFI  [EUENT-I!; 

{determine  recognition  efi  vs  observation  time) 
EUENTINC  !”  (2  *  FLASHLEN )  DIU  POINTSt  {difference:  end  -  start) 
IF  EUENTINC  >  0  THEN  {slide  two  period  window  down  to  left) 
BEGIN 

FOR  EUENT  :=  NUMTIMES  DOWNTO  1  DO 
BEGIN 

I  EUENT  -  EUENTINC  +  1 t 
IF  I  >  0  THEN 

BEGIN  {got  min  efi  in  this  2  period  window) 

FOR  K  :=  I  TO  EUENT  DO 

IF  MINEFI  [K!  <  MINEFI  [EUENT!  THEN 
MINEFI  [EUENT!  MINEFI  [K!i 

END 

ELSE 

MINEFI  [EUENT!  :=  0; 

ENDi 

END; 

FOR  EUENT  2  TO  NUMTIMES  DO  {slide  maximum  up  to  right) 

IF  MINEFI  [EUENT-l!  >  MINEFI  [EUENT!  THEN 
MINEFI  [EUENT!  MINEFI  lEUENT-H; 


END;  {WITH  WINDOW) 


PROCEDURE  INCREMENT_PROBS  (UAR  PROB  :  PROBRECORD; 

UAR  WIN  :  DATAWINDOWS  >t 

{Increments  the  POD  and  POR  arrays  in  PROB  depending  on  the  effective 
intensity  values  in  WIN. MAXEFI  and  WIN. MINEFI  arrays.  POD  is  the  probability 
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of  detecting  at  least  part  of  a  flash  cycle  given  mariner  at  specified 
distance  looks  in  direction  of  buoy  for  a  specified  length  of  time.  POR 
is  defined  as  the  probability  a  mariner  sees  two  complete  and  consecutive 
flash  cycles.  The  minimum  EFI  in  every  pair  of  adjacent  periods  was  found  in 
GET_E'JENT_EF I .  The  minimum  gives  distance  that  both  those  periods  are  seen. 
The  maximum  of  these  minima  was  then  selected.  This  gives  the  distance  that 
at  least  one  pair  of  periods  in  OBSTIME  time  increments  will  be  seen.  It  was 
assumed  that  if  one  can  detect /recognize  a  flash  at  distance  X  in  T  time,  then 
one  can  also  de t ec t /recogni ze  same  flash  in  time  greater  than  T.  Same 
flash  will  also  be  detected/recognized  at  shorter  distances,  but  may  require 
less  flash  periods  at  these  distances.! 

WAR 

PREVIOUS  :  REAL;  <efi  in  previous  period,  used  for  recognition! 

EFIDETEC  :  REAL;  {efi  for  detection  probability! 

EFIRECOG  :  REAL;  -Cefi  for  recognition  probability! 

EFITEST  :  REAL;  tminimum  efi  of  current  two  periods! 

DIST  :  INTEGER;  <distance  number! 

OBS  :  INTEGER;  {observation  time  numbers! 

P  :  INTEGER;  {period  counter! 

BEGIN 
PREUIOUS 
EFIDETEC  := 

FOR  OBS  :=  1 

BEGIN 

DIST  :  = 


UIN.MINEFI  [  I  ]; 
0; 

TO  NUMTIMES  00 


WHILE  (WIN.MAXEFI  [OBSI  >-  PROB . EFI C DIST ]  )  AND  (DIST  <-  NUMDISTS)  DO 
BEGIN 

PROB. POD  [DIST,  OBSI  :=■  1.0  +  PROB. POD  [DIST,  OBSI; 

DIST  :=  DIST  +  1 ; 

END; 

DIST  :=  1 ; 

WHILE  (WIN.MINEFI  COBS]  >=  PROB.EFICDIST] )  AND  (DIST  <-  NUMDISTS)  DO 
BEGIN 

PROB. POR  CDIST,  OBSI  1.0+  PROB. POR  [DIST,  OBSI; 

DIST  :=  DIST  +  1 ; 

END; 

END; 


I 

i 

4 


{**♦*♦**♦*•»*♦****♦♦******♦*♦*•■»*#*♦#****♦*♦♦♦♦*****♦**♦♦*******♦*********♦**«} 

{♦♦*****♦*♦***♦***«**♦*•»*•*«**#**»••*♦*****«*****♦♦#************♦*•*•»****♦***> 

PROCEDURE  N0RMALIZE_PR0BS  ( UAR  PROB  :  PROBRECORD; 

NUMWINDOWS  :  INTEGER); 

{normalizes  all  the  probabilities  by  the  number  of  windows  that  were  examined 
in  the  data  record! 

UAR 

DIST,  T  :  INTEGER; 

BEGIN 

WITH  PROB  DO 

FOR  DIST:-1  TO  NUMDISTS  DO 
BEGIN 

FOR  T:-1  TO  NUMTIMES  DO 
BEGIN 

POD  [DIST,  T]  :=  POD  CDIST,  T]  *  100  /  NUMWINDOWS; 

POR  CDIST,  T]  POR  CDIST.  T]  *  100  /  NUMWINDOWS; 

END; 
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ENDi 


END! 

<(4444444444444444444444444444#44444444444444444444444444444444444444444444444> 

•{4444444444444444444444444444444444444«44444444444444444444444444444444444444} 

PROCEDURE  PRINT_PROBABILITIES  ( UAR  OUTDISK  :  PROBFILEi 

PROB  :  PROBRECORDi 
HEADER  :  HEADRECORD: 

FLASHCHAR  :  CHARACTERISTICS)! 

lUlrites  Probabilities  to  floppy  disk  and  printer.  If  first  Prob  record  then 
procedure  determines  disk  file  name  from  header  data  and  opens  disk  file  for 
writing.  If  last  Prob  record  it  closes  disk.  Roll  period  is  printed  for 
simulated  data.  Simulated  vs  actual  is  determined  by  whether  or  not 
rol lamp  1 1 1 ude  equals  zero.) 

UAR 


P  :  1  .  .2; 

I,  D,  T  :  INTEGER! 

FILENAME  :  STRINGC20]! 

S  ;  STRING  1801! 

BEGIN 

IF  FLASHCHAR  =  FL6  THEN  {get  file  name  and  open  disk  file) 
BEGIN 

FILENAME  :=  ’#4:D’! 

STRWRITE  ( FILENAME  .5 .T .HEADER . MONTH : 0  )  i 
FILENAME  :=  FILENAME  +  i 
STRWRITE  ( FILENAME, T+1  .T .HEADER.DAY:0  )  i 
FILENAME  :=  FILENAME  +  ’H’j 

STRWRITE  ( filename, T+)  ,T, HEADER. STARTTIME  DIU  360000  :0)| 
REWRITE  (OUTDISK,  FILENAME)! 

ENDi 

WRITE  (OUTDISK  ,  PROB  )  i 

IF  FLASHCHAR  =  0CC4  THEN  {close  if  last  prob  record) 

CLOSE  (OUTDISK  ,  ’SAUE’  >! 

FOR  P:=1  TO  2  DO 
BEGIN 

WRITELN(  PRIN  ,  ESC  ,  ’ 8rk3S  ’  ) !  {expanded-compressed  10.7  cpi) 
WRITELN(PRIN)! 

WRITELN(PRIN)i 
IF  P=I  THEN 

WRITELN( PRIN , '  ’ : 22  , ' Probab 1 1 1 1 y  of  Detection') 

ELSE 


WRITELNCPRIN,’ 
CASE  FLASHCHAR  OF 


FLG 

WRITELN 

FL4 

WRITELN 

FL25 

WRITELN 

QKFL 

WRITELN 

GPFL5 

WRITELN 

GPFL6 

WRITELN 

IQKFL 

WRITELN 

MOA 

WRITELN 

EINTG 

WRITELN 

0CC4 

WRITELN 

END!  {CASE) 

WRITELN  (PRIN)! 
WRITELN  (PRIN)! 


:21  .’Probability  of  Recognition'  )i 


(PRIN  .  ’ 

’  :29. 

’FL. 

6  (0 

.6  ) 

’  )! 

(PRIN,’ 

’  :29, 

’FL. 

4  (0 

.4  ) 

’  )i 

(PRIN,’ 

’  :28  . 

’FL. 

2.5 

(0. 

3  )’  ) 

(PRIN,’ 

’  :31  , 

’QK. 

FL.  ’ 

)l 

(PRIN,’ 

’  :30, 

’GP. 

FL. 

5’  ) 

1 

(PRIN,’ 

’  :30. 

'GP. 

FL. 

G’  ) 

y 

(PRIN,’ 

’  :30  , 

’I. 

QK.  FL.  ’ 

)i 

(PRIN,’ 

’:3'  . 

'MO. 

( A  )’ 

)i 

(PRIN,’ 

’  :30  . 

’E. 

Int  . 

6’  ) 

y 

(PRIN,’ 

’=32. 

’  Occ 

.  4’  ) 

f 
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URITELN  (PRIN)! 

URITELN  (PRIN,  ESC,'&k0S’);  {normal  print  12cpi> 

WRITELN  ( PRIN  Di stance ’  ’ :25  , ’Observation  Tine  (s)’); 

WRITE  (PRIN,’  ’;l  ,’Nt.Mi.  ’  ); 

WRITE  (PRIN,  ESC,’&k2S’);  tconpressed  print  21.3  cpi) 

FOR  I:  =  1  TO  NUMTIMES  DO  WRITE  ( PRIN  ,  ( I *OBSTIMEINC) : S  ,  ’  ’:3)i 
WRITELN  (PRIN); 

WITH  PROB  DO 

FOR  Dr-I  TO  NUMDISTS  DO 
BE6IN 

WRITE  (PRIN,  ESC,’8,k0S’  );  Inornal  print) 

WRITE  (PRIN,’  ’ :2  ,(D  *  DISTANCEINC):3: 1  ,’  ’:3); 

WRITE  (PRIN,  ESC,’&k2S’);  {compressed) 

WRITE  (PRIN,’  ’;6); 

IF  P-1  THEN 
BEGIN 

FOR  I:-1  TO  NUMTIMES  DO  WRITE  (PRIN, POD  ID,!]  :7:0,’X’) 

END 

ELSE 

BEGIN 

FOR  I:-1  TO  NUMTIMES  DO  WRITE  (PRIN,POR  10,11  :7:0,’%’) 
END; 

WRITELN(PRIN); 

WRITELN(PRIN); 

END; 

WRITE  (PRIN,  ESC.’&k0S’);  {normal  print) 

FOR  I:-1  TO  4  DO  WRITELN  (PRIN); 

WITH  HEADER  DO 
BEGIN 
S  :  =  ’  ’  ; 

FOR  I:-1  TO  80  DO  STRWRITE  (S,I,T,’  ’); 

STRWRITE  (S,5,T,’Date:  ’  ); 

STRWRITE  (S,1 1  ,T,M0NTH:2); 

STRWRITE  (S,13.T.’/’  ); 

STRWRITE  ( S  ,14  ,T  ,DAY:2  ); 

STRWRITE  (S,16,T.’/’  ); 

STRWRITE  (S  ,17  ,T  ,YEAR:2 >; 

STRWRITE  (S, 45, T. ’Time:  ’  ); 

STRWRITE  ( S, 51  .T. STARTTIME  DIU  360000  :2); 

STRWRITE  (S,53,T,’:’  ); 

STRWRITE  ( S, 54  ,T.( STARTTIME  MOD  360000)  DIU  6000  :2); 

STRWRITE  (S,56  ,T,’ : ’  ); 

STRWRITE  ( S, 57, T,(( STARTTIME  MOD  360000)  MOD  6000)  DIU  100  :2) 
WRITELN  (PRIN,  S); 

FOR  I:-1  TO  80  DO  STRWRITE  (S.I.T,’  ’); 

STRWRITE  (S. 5, T, ’Sample  Rate:  ’); 

STRWRITE  ( S  ,  1 8  , T  .SAMPLERATE : 2 : 2  ) ; 

STRWRITE  (S,21  .T.’Hz’  ); 

STRWRITE  (S, 45, T, ’Lantern  File:  ’); 

STRWRITE  (S,59,T,LANTERNFILE); 

WRITELN  (PRIN,  S); 

FOR  I:-1  TO  80  DO  STRWRITE  (S.I.T,’  ’); 

IF  ROLLPERIOOO0  THEN  {data  is  simulated) 

BEGIN 

STRWRITE  (S,5,T,’Roll  Amplitude:  ’); 

STRWRITE  (S,21  ,T .ROLLAMPLITUDE: 2 :0  ) ; 


p 


I 


1 


y 


STRWRITE  (S,24.T,’  deg.’ ); 

STRUIRITE  (  S  ,45  ,T  Sinusoidal  Roll  Period:  ’  )i 
STRWRITE  (S  ,G9  ,T  .R0LLPERI0D:2: 1  ); 

STRWRITE  VS.73.T.’  sec’ ); 

END 

ELSE 

BEGIN 

STRWRITE  <S,5.T.’RMS  Roll  Amplitude:  ’); 

STRWRITE  <S.25.T.R0LLAMPLITUDE:2:l  >! 

STRWRITE  <S.29,T.’  deg.’  ); 

STRWRITE  (S  .45. T, ’Mean  roll  angle:  ’  ) ? 

STRWRITE  ( S  .62  ,T  .MEANR0LLAN6LE:2: 1  ); 

STRWRITE  (S.67.T.’deg. *  ); 

END; 

WRITELN  (PRIN  .  S ); 

END; 

IF  FLASHCHAR  =  FLG  THEN 
BEGIN 

FOR  I:-l  TO  80  DO  STRWRITE  (S.I.T,’ 

STRWRITE  <S. 5. T. 'Output  Data  File:  ’ >; 

STRWRITE  (S. 23, T, FILENAME); 

WRITELN  (PRIN.  S); 

END; 

PAGE  (PRIN); 

END; 

END; 

{»*»»•«*•»*«««»••••*•**«»*««••»**«•«•««»*•«•****«»**««•««»•«««•*«•*•«•«**«»•*} 

{f*************************************************************'**************} 

<*««*««»»*«***«»»»»»»»»»•*»»»»•»»»•*•»»«*»»»•»»*»«»•«»»«•«»»»••««**««*««««««»} 

BEGIN  <BU0Y_PR08ABILITIES> 

REWRITE  (PRIN.  ’PRINTER:’  ); 

INPUT_USER_OATA  (FILENAMES.  LANTERNFILE  ) : 

I  :=  1;  {counter  for  data  blocks  on  disk! 

WHILE  FILENAMES  III  <>  ’NONE’  DO  {repeat  for  all  data  sets) 

BEGIN 

STOREDATA_IN_RAMFILE  (FILENAMES  111.  LANTERNFILE,  HEADER,  NUMREADINGS ) ; 
RANGE  :=  (NUMREADINGS  -  WINDOWLENGTH )  DIU  2;  {for  random  pos> 

WINDOW. SAMPLERATE  :»  HEADER . SAMPLERATE  *  RATEMULTIPLIER ; 

FOR  FLASHCHAR  :-  FLG  TO  0CC4  DO  {repeat  for  all  characteristics) 

BEGIN 

SEED  :*■  1;  {seed  for  random  number  generator) 

INITPR0B_REC0RD  (PROBABILITY);  {zero  prob  array) 

FLASH_PROPERTIES  (FLASH.  FLASHLEN  ,  FLASHCHAR,  WINDOW . SAMPLERATE  ) ; 

FOR  W  :-  1  TO  NUMWINDOWS  DO 
BEGIN 

FILEPOS  :-  2  *  RAND  (SEED.  RANGE)  +  1; 

GET_DATA_WINDOW  (WINDOW,  FILEPOS); 

WITH  WINDOW  DO  {zero  maxefi  and  minefi  arrays) 

FOR  TIME  :-  1  TO  NUMTIMES  DO 
BEGIN 

MAXEFI  [TIME]  :=  0; 

MINEFI  [TIME]  :»  0; 

END; 

STARTPOS  :-  RAND  (SEED,  FLASHLEN); 

GET_FLASHEO_INTENSITIES  (WINDOW.  FLASH.  FLASHLEN,  STARTPOS); 
GET_EUENT_EFI  (WINDOW,  FLASHLEN); 


S.S.T.’RMS  Roll  Amplitude:  ' 
S  .25 .T.ROLLAMPLITUDE:2:  1  >; 
S.2g,T.’  deg.’  ); 

S. 45, T, ’Mean  roll  angle:  ’)i 
S .62 ,T .MEANROLLANGLE:2: 1  ); 
S.67,T.’deg.  ’  ); 


DO  STRWRITE  (S.I  .T,’  ’  ); 
T, 'Output  Data  File:  ’  ); 
,T .FILENAME  ); 

S); 


d 
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{PLOT_WINDOW  (WINDOW)?}  {for  testing) 

INCREMENT_PROBS  (PROBABILITY.  WINDOW); 

END; 

NORMALIZE_PROBS  (PROBABILITY,  NUMWINDOWS); 

PRINT_PROBABILITIES  (OUTDISK,  PROBABILITY,  HEADER.  FLASHCHAR  ) ; 
END;  {for  flashchar) 

I  I  +  1  ; 

END;  {while  datalocat ion) 

CLOSE  (PRIN); 


Program  LanternProfileMaker 


PROGRAM  LanternProfileMaker 
{Generates  lantern  bean  profiles 
vertical  divergence  ( FWHM  ) , 
lanternfile  used  in  PROGRAM 
CONST 

lanternfile  "  ’t3:LANT.REAL’ i 
DegToRad  •  0.017453293; 

TYPE 

filenameType  •  STRING  [201; 

UAR 

P  :  TEXT; 

filenane  :  f i lenaneType; 
flux  :  REAL; 
fullWidth  :  REAL; 
linitLevel  :  REAL; 
k  :  REAL; 
integral  :  REAL; 
peak  :  REAL; 
s  :  STRING  M0]; 
i  :  INTEGER; 


a  user  specified 
an  existing 
Brown,  4/2/87> 


{nane  of  lantern  file  to  integrate) 
{total  flux  in  integrated  profile) 

{full  width  half  nax  in  degrees) 
{fraction  of  peak  at  linits  of  integ) 
{constant  for  gaussian  exponential) 
{integral  of  gaussian) 

{peak  of  gaussian) 


(INPUT.  OUTPUT >; 

with  a  Gaussian  shape, 
with  the  sane  total  flux  as 
BUOY_PROBABILITIES.  by  Dan 


. . . . . . 

. . . 

PROCEDURE  IntegrateLanternFi le  (filenane  :  f i lenaneType ; 

linLevel  :  REAL; 

UAR  flux  :  REAL); 

TYPE 


lantRecord  ■  RECORD 

angle  :  REAL; 
intensity  ;  REAL; 

END; 

lantArray  ■  ARRAY  C0..4S5]  OF  lantRecord; 


) 

) 


UAR 

lanternOata  :  lantArray; 
naxintensity  :  REAL; 
anglel  nc  :  REAL; 
f  :  FILE  OF  REAL; 
i  :  INTEGER; 

. . . . * . * . .*.*.) 

PROCEDURE  GetLanternData  (UAR  lanternData  :  lantArray; 

UAR  naxlnten  :  REAL); 

{reads  lantern  file  on  disk) 

UAR 

i  :  integer; 

BEGIN 

lanternData  C0]. angle  -95.0; 
lanternData  C0] . intensity  :■  0.0; 
naxlnten  :■  0.0; 
i  :  -  1  ; 

WHILE  (NOT  EOF(f))  AND  ( i  <  455)  DO 
WITH  lanternOata  C 1 1  00 
BEGIN 

READ  (f,  angle,  intensity); 

IF  intensity  >  naxlnten  THEN 
naxlnten  intensity; 

:■  i  +  t  ; 


i 


i 


i 


u 


END; 

WHILE  i  <=  455  DO 
BEGIN 

lanternData  [il.a.Tgle  95; 
lanternData  C  i  ] .  int ens i t y  :=  0.0; 

1  :  =  1  +  1  ; 

END; 

END; 

PROCEDURE  PnntLanternData  <  UAR  lanternData  :  lantArray); 

UAR 

p  :  TEXT; 

1  :  INTEGER; 

BEGIN 
i  0; 

REWRITE  (p  ,  'PRINTER: ’  ); 

WHILE  lanternData  li  I. angle  <  90  DO 
BEGIN 

1  :  =  1  +  1  ; 

WRITELN  <p.  lanternData  Cil. angle  :5:1,  lanternData  [ 1 1 . int ens i ty : 1 0 : 2 ) ; 
END; 

END; 

BEGIN  < IntegrateLanternPi le> 

RESET  (f,  filenapne); 

GetLanternData  (lanternData,  naxlntonsi ty  ) ; 

CLOSE  <  f  ); 
i  :  =  )  ; 

WHILE  lanternData  Cil. angle  <  0  00  i  :»  i  +  I; 
linLevel  :  •=  limLevel  ♦  max  Intensi  ty ; 

WHILE  lanternData  C  i  I . intensity  >  linLevel  DO  i  : *  i  -  1 ; 

angleinc  :•  lanternData  Ci+ll. angle  -  lanternData  til. angle;  tin  deg.} 

WITH  lanternData  Cil  DO 

flux  intensity  ♦  cos  (angle  *  DegToRad); 
i  :  =  1  +  1  ; 

WITH  lanternData  Ii]  DO 

flux  flux  +  4  *  intensity  *  cos  (angle  *  DegToRad); 
i  :  =  1  +  I  ; 

WHILE  lanternData! 1 1 . intensity  >  linLevel  DO 
BEGIN 

WITH  lanternData  [il  DO 

flux  :  =■  flux  +  2  *  intensity  •  cos  (angle  *  DegToRad); 

1  :  =  i  +  1  ; 

WITH  lanternData  Cil  DO 

flux  :■  flux  +  4  *  intensity  *  cos  (angle  *  DegToRad); 
i  :  •  i  +  1  ; 

END; 

WITH  lanternData  ! i 1  DO 

flux  flux  +  intensity  *  cos  (angle  *  DegToRad); 
flux  :=  angleinc  *  DegToRad  *  flux  /  3; 

END; 

FUNCTION  Gauss  (  k  ,  x  :  REAL)  :  REAL; 

BEGIN 

Gauss  :«  exp  (  k  *  x  *  x  ); 
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i 

I 


a 


{****-»-»**'4******4**-»***1Ht*****i 

{**44**44*4**4*****4M*4*4*4***^ 


t4**4****44****1 


^^^i^**4***^^*44*****^^*****4***4*4****^^} 


PROCEDURE  IntegrateGaussian  (  hwht*i ,  {half  width  half  max) 

lim  :  REAL;  {%  peak  at  limit  of  integrate) 

UAR  k  :  REAL;  (const  for  exponential) 

UAR  SUM  :  REAL); 

{finds  total  flux  in  a  gaussian  beam  profile) 

CONST 

nMax  =  400; 

UAR 

n  :  INTEGER; 
middle  :  INTEGER; 
angle  :  REAL; 

angleLimit  :  REAL;  {limits  of  integration  of  gaussian) 

anglelnc  :  REAL; 

BEGIN  { IntegrateGaussian) 
k  :=  LN  (0.5)  /  hwhm  /  hwhm  ; 
angleLimit  :=  SQRT(  LN(  lim  )  /  k  ); 
anglelnc  :=  2  »  angleLimit  /  nMax; 
angle  -angleLimit; 

sum  :=  Gauss<k,  angle)  *  cos  (angle  *  DegToRad); 
angle  :=  angle  +  anglelnc; 

sum  :=  sum  +  4  •  Gauss(k,  angle)  *  cos  (angle  *  DegToRad); 
middle  :»  nMax  DIU  2; 


WHILE  n  <  nMax  DO 
BEGIN 

angle  :■  (n  -  middle)  *  anglelnc; 

sum  :=  sum  +  2  •  Gau33(  k,  angle  >  *  cos  (angle  *  DegToRad); 
n  :  =  n  +  1  ; 

angle  :=  (n  -  middle)  *  anglelnc; 

sum  :=  sum  +  4  »  Gauss(  k,  angle  )  *  cos  (angle  *  DegToRad); 
n  :  =  n  +  )  ; 

END; 

angle  :=  (n  -  middle)  •  anglelnc; 

sum  :=  sum  +  Gau35(  k,  angle  )  *  cos  (angle  *  DegToRad); 
sum  :=  sum  ♦  anglelnc  *  DegToRad  /  3; 


{*****♦»♦**♦*****■»#******»***♦•*»****#***♦♦**♦♦< 

■{*•**♦♦♦•**•♦**♦«•♦*♦******«■»■»*******•»*♦****#*#*< 

PROCEDURE  StoreGaussProf i le  (k,  peak,  fullWidth 
CONST 

anglelnc  =  0.2;  (degrees) 

UAR 

filename  :  f i lenameType ; 
f  ;  FILE  OF  REAL; 
angle  :  REAL; 

1  :  INTEGER; 

BEGIN 

filename  ’#3:GAU’; 

5TRWRITE  (filename,  7,  i,  R0UND( ful IWidth  )  :0  ); 
filename  filename  +  ’.REAL'; 

WRITELN  (’File  name  is:  ’.filename); 

REWRITE  (f,  filename); 

1  ROUND  (-45  /  anglelnc  ); 


(**«*}. 


REAL); 
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'I*,. 


while  angle  45.0  DO 
BEGIN 

angle  :=  i  ♦  angleinci 

WRITE  <f,  angle,  peak  *  6aus5< k  .angle )  )i 

1  :  =  1  +  1  i 

ENOi 

CLOSE  ( f  .  ’SAUE’  >! 

END; 

{**»**««««*«**««*****•«***«*«»•••*«»««»*»««»«»•*****••**«*«*****«««*•*•««»**} 

BEGIN  {MAIN) 

REWRITE  (P  .  'PRINTER: ’  ); 

WRITELN  (’Enter  file  name  (including  volume  name)  of  lanternf i le ’  ) ; 

WRITELN  ('to  match  output  flu.-^  to.’); 

READLN  (filename); 

WRITELN  (  ’Enter  fraction  of  peak  intensity  at  limits  of  integration’  ); 

WRITELN  (  ’default  is  0.01  ’  ); 

1 imi tLevel  : =  0.01; 
s  :  =  ’  ’  ; 

READLN  ( 5  ); 

IF  3  <>  ’’  THEN 

STRREAD  ( s  ,  1  , i , 1 imi tLevel ) ; 

WRITELN  ( 1 imi tLevel  ) ; 

WRITELN  (’Enter  full  vertical  width  (in  degrees)  of  generated  profile’); 
WRITELN  (’at  50%  intensity’); 

READLN  ( fullWidth  ); 

I ntegrateLanternFi le  (filename,  limitLevel,  flux); 

WRITELN  (  ’ flux  "  ’  ,flux  ); 

IntegrateGauss ian  (fuiiWidth  /  2,  limitLevel,  k,  integral); 
peak  :=  flux  /  integral; 

WRITELN  (’PEAK  INTENSITY  OF  GAUSSIAN  =  ’.peak); 

WRITELN  {  ’k  =  ’  .k  ); 

StoreGaussProf  i  le  (V.,  peak,  fullWidth); 

END. 


Program  FastFourierTransform 

PROGRAM  FaatFourierTranaformC INPUT  .OUTPUT  )i 

{Performs  a  Fast  Fourier  Transform  on  any  number  of  data  points  equal  RAOIX 
to  power  of  GAMMA .  User  inputs  approx .  number  of  data  points.  Program 
calculates  global  variables  GAMMA  and  N  (number  of  data  points).  User  may 
convolve  data  with  a  five  number  kernel  prior  to  FFT  to  remove  noise, 
differentiate,  etc.  Five  numbers  (or  weights)  are  user  input.  Convolution 
15  a  running  average  of  current  point  and  two  points  on  either  side  of  curre-  t 
with  weighting  defined  In  kernel.  Windowing  of  time  domain  is  also  allowed 
after  done  convolving.  Four  windows  are  provided:  rectangular,  Hamming, 
Hanning,  and  Hanning  Squared.  User  may  also  take  average  of  several  FFT’s  t  . 
specifying  that  more  than  one  data  files  are  to  be  transformed. 

BY  DAN  BROWN,  April  1987) 

IMPORT 

lODECLARATIONS  . 

SYSDEUS  , 

DGL_LIB . 

DGL_INQ; 

CONST 

RADIX  =■  2;  {base  or  radix  of  FFT  routine) 

NMAX  -  1G384| 

PI  -  3. 14l592G54i 

TYPE 

COMPLEX_NUMBERS  -  RECORD 

RE  :  REAL; 

IM  :  TEAL; 

END; 

OATA_ARRAY  -  ARRAYI0. .NMAX ]  OF  COMPLEX_NUMBERS ; 

POINTER  *  •'DATA_ARRAY; 

FILESTRING  -  STRINGI201; 

PLOT_TYPES  -  (  MAGNITUDE.  PHASE.  TIMESERIES  ); 
messageType  »  STRING  CG0]; 

UAR 

GAMMA  :  INTEGER; 

N  :  INTEGER:  {number  of  data  points  is  RADIX  to  the  power  of  GAMMA) 

X  :  POINTER;  {points  to  data  array) 

HEAP  :  POINTER;  {used  to  restore  heap  memory) 

PTR  :  TEXT;  {file  name  for  printer) 

numfiles  :  INTEGER; 

FILENAME  '■  FILESTRING;  {file  name  for  data  file  on  disc) 

TIMEINC  :  REAL;  {time  in  seconds  between  data  samples) 

VARIANCE  :  REAL:  {mean  square  amplitude) 

MAXFREQ  :  REAL; 

IMAX  :  INTEGER; 

1  :  INTEGER; 

{•♦•**#*#**********»****#***********************«***»*J****************«***-; 

. . * . * . * . 

PROCEDURE  CLEARSCREEN; 

{clears  CRT  of  all  text  and  spaces  cursor) 

BEGIN 

WRITELN( CHR(  I  2  )  ) ;  {homes  cursor  and  clears  screen) 

WRITELN; 

WRITE  (  ’  '  ); 

END; 


■{*******♦**♦*♦*******»«*♦******»*********•**♦♦♦*♦****♦♦♦•**»***»*«•»«♦«♦♦♦»♦*) 

FUNCTION  AnsuierYes  (message  :  messageType)  :  BOOLEAN! 

{returns  TRUE  if  anwer  with  a  Y  (yes),  returns  FALSE  if  uith  N  (no)} 

UAR 

ans  :  CHAR; 
ansUalid  :  BOOLEAN! 

BEGIN 

CLEARSCREEN! 
ansUalid  :=  FALSE: 

REPEAT 

WRITE  (message): 

WRITE  (  '  ( Y  or  N)  ’  )! 

READ  ( ans  ) ; 

CLEARSCREEN! 

CASE  ans  OF 

’Y’  .’y’  .’N’  ,’n’  :  ansUalid  :=  TRUE; 

OTHERWISE  END!  {CASE} 

UNTIL  ansUalidi 
CASE  ans  OF 

’Y’.’y’  :  Answer  Yes  .'=•  TRUE; 

’N’.’n’  :  AnswerYes  :•  FALSE: 

OTHERWISE 

END: 

ENDi 

{*•»***«*♦********♦****•********•»*»•**♦*♦****•**♦*********♦********»*********} 

-(**«*««***««««««•«*«««•««*«««««*«»«••«*•«•*«*««•««*«*•*««•«*««••••««•«••«••«} 

PROCEDURE  WriteMessage  (message  :  messageType > ! 

UAR 

1  :  INTEGER: 

BEGIN 

CLEARSCREEN: 

FOR  i  :=  1  TO  10  DO  WRITELN  (’ 

WRITELN  (message): 

END: 

{»**••*»*••****••«•»*••••««««•««•*«•«•«••*«*««*•*»•«*••*«««««»*««*««•«*•****} 

{*•**•*««««««*««**«*««««•*««»•«««*»««««»«««***•«««*«•«««•«•«««««**«««»««««•«} 

FUNCTION  DIGIT_REUERSE(  Z ,  GAMMA  :  INTEGER)  :  INTEGER: 

{forms  the  bit  reversed  number  of  Z.  Example:  0010  ->  0100.} 

UAR 

I,  Y  :  INTEGER: 

BEGIN 
Y :  =0 : 

FOR  I:»t  TO  GAMMA  DO  {gamma  is  also  number  of  digits} 

BEGIN 

Y  :=  RADIX  *  Y  +  <Z  MOD  RADIX):  {shift  left,  add  rightmost  digit} 

Z  :=  Z  DIU  RADIX:  {shift  right} 

END: 

DIGIT_REUERSE:-Y: 

END: 

{*♦***♦***♦*♦*********•*«***»*******■»********•********♦*****♦******♦****♦***} 

{*****■*»**♦*********♦**♦♦*«****#*##♦*#**«*******#♦***♦#**«***#*•»•»**«•»##***#♦} 

PROCEDURE  LOAD_DATA_ARRAY<  UAR  X  :  POINTER: 

N  ;  INTEGER: 

UAR  TIMEINC  :  REAL: 

UAR  UARIANCE  :  REAL: 

FILENAME  :  FILESTRING  ): 


{reads  in  data  from  a  disc  file  and  stores  data  into  array  X.  Procedure 
assumes  all  data  in  file  is  real,  no  imaginary  part.  Thus  zeros  are  put 
in  imaginary  part  of  array.  Data  file  can  have  more  or  less  than  N  points.} 

S  :  5TRINGC80]i  {data  string  from  ASCII  file) 

F  :  TEXT) 

I ,  T  :  INTEGER) 

NUMLINES  ;  INTEGER) 

MEAN  :  REAL) 

BEGIN 

RESET(F. FILENAME); 

WRITE  (’number  of  lines  of  data  file  to  skip?  ’); 

READLN  (NUMLINES); 

FOR  I:=l  TO  NUMLINES  DO 
BEGIN 

READLN  (F.  S); 

WRITELN  (S); 

END) 

WRITELN; 

WRITELN) 

WRITE  (’sampling  rate  of  data  (Hz)?  ’); 

READLN  (TIMEINC); 

TIMEINC  :•  1 /TIMEINC;  {no.  seconds  between  points} 

MEAN  0) 

UARIANCE  0; 

FOR  I:-0  TO  N-1  DO 
BEGIN 

IF  NOT  EOF(F)  THEN 
BEGIN 

REAOLN(F.X''[I].RE)) 

UARIANCE  :»  UARIANCE  +  Xnn.RE  *  X-'CII.RE; 

MEAN  MEAN  +  X''Cn.RE) 

END 

ELSE 

X'’[I].RE  :»  0.0;  {for  files  with  less  than  N  points} 

X-'Cn.IM  :=  0.0; 

END; 

CL05E(F  )) 

CLEARSCREEN; 

WRITELN  (PTR,  'TOTAL  ENERGY  IN  TIME  DOMAIN  -  ’.UARIANCE); 

UARIANCE  •-  UARIANCE  /  N; 

WRITELN  (PTR,  ’UARIANCE  OF  WAUE  AMPLITUDE  -  ’, UARIANCE)) 

WRITELN  (PTR,  ’MEAN  OF  WAUE  AMPLITUDE  »  ’.MEAN  /  N)) 

END) 

{*•««•*•**««•*•••*•«**«••«•••»«••«»•«•»•••«*•*««*«**««**«»««««••«•••«*•»•«»«> 
{*#*#♦  **♦****#**♦**********♦*♦**♦****■»«*********«***♦*♦*♦#•»****  *******♦****♦> 

PROCEDURE  CONUOLUE  (UAR  X  ;  POINTER;  N  :  INTEGER); 

{Convolves  data  with  a  five  number  kernel} 

UAR 

kernel  :  ARRAY  C-2..+2]  OF  REAL; 
buffer  :  ARRAY  [-2..0]  OF  REAL; 
i  .  j  ,  k  :  INTEGER) 
sum  :  REAL) 

F  :  TEXT; 

BEGIN 


WRITELN  (’input  convolution  Kernel’); 

FOR  1 :=  -2  TO  2  DO 
BEGIN 

WRITE  ( 1  ->’  ); 

READLN  (kernel  Cil); 

END; 

WRITELN  (PTR  . ’CONUOLUTION  KERNEL:’); 

FOR  i:=  -2  TO  2  DO  WRITELN  (PTR,  kernel  Ill): 

FOR  k  :=  -2  to  0  DO 

buffer  Ik]  :=  0.0;  Iinit  buffer) 

WriteMessage  (’WAIT  '  CONOOLOIN6  DATA’); 

FOR  i:=  0  TO  N-3  DO 
BEGIN 
sun  :=  0; 
j  3; 

WHILE  (j  >  -2)  AND  ( i  +  j  >  0 )  DO 
BEGIN 

J  :=  J  -  1 ; 

sun  :=  sun  +  kernel  Cj]  *  X^Ii+jl-RE; 

END; 

IF  i-2  >  0  THEN 

X''[i-3].RE  :=  buffer  [-21;  Irenove  last  entry  fron  buffer) 

FOR  k:*=  -2  TO  -1  DO 

buffer  Ck]  buffer  Ck+ll;  <3lide  entries  down  one) 

buffer  [0]  :=  sun; 

END; 

FOR  k:— 2  TO  0  DO 

X''[N-3+k  I  ,RE  buffer  Ckl;  {empty  buffer) 

CLEARSCREEN; 

{REWRITE  <F,  ’#3:C0NU0L.TEXT’  >; 

FOR  k:=0  TO  N-l  DO 

WRITELN  (F,  X'CkI.RE); 

CLOSE  (F  .  ’SAVE’  ); ) 


END; 

FUNCTION  Signif icantOigits  (nun  .’REAL;  digits  :  INTEGER)  :  REAL; 
{returns  ''nun'"  rounded  to  "digits"  significant  digits) 

VAR 

comp  :  INTEGER; 

BEGIN 

comp  :=  10;  {integer  with  "digits"  significant  digits) 

WHILE  digits  >  1  DO 
BEGIN 

comp  :■  comp  ♦  10; 
digits  digits  -  1 ; 

END; 

digits  :*  0;  {number  of  digits  to  shift  num  after  rounding) 

WHILE  num  (  comp  DO 
BEGIN 

num  :■  num  ♦  10; 
digits  : ■  digits  -  1 ; 

END; 

WHILE  num  >  comp  DO 
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urn 


REAL; 

STRIN6C80] 

REAL; 

REAL; 

INTEGER; 


CHARWIOTH,  CHARHEIGHT 
TITLE 

X  ,  DX  .  Y  ,  DY 
XLABEL,  YLABEL 
I  ,  J  .  K 
BEGIN 

CHARWIDTH  :=  < XMAX  -  XMIN>  /  50; 
CHARHEIGHT  :=  ( YMAX  -  YfllN)  /  25; 
SET_CHAR_SIZE( CHARWIDTH .CHARHEIGHT ) ; 
CASE  PLT  OF 


<5ets  constant  relative  size) 


TIMESERIE5 

MAGNITUDE 


BEGIN 

SET_TEXT_ROT  (0.1  >; 

MOUE(XMAX .0 )i 
LINE(0.0>; 

X  :  =•  0 ; 

DX  (XMAX  -  XMIN  )  /  40;  {Label  :<  axis) 

FOR  I:=1  TO  4  DO  {major  tick  marks,  horizontal  axis) 

BEGIN 

FOR  J:=l  TO  10  DO  {minor  tick  marks) 

BEGIN 

X  :=  X  +  DX; 

M0UE{X .0  ); 

LINE(X  .-0.5*CHARHEIGHT ); 

END; 

LINE(X. -CHARHEIGHT); 

TITLE:"*’;  {label  major  ticks) 

IF  PLT  "  TIMESERIES  THEN 
BEGIN 

XLABEL  X  *  TIMEINC;  {time) 

XLABEL  Sign! f icantOigits  (XLABEL,  3); 

STRWRITE( TITLE, 1  .K  .XLABEL : 3 : 2 > ; 

END 

ELSE 

BEGIN 

XLABEL  ;=  X  /  TIMEINC  /  N;  {frequency) 

XLABEL  :=  Signif icantDigits  (XLABEL,  3); 

STRWRITE(TITLE.I  .K  .XLABEL : 4 : 2  )  i 

END; 


MOUE( X  .CHARHEIGHT ); 
GTEXT( TITLE  ); 

END; 

SET_TEXT_ROT  (  1  ,0); 
MOUE(0.YMAX ); 
LINE(0,YMIN); 
LINE(-CHARWIDTH,YMIN); 

Y  YMIN; 

DY  (YMAX  -  YMINl/IB; 
FOR  I:-1  TO  10  DO 
BEGIN 

Y  ;=  Y  +  DY; 
MOUE(0,Y); 
LINE(-CHARWIDTH,Y)| 
END; 


{vortical  axis  tic  marks) 


YLABEL  :=  S igni f icantOigi t s  ( YMAX  ,  2)\ 

STRWRITE( TITLE  .t  ,K  .YMAX ; 4:2  )  i 
MOUE( -STRLEN( TITLE  )*CHARUIDTH  ,  YMAX  )  i 
6TEXT(TITLE  ); 

ENDi 

PHASE  :  BEGIN 

SET_TEXT_ROT  <0.1); 

MOUE(XMAX  ,0  ); 

LINE(XMIN,0); 

MOLlE(0.YriIN); 

LINE(0,YMAX  ); 

X  XMIN; 

DX  (XMAX  -  XMIN  )  /  80;  {Label  x  axis) 

FOR  I:“l  TO  8  00  {major  tick  marks) 

BEGIN 

FOR  J:=l  TO  10  DO  {minor  tick  marks) 

BEGIN 

X  X  +  DX; 

MOUE<X  .0  ); 

LINE<X .-0.5*CHARHEIGHT  >; 

END; 

LINE(X.-CHARHEISHT); 

TITLE:=’';  {label  major  ticks) 

XLABEL  X  /  XMAX  /  TIMEINC  /  2;  {frequency) 

STRWRITE< TITLE. 1  .K  .XLABEL :0 : 2  ) ; 

MOOE<X .CHARHEI6HT); 

GTEXT< TITLE  ) 

END; 

SET_TEXT_ROT  ( 1  .0  ); 

Y  YMIN;  {label  y  axis) 

DY  PI/18; 

MOUE(-CHARUIDTH  .YMIN); 

LINE(CHARWIDTH.YMIN); 

TITLE  ’-90  deg’ ; 

6TEXT<  TITLE); 

FOR  I:-l  TO  18  DO 
BEGIN 

Y  Y  +  DY; 

M0yE(-CHARWIDTH/2.Y)i 

LINE(CHARWI0TH/2.Y); 

END; 

LINE(CHARUIDTH .YMAX  ); 

TITLE:-’+g0  deg’; 

GTEXT( TITLE); 

END; 

OTHERWISE 
END;  {CASE) 

CASE  PLT  OF 

TIMESERIES  :  TITLE  'Time  Series’; 

MAGNITUDE  :  TITLE  ’Magnitude’; 

PHASE  :  TITLE  'Phase  angle’; 

END;  {CASE) 

MOUE  <  (XMAX  -  XMIN  -  STRLENl TITLE  )*CHARWIDTH  )/3  +  XMIN.  YMAX  ); 
GTEXT( TITLE); 

END; 


PROCEDURE  SET_U5ER_SCALE  ( UAR  P  :  POINTER; 

PLOT  :  PLOT_TYPES; 

XMAX  :  INTEGER; 

TIMEINC  ;  REAL; 

UAR  MAXUALUE  :  REAL); 

UAR 

UALUE  ;  REAL; 

I  :  INTEGER; 
naxPoint  INTEGER; 

XMIN,  YMIN,  YMAX  :  REAL; 

ENERGY  :  REAL; 

BEGIN 

MAXUALUE  :=  0; 

ENERGY  :=  0; 

FOR  I:=0  TO  N- I  DO  {find  Maximum  value  in  data  arrayT 

BEGIN 

UALUE  P''[I].RE*P''Cr].RE  +  P"!  11 .  IM»P'' C 1 1 .  IM; 

ENERGY  ENERGY  +  UALUE; 

IF  UALUE>MAXUALUE  THEN 
BEGIN 

MAXUALUE: -UALUE; 
maxPoint  I; 

END; 

END; 

MAXUALUE  SQRT( MAXUALUE  ) ; 

WRITELN  (PTR.  ’MAX  UALUE  =  ’.MAXUALUE,’  AT  I-’ .maxPoint  )  i 
CASE  PLOT  OF 

TIMESERIES  :  BEGIN 

{WRITELN  (PTR.  ’TOTAL  ENERGY  IN  TIME  DOMAIN  -  ’.ENERGY);) 

XMIN:-0.0i 

YMIN:-  -MAXUALUE; 

YMAX:-  MAXUALUE; 

END; 

MAGNITUDE  :  BEGIN 

{WRITELN  (PTR, ’TOTAL  ENERGY  IN  FREQUENCY  DOMAIN  =  ’. 

ENERGY  /  N>;> 

XMIN:-  0.0; 

YMIN;-0.0; 

YMAX: -MAXUALUE; 

END; 

PHASE  :  BEGIN 

XMIN:— XMAX; 

YMIN:— PI/2; 

YMAX:-PI/2; 

END; 

END;  {CASE) 

SET_WINDOW( XMIN  .XMAX .YMIN .YMAX  );  {Sets  user  coord,  scale) 

LAB£L_PLOT  < PLOT  .XMIN , XMAX .YMIN, YMAX); 

END; 

{«««•««««•«««««*«««««««*«««««•«««««««••«««****«««****•*»»»**•••*»•*»**»••«»•} 

BEGIN  {PLOT_DATA) 

INITIALIZE_PLOT; 

IF  PLOT-TIMESERIES  THEN 

BEGIN  <.****#  time  SERIES  PLOT  *****) 

SET_USER_SCALE  (P .TIMESERIES  .IMAX  .TIMEINC  .MAXMAG  ); 

MOUE(0,  P-'C0].RE); 
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FOR  I:=1  TO  IMAX  DO 
LINEC  I  ,P"[  n.RE  )i 

IF  AnsuerYes  (’HARD  COPY  OF  PLOT*?  ’  )  THEN 
CALL{DUMP6RAPHICSH00K  >? 

CLEAR_DISPLAY; 

END 

ELSE 

BEGIN 

{*♦.**  MAGNITUDE  PLOT  *•***> 

SET_USER_SCALE  ( P  .MAGNITUDE  . IMAX  .TIMEINC  .MAXMAG  )  i 
MOUE(0,SQRT(  P-'Cai  .RE*P''[0].R£  +  P-C  0 1 .  IM*P'' [  0 1 .  IM  )  ) ; 

FDR  I:=1  TO  IMAX  DO 
BEGIN 

Y  :>  SQRT(  P'l  II  .RE*P''C  II  .RE  +  P*  [  II .  IM*P' C I  ] .  IM  >i 
LINEd  .Y)? 

END; 

IF  AnsuerYes  (’HARD  COPY  OF  PLOT?’)  THEN 
BEGIN 

CALL( DUMPGRAPHICSHOOK  ) ; 

PAGE  (PTR); 

END; 

CLEAR_DISPLAY; 

<##**»  phase  PLOT  *##*♦> 

SET_USER_SCALE  ( P .PHASE .IMAX .TIMEINC  .MAXMAG  ) ; 

HALF_N  N  DIU  2; 

I:“N  -  IMAX;  {max  negative  freq  starts  at  N/2  or  oreater> 

FIRSTPOINT  TRUE; 

REPEAT 

X  I  -  (I  DIU  HALF_N)*N;  <I-N  when  I>N/2.  I  when  I<N/2> 

IF  P''CI].IM*P''n].IM  +  P‘'n].RE*P''CI].RE  >  0.01*MAXMA6  THEN 
BEGIN  <non  zero  phase) 

IF  Pd  11. RE-0  THEN  <90  deg  phase) 

BEGIN 

IF  PdI].IN<0  THEN  Y;-  -PI/2 
ELSE  Y:-  PI/2; 

END 

ELSE  <  phase  <  90  ) 

BEGIN 

Y  ARCTAN<P''Cn.IM/P’’[I].RE); 

END; 

END 

ELSE  {zero  phase) 

BEGIN 

Y:-0; 

END; 

IF  FIRSTPOINT  THEN 
BEGIN 
MOUE(X  .Y  ); 

FIRSTPOINT  FALSE; 

END 

ELSE 

LINE(X  .Y); 

I  (1  +  1  )  MOD  N;  {start  midway  and  wrap  around) 

UNTIL  I-IMAX; 
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4. 


IF  ApswerYes  (’HARD  COPY  OF  PLOT?’)  THEN 
BEGIN 

CALL( DUMPGRAPHICSHOOK  >; 

PAGE  (PTR); 

END) 

CLEAR_DISPLAYi 
END-, 

GRAPHICS_TERN! 

ENDi 

{*«•»*»«*«••««*•••««*«««*••*•••**••••«»•••««**•«* 
■(•»*♦***♦*♦*♦*♦♦*♦♦****■»***♦♦*♦***********•♦♦***** 

PROCEDURE  ORDER_DATA_ARRAY(UAR  X  :  POINTER;  GAMMA 
{unscrambles  the  data  array) 

UAR 

I  :  INTEGER;  {current  node) 

J  :  INTEGER;  {node  to  be  exchanged) 

NODE  :  COMPLEX_NUMBERS; 

BEGIN 

FOR  I:=0  TO  N-1  DO 
BEGIN 

J  DI6IT_REUERSE( I .  GAMMA); 

IF  J<I  THEN  {exchange  only  those  nodes  not  yet  exchanged) 


BEGIN 

NODE 

-  x'-in 

X"!  I  ] 

:=  X'CJ 

X'lJ  ] 

:=  NODE 

END; 

END; 

END; 

{*•«•««»**««••***»«*•«•«•***«««•*•««•««**•**•*««*••*«*«***«*«***«««*«*«•••««} 

{«*«*•«•******»•*«*****«««*•*««•«•*••«•*****•*«•*«««***««««*««««««•***•*«•**} 

PROCEDURE  WINDOW_DATA  (UAR  X  :  POINTER;  N  :  INTEGER >; 

CONST 

RETURN  =  CHR<32); 

UP  =  CHR( 31  ); 

DOWN  -  CHR( 10  ); 
numSelects  =  4; 

UAR 

C  :  CHAR; 

I  ,  T  :  INTEGER; 

ITEM  :  0. .4; 

SELECTION  :  (RECTANGLE,  HAMMING,  HANNING,  SOHANNING); 

NUM  :  REAL; 

{•*«*4*«*»*«***«*******«*»*««»«*««*«»«««»*«**«««*««**»«*»*«««**«««***««««*«*} 

PROCEDURE  SETCURSOR(J  :  INTEGER); 

{positions  cursor  next  to  Jth  item  in  selection) 

UAR 

I  :  INTEGER; 

BEGIN 

URITE( CHR( 8  )  ) ;  {move  cursor  left  one  position  if  possible) 

WRITE( ’  ’  ); 

WRITE( CHR( 1  ) ) ;  {home  cursor) 

FOR  I:  =  1  TO  3  DO  WRITE( CHR( 1 0  ) ) ; 

FOR  I:-1  TO  2B  DO  WRITE( CHR( 28  )  > ; 

WHILE  JO0  DO 


*•******«••••«««**••*••*•••). 

,  N  :  INTEGER); 


{move  cursor  down) 
{move  cursor  to  right) 


BEGIN  {position  cursor  next  to  selection) 

! 

WRITE( CHR( I  0  )  ) ;  {move  cursor  down  one  line) 

END  I 

WRITE<’<’)i  {write  the  pointer) 

END; 

BEGIN 

CLEARSCREEN; 

WRITELN(  CHR< 1  )  ); 

WRITELN(  'Select  a  data  filter  window  by  setting  pointer  with  up  or  down’  ); 
WRITELN(  'arrow  keys.  Then  press  <ENTER>  key  to  make  selection'  ); 

WRITELN; 


WRITELN( ’* 

Rectangle 

*’  ); 

WRITELN( ’* 

Hanning 

*’  ); 

WRITELN(  ’* 

Hanning 

*’  ); 

WRITELN( ’* 

Hanning  Squared 

*’  ); 

ITEM  t; 

WRITELN; 

SETCURSOR( ITEM  >; 

REPEAT 

READ(C); 

CASE  C  OF 

RETURN  :  CASE  ITEM  OF 

1  :  SELECTION  RECTANGLE; 

2  :  SELECTION  HAMMING; 

3  :  SELECTION  HANNING; 

4  :  SELECTION  SQHANNING; 

END; 

UP  :  BEGIN 

ITEM  :*  < { ITEM-2+nunSelects  )  MOD  nunSelects >+l ; 

URITE( CHR( 1 0  ) ) I  {negate  key  press  with  down  arrow) 

END; 

DOWN  :  BEGIN 

ITEM  :*  (ITEM  MOD  nunSelects)  +  1; 

URITE( CHR< 31  ) ) ;  {negate  key  press  with  up  arrow) 

END; 

OTHERWISE 
END;  {CASE) 

SETCURSOR< ITEM  ) ;  {position  cursor) 

UNTIL  C-RETURN; 

CLEARSCREEN; 

CASE  SELECTION  OF 
HAMMING  :  BEGIN 

WRITELN  ('WAIT  !  WINDOWING  WITH  HAMMING'); 

FOR  T:-0  TO  N-t  00 

yiTI-RE  X'lTI.RE  *  ( 0 .08+0. 4B-0 . 46*C0S( 2*PI *T/N  ) ) ; 

END; 

HANNING  :  BEGIN 

WRITELN  ('WAIT  !  WINDOWING  WITH  HANNING’); 

FOR  T:>-0  TO  N-1  00 

X'CTI.RE  X-CTJ.RE  *  (0.5-0.5*COS(2*PI*T/N ) ); 

END; 

SQHANNING  :  BEGIN 

WRITELN  (’WAIT  !  WINDOWING  WITH  HANNING  SQUARED’); 

FOR  T:-0  TO  N-t  DO 


RECTANBLE  : 

END;  <CASE> 
CLEARSCREEN; 


BEGIN 

NUM  :=  0.5-0.5«COS(2*PI«T/N)( 
X“[T].RE  :=■  X-'Cn.RE  *  NUM  *  NUM; 
END; 

END; 

BEGIN  <no  change} 

END; 


PROCEDURE  FFT_RADIX_2  < UAR  X  :  POINTER;  GAMMA,  N 
{performs  radix  2  FFT  on  data  pointed  to  by  X> 
CONST 

PI  =  3.14t592BS4; 


INTEGER  ); 


m 


1 

KSi 


ARRAY_NUM  :  INTEGER;  {index  o 

K  ;  INTEGER;  {index  o 

K_DUAL  :  INTEGER;  {index  o 

DN0DE_SPACING  :  INTEGER;  {index  s 

FIRST_DNODE  :  INTEGER;  {beginni 

P  :  INTEGER;  {integer 

ANGLE  :  REAL;  {angle  a 

ANGLE_INC  :  REAL;  {angle  i 

COSINE  :  REAL;  {value  o 

SINE  :  REAL;  {value  o 

WX_REAL  :  REAL;  {real  pa 

WX_IMAG  :  REAL;  {imagine 

BEGIN 

ANGLE_INC  -2*PI/N; 

DN0DE_SPACING  N; 

FOR  ARRAY_NUM:-1  TO  GAMMA  00 
BEGIN 

DNODE_SPACING  :=  0N0DE_SPACING  DIV  2; 
K  :=  0; 

WHILE  K  <  N  DO 
BEGIN 

FIRST  ONODE  K  +  DNOOE  SPACING; 


{index  of  array} 

{index  of  data  node  in  array} 

{index  of  the  dual  node  of  node  K} 

{index  spacing  between  dual  nodes  in  array} 
{beginning  index  of  next  set  of  dual  nodes} 
{integer  power  of  complex  exponential} 

{angle  argument  in  complex  exponetial} 

{angle  increment} 

{value  of  the  cosine  of  the  angle} 

{value  of  the  sine  of  the  angle} 

{real  part  product  of  node  and  complex  exp.} 
{imaginary  part  product  of  node  and  exp.} 


{reduced  by  two  for  each  array} 
{start  at  top  of  array} 


FIRST_0N0DE  K  +  DN00E_SPACING ; 

{determine  power  to  W} 

P  DIGIT_REUERSE<  K  DIU  DN0DE_SPACING 
ANGLE  P  *  AN6LE_INC; 

COSINE  :=  C0S(  ANGLE  ); 

SINE  SIN(  ANGLE  ); 

WHILE  K  <  FIRST_DN0DE  DO 
BEGIN 

K_DUAL  K  +  DN0DE_SPACING ; 

WX  REAL  X'lK  DUAL]. RE  *  COSINE  - 


GAMMA  >; 


WX_REAL  X " I K_DUAL ] . RE  *  COSINE  -  X'-[K_DUAL] .  IM  *  SINE: 
WX_IMAG  :»  X''CK_DUAL].IM  *  COSINE  +  X " C K_DUAL ] . RE  *  SINE; 
X''[K_DUAL3.RE  :=  X'CKl.RE  -  WX_REAL ;  {calculate  dual  nodes} 
X-[K_DUAL].IM  :=•  X^IKl-IM  -  WX_IMAG; 

X-IKl.RE  X-'CKl.RE  +  WX_REAL;  {calculate  K  nodes} 

X'CKl.IM  X''[K].IM  +  WX_IMAG; 

K  K  +  1  ; 

END: 

■  K  +  DN0DE_SPAC ING ;  {skip  dual  nodes,  they're  done  already} 
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END; 


w 

k- 

m  * 

I 


> 


} 


END; 

END; 

. . . . . . . . . . 

{•*•»••**•««*•••*««•••*•«•«•••••«««••«*••*«•••******•**««••••••»••*••«*•••*«} 

PROCEDURE  Derivative  ( UAR  X  ;  POINTER;  TIMEINC  :  REAL;  N  ;  INTEGER): 
{multiplies  FFT  by  imaginary  radian  frequency.  Done  when  transforming 
derivative  of  a  function  and  not  taking  derivative  until  after  transform) 

UAR 

j  :  INTEGER; 

cnumber  :  COMPLEX_NUMBERS ; 
freq  :  REAL; 

BEGIN 

FOR  j  :=  0  TO  N  DIU  2  DO 
BEGIN 

freq  :=  j  *  2  *  PI  /  TIMEINC  /  N; 
cnumber  : =  X "  [ j  ] ; 


X ' [ J ] . RE  : =  freq  • 

cnumber . IM ; 

X  ■'  [  J  ]  .  IM  :  =  -  freq 

*  cnumber. RE; 

END; 

FOR  J  (N  DIV  2 )  +  1 

TO  N-1  DO 

BEGIN 

freq  : «  ( j  -  N )  *  2 

*  PI  /  TIMEINC 

cnumber  : =  X " C  j  ] ; 

X  ^  C  J  1 . RE  : »  freq  * 

cnumber . IM; 

X'l J 1, IM  -  freq 

*  cnumber. RE; 

END; 

END; 

. . 

{' 

PROCEDURE  PouerSpectralDensity  ( UAR  X  :  POINTER; 

TIMEINC  :  REAL; 

N  :  INTEGER; 

VARIANCE  :  REAL); 

{places  magnitude  of  FFT  squared  in  real  part  of  complex  array) 

VAR 

1  :  INTEGER; 
energy  :  REAL: 
norm  :  REAL: 

BEGIN 

energy  :»  0; 

FOR  i  0  TO  N-1  DO 
BEGIN 

X'Cil.RE  X'lil.RE  *  X'Cil.RE  +  X-'IiI.IM  •  X'til.IM; 

X'Cil.IM  01 

energy  :■  energy  +  X^Cil-REi 
END: 

norm  :■  VARIANCE  ♦  N  *  TIMEINC  /  energy; 

FOR  1  0  TO  N-?  00 

X^Iil-RE  :■  norm  ♦  X^Cil-REi 

END; 

. . . . . . . . ****** . . 

. . * . ***************** . ********** . . 

PROCEDURE  AverageFFTs  (VAR  X  :  POINTER;  N,  iteration,  numfilas  :  INTEGER); 
{iteration  is  one  if  first  FFT  and  indicates  to  store  values  on  RAM  file 
without  adding) 


FT  :  FILE  OF  COMPLEX_NUMBERS i 
j  :  INTE6ER! 

cnurnber  :  COMPLEX_NUMBERS ; 

BE6IN 

IF  iteration  =  I  THEN 
BEGIN 

REWRITE  (FT,  ' RAM : AUETRANS ’  )  f 
FOR  j  :=  0  to  N-1  DO 
WRITE  ( FT  ,  X"! J  ]  >i 

END 

ELSE 

BEGIN 

RESET  (FT.  ’RAM:AUETRANS’ ); 

FOR  J  :=  0  TO  N-l  DO 
BEGIN 

READ  (FT,  cnumber); 

X'ljI.RE  :=  X''[j].RE  +  cnunber .  RE  i 
y'ljI.IM  ;=  X''Cj].IM  +  c  number .  IM ! 

ENDt 

CLOSE  (FT)i 

REWRITE  (FT,  ’ RAM : AUETRANS ’  ) ; 

FOR  J  :•  0  TO  N-l  DO 
WRITE  (FT.  X-Cj]  )i 

ENDi 

CLOSE  (FT  .  ’SAVE’  )  i 

IF  iteration  *  nunfiles  THEN 
BEGIN 

FOR  j  :=  0  TO  N-l  DO 
BEGIN 

X-CjI-RE  :=  X-IjI-RE  /  numfiles; 

X"[J].IM  :=  X^IjI-IM  /  nunfilesi 
ENDi 

ENDi 

ENDi 

{**♦**♦*♦■*♦**♦♦*♦********♦*♦******•»***•***♦***♦*♦♦•»♦*■»*****♦*♦«♦***♦♦*♦*♦♦***). 
{**«******«*«*•««**••**••«**«•«•«*»»•««««•»**««««•«***«**««*« 

BEGIN  {FFT} 

REWRITE(PTR. ’PRINTER:’  ); 

MARK(HEAP)i  <nark  top  of  heap  to  later  restore  it  to  rnen.I 
NEW(X  )i 

WRITE  (’number  of  data  files  to  transform  and  average?  ’>; 

READLN  ( numf i los  )  i 

WRITE  ( ’number  of  data  points  to  transform?  ’  )j 
READLN  (N)i 

GAMMA  ;=  ROUND!  LNIN)  /  LN( RADIX)  )t 
N  R0UND(  EXP  ( GAMMA*LN< RADIX  ) )  )i 
IF  N  >  NMAX  THEN 
N  :=  NMAXi 


FOR  i  :•  1  TO  numfiles  DO 
BEGIN 

CLEARSCREENi 

WRITE! ’DATA  FILE  NAME  IS:  ’)( 


READLN( FILENAME )  I 

WRITELN  (PTR,  ’FOURIER  TRANSFORM  OF  DATA  FILE:  ’.FILENAME); 

REPEAT 

LOAD_DATA_ARRAY< X .  N.  TIMEINC,  UARIANCE,  FILENAME); 

PLOT_DATA  (X.  N-1  ,  TIMEINC,  N.  TIMESERIES); 

IF  AnswerYes  (’Convolve  data")”  )  THEN 
BEGIN 

CONUOLUE  ( X .  N ); 

PLOT_OATA  (X.  N-1.  TIMEINC.  N,  TIMESERIES); 

END; 

IF  AnswerYes  (’Window  data?’)  THEN 
BEGIN 

WINDOW_DATA  (X,  N); 

PLOT_DATA  (X,  N-1,  TIMEINC.  N,  TIMESERIES); 

END; 

UNTIL  AnswerYes  (’Finished  with  tine  domain?’); 

WriteMessage  (’WAIT  !  CALCULATING  THE  FFT’); 

FFT_RADIX_2  (X,  GAMMA.  N); 

ORDER_DATA_ARRAY  (X,  GAMMA.  N)i 
CLEARSCREENi 

{if  determining  transfer  function  from  step  response  data) 

IF  AnswerYes  (’transform  step  response  to  impulse  response?’)  THEN 
Derivative  (X,  TIMEINC.  N); 

<if  plotting  power  spectral  density,  square  FFT> 

IF  AnswerYes  ('plot  power  spectral  density')  THEN 
PowerSpectralDensity  <X.  TIMEINC,  N.  UARIANCE); 

IF  numfiles  >  1  THEN  {add  new  to  old  FFT> 

AverageFFTs  (X,  N,  i,  numfiles); 

REPEAT 

WRITE  (’max  frequency  to  plot?  ’); 

READLN  (MAXFREQ); 

CLEARSCREEN; 

IMAX  ROUND  (MAXFREQ  *  TIMEINC  *  N); 

IF  IMAX  >  N  DIU  2  THEN 
IMAX  N  OIU  2; 

PLOT.DATA  (X.  IMAX.  TIMEINC.  N.  MAGNITUDE); 

UNTIL  AnswerYes  (’Is  plot  OK?’); 

END; 

DISPOSE  (X); 

RELEASE( HEAP  ) ;  (restores  heap  to  original  state  to  use  again) 


